TRANSLATION FROM RUSSIAN TO ENGLISH THE BOOK “BLAST EFFECTS 
CAUSED BY EXPLOSIONS” AUTHORED BY B. GELFAND AND M. SILNIKOV 


Final Technical Report 
by 


Gelfand Boris 
(April 2004) 


United States Army 
EUROPEAN RESEARCH OFFICE OF THE U.S. ARMY 
London, England 
CONTRACT NUMBER N62558-04-M-0004 


Professor Boris Gelfand 


Approved for Public Release; distribution unlimited 


. Form Approved 
Report Documentation Page OMB No, 0704-0188 


Public reporting burden for the collection of information is estimated to average | hour per response, including the time for reviewing instructions, searching existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, 
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington 
VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penaity for failing to comply with a collection of information if it 
does not display a currently valid OMB control number. 


1. REPORT DATE 2. REPORT TYPE 3. DATES COVERED 


00 APR 2004 N/A - 


4, TITLE AND SUBTITLE 5a. CONTRACT NUMBER 
Translation from Russian to English the Book "Blast Effects Caused by | 5. Ga ane ntipER 

Explosions" Authored by B. Gelfand and M. Silnikov, Final Technical | y69558 _94-M-0004 
Report 


5c. PROGRAM ELEMENT NUMBER 


6. AUTHOR(S) 5d. PROJECT NUMBER 


Se. TASK NUMBER 


5£. WORK UNIT NUMBER 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT 
European Research Office of the U.S. Army, London, England NUMBES 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR’S ACRONYM(S) 


11. SPONSOR/MONITOR’S REPORT 
NUMBER(S) 


12, DISTRIBUTION/AVAILABILITY STATEMENT 
Approved for public release, distribution unlimited 


13. SUPPLEMENTARY NOTES 
The original document contains color images. 


14. ABSTRACT 


15. SUBJECT TERMS 


16. SECURITY CLASSIFICATION OF: 17. LIMITATION 18. NUMBER | 19a. NAME OF 
OF ABSTRACT OF PAGES RESPONSIBLE PERSON 
a. REPORT b. ABSTRACT c. THIS PAGE UU 249 


unclassified unclassified unclassified 


Standard Form 298 (Rev. 8-98) 
Prescribed by ANSI Std Z39-18 


Abstract 


This report contains the translation from Russian to English of the book “Blast effects 
caused by explosions” authored by B. Gelfand and M. Silnikov. 
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INTRODUCTION 


The expected development of criminal situation and the analysis of terrorist activity 
over the recent years give rise to the necessity of in-depth study of modern methods of 
application, search, classification, and neutralization of explosive devices (ED). Particular 
attention must be paid to the prevention of the non-authorized operation and the reduction 


of destructive consequences of ED [1... 16]. 


The reason of ED becoming widespread is the availability of the information about 
the technology of ED manufacturing and application, as well as the data on initial materials 
and accessories (blocks) for their manufacture in home-craft conditions. For example, 
originally published in 1971 (and repeatedly republished) in USA the book “The Anarchist 
Cooking Book” contains the comprehensive technological data and instructions for 
manufacturing of the simple ED in a kitchen! Criminalistic analysis showed, that 
significant part of high explosives (HE) and ED used by the terrorists and criminals were 
really produced at home conditions, and necessary ingredients were freely b ought in the 


popular shops of chemicals. 


The following facts can serve as the evidence of wide use of mine-explosive 
devices in terroristic acts. According to the report of the Ministry of Justice more than 
3000 cases of illegal application of ED were registered in 1995 in USA. In Great Britain 
after covenant with the Irish republican army (IRA) on the discontinuance of fire in Ulster 
the amount of terroristic acts with the use of ED was not reduced, but even has increased 


approximately by 15 percents. 


The number of criminal explosions in Russia sharply increases after the 


disintegration of USSR. 


The statistical data show, that ED with small and average energy of explosion (up 
to 1 kg of TNT equivalent) are mainly used in Western Europe and Russia; in USA - more 
powerful ED (from 3 up to 5 kg of TNT equivalent). In recent years the terrorists more 


often use the explosive devices with large weight of HE. 


The most known terroristic act of this category is the undermining of office 
buildings in Oklahoma city (USA, 1997), in which ED containing more than 400 kgs of 
TNT equivalent charge was used. The number of victims can illustrate the scale of this 


crime: 168 men were died and more than 500 were heavy injured. 


The tendency of increasing the ED power, used by terrorists in Russia, is distinctly 
traced since 1995. Suffice it to say about the undermining of apartment houses of the 
families of Russian frontier guards in Caspyisk, explosion of railway station in Pyatigorsk, 
terroristic act in the house belonged to the mayor of Makhachkala, where 17 men were 
died and 26 men were heavy wounded. 30 buildings were destroyed and damaged. ED 
power exceeded 100 kg of TNT equivalent (1998). In September 1999 explosions in 
Buinaksk, Moscow, and Volgodonsk took place one behind another. The explosion of the 
apartment house in the Gurianov Street in Moscow of September 9, 1999 resulted in the 


death of 23 men; about 250 were injured, including 73 heavy wounded. 


The explosive device at the center of Moscow was left near the stall with the 
cosmetic goods in pedestrian subway under Pushkin Square (2000). The saleswoman has 
noticed the suspicious plastic bag, left near her workplace, and has addressed to a safety 
service of a mall. The personnel responsible for the safety started to discuss what to 
undertake. Meanwhile, about 5 minutes has passed, and explosion took place. 13 men were 
died, 118 wounded. The power of the explosive device was about 500... 600 g of TNT 


equivalent. 


After analysis of the situation, the experts came to a conclusion: the tragedy could 
be avoided, if the mall safety service had ready for use one of the known anti-explosion 
protection devices. The conclusion is applicable to a set of terrorist explosions occurred not 


only in Russia, but also in other countries. 


The efforts of law machinery should be directed not only on investigation of 


previous crimes, but first of all also on their prevention. 


Major element in accident prevention is the development of mobilization measures 


in case of terrorist explosion threat carried out by law machinery in cooperation with safety 


services of firms and organizations. In Great Britain, for example, “Rules of maintenance 
of health and the safety at work” demand the obligatory planning the mobilization 
measures. They include the use of explosions covers, the application of anti-fragment film 


in window glasses, elaboration of the plans of personnel evacuation, and other measures. 


At ED terroristic threat the counter-measure should estimate the degree of danger 
by the analysis of risk and vulnerability. The experts carrying out the analysis should have 
the experience in such fields as HE chemisrty, structural strength of building, organization 
of safety systems. It is necessary to develop a system of measures for reduction of 
vulnerability of buildings and personnel. However, it is necessary to take into 
consideration, that the measures, which limit the freedom of action too much, will not be 


maintained by the personnel and consequently will not ensure the safety. 


The interest to searching the methods of destructive explosion suppression is not 
weakening for a long time. The need for the decision of this problem exists in army, in law 


machinery activity, and also at realization of industrial blasting operations. 


The intensive development of technological processes involving the energy of 
explosion essentially expands the area of application of explosive substances in industry. 
Thus, considerably increased is the specific volume of technological blasting operations, 
which should be carried out not only at specially equipped testing areas and explosive 
chambers, but also at urban conditions, the workshop rooms etc. The HE charges weight 
restrictions developed for conducting of blasting operations in open-cast mines and testing 
areas, are generally inapplicable, as they do not satisfy to existing ecological safety 
regulations of the energy sources application at urban conditions, on water areas of the seas 


having economic importance, and in other similar situations. 


In the process of development of accident protection devices and safety regulations 
for the technological operations, which use the explosion energy, the allowable effect of 
the shock wave on the object should be determined by criterion of acceptable risk but for 
the criterion of absolute safety. This can be done by comparison, for example, with the risk 
from the effect of the similar natural factors or other kinds energy sources, used on 


practice. 


The most hazardous factor of condensed HE detonation is the shock (blast) wave. 
The analysis of experimental data of the shock wave loads onto the objects testifies about 
the probabilistic character of their damage at dynamic loadings. The probability of damage 
depends both on the individual characteristics of the object, and also on uncontrollable 


factors of dynamic impacts in particular conditions. 


In the book, the modern concepts for the determination of the parameters of the 
blast waves from explosion of HE charges in air and water are presented. The given 
information allows constructing the diagrams of a damage of bio-objects and typical 
structural elements in relation to the parameters of the blast loadings. The examples of the 
damage diagrams in coordinates “pressure — impulse” are shown. Some technical solutions 
directed on the suppression of the blast loadings from HE explosions are considered both 
for the destruction of criminal explosive objects, and for industrial use of explosive 
technologies. In particular, we present the analysis of efficiency of attenuating of the blast 
loadings (pressure and impulse characteristics) with the help of the perforated screens and 


by placing of the explosive sources in elastic envelopes containing gaseous substances. 


The book is constructed on the basis of the original and detailed analysis of results 
of the research works executed at practical application of proposed methods of blast 
loadings suppression by the scientists and the experts of Paton's Institute of Electrical 
Welding, Semenov's Institute of Chemical Physics, Lavrent’ev’s Institute of 
Hydrodynamics, NPO Special Materials, and also on the basis of a number of the foreign 


publications and original publications of the authors. 


The choice of simple and effective ways of blast protection created by explosive 
substances remains an urgent task for fundamental and applied researches. The authors 
analyze the ways of reduction of the blast loadings by placing the sources (HE charges) in 
various explosion protection envelopes or immersing in various explosion absorption 
media. The chosen ways of suppression of the blast waves have found the direct technical 
solution in a series of devices for operative and reliable protection of the equipment and 


personnel at preservation and/or liquidations of explosive objects and devices. 
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CHAPTER 1. EXPLOSION AND EXPLOSIVE PHENOMENA 
1.1. Basic concepts 


The energy of the substance (or the energy of any process) determines the ability to 
perform a mechanical work or to release the known amount of heat. The energy of natural 


processes changes over the wide range represented in Table 1.1. 


Table 1.1 Energy of some processes 


Start of the space shuttle 


Nuclear bomb explosion 


Earthquake 


Star explosion 10°" 


The concept of explosion is closely related to the concept of energy. As a rule, 
explosion is brought into correlation with noise/sound and thermal effects, as well as 
damages. Generally explosion can be considered as the sudden physical or chemical 
change of the state accompanied by thermal and kinetic energy release. The explosion in 


the atmosphere represents the energy release in rather limited volume during a short time 


interval resulting in occurrence of an explosive (blast) wave with finite amplitude [1.1... 


1.14]. The pressure wave attenuates with distance from the source of explosion. 


Different types of explosions are known: 


e Explosion of the boiling devices [1.1]; 

e Nuclear explosion [1.2,1.3]; 

e Explosions in the internal combustion engines; 
e Electrical discharges [1.3]; 

e Self-accelerating chemical reactions [1.1,1.3]; 


e Dust explosions in mines and elevators [1.17]. 


Some explosions have a natural origin: the volcanoes, the lightning. Others are 
intentional ones and therefore controllable: explosions of pyrotechnic substances, industrial 
HE detonation. It is also possible to consider the explosions at industrial man-caused 


accidents [1.1,1.3] as an individual class of explosions. 


Depending on the source of energy it is possible to distinguish physical and 
chemical explosions. At physical explosions there is no heat release, and the source of 
explosion is the energy of compressed gas, dust or two-phase gas-liquid system. The 


example of physical explosion is the bursting of a vessel containing compressed gas [1.6]. 


The chemical explosions are accompanied by the thermal energy release due to the 
chemical reactions. Formally, one can consider here such phenomena as nuclear 
explosions, HE detonations, gas and dust explosions. The most dangerous consequences of 
explosion is the generation of the dynamic impulse in the form of the blast wave resulting 
in the destruction of equipment, buildings and the damage of alive organisms. The thermal 
radiation of explosion products can promote fires [1.1,1.17]. The thermal impulse may 


cause the damages of surfaces, thermal deformations, break of carrying elements, roof fall. 


Table 1.2 Pressure in the products of explosion 


Substance Explosion pressure, kPa 


Milk-sugar 


Sulfur 


Hydrogen 


Methane 


Aluminum 


PETN 


It is useful to remind the magnitudes of pressure in the explosion products for 


various substances, including condensed HE (for example, PETN [1.17]) - see Table 1.2. 
1.2. Energy distribution at explosion 


The parameters of explosion depend on the energy distribution in the field of 
explosion and its redistribution as the blast wave propagates away from the source. Initially 
the energy is stored in the source of explosion in the form of potential energy. At the 
instant of explosion this energy transforms into the kinetic energy and the thermal energy 
of the system, which includes now all substance inside propagating shock wave. System is 
non-steady owing to the increase of the mass of the substance involved in the movement. 
Another reason of unsteadiness is the proceeding redistribution of energy both in the 


explosion products and in the gas subjected to shock compression [1.1]. 


For the analysis of energy distribution we shall accept the following idealization: 1) 
explosions have the strictly spherical form in originally homogeneous unconfined gas 
environment; 2) sources of explosion consist both of energy containing material, and of 
inert confining material, and during the explosion process these substances do not mix up 
with each other or with an environment atmosphere; 3) shock wave is the only factor 
resulting in energy dissipation. In such idealized system the potential energy of a source 
transforms into thermal and kinetics energies of various parts of the system, as well as into 


the energy of radiation. 


Blast wave energy 


The energy of the blast wave, which here is understood as a driven part of gas, 


contains the thermal energy 


and the kinetic energy 
ee [ awa 
eon 


Where: 

pis the density; 

Cy is the specific heat capacity of a gas at constant volume; 
Gis the temperature; 

& is the reference temperature in fresh mixture; 

u is the mass speed of a gas; 


V is the gas volume. This volume does not include the volume occupied by the 


explosion products or by the fragments of confinement. 


At late stage of the process development, when the kinetics energy of the source 
and confining material approaches zero, and the amplitude of the blast wave is small (so 


the energy of dissipation can be neglected) the total wave energy Ey, = Ep + Ex remains 


constant with time. The constancy of the wave energy at the stage of weak explosion is the 


characteristic feature for explosive processes [1.15]. 


Residual energy in the atmosphere (« waste energy ») 


As the explosions are accompanied by the shock waves, and the atmosphere is 
subjected to non-isentropic compression, then after returning of the central part of a system 
to initial pressure, there will be a residual temperature rise. The residual energy (« the 


waste energy ») [1.2] becomes of a constant value at the stage of weak explosion [1.16]. 
Kinetic and thermal energies of the fragments 


Initially the material of confining envelope will move with acceleration and to be 
heated due to heat transfer, friction etc. Then the speed of the fragments will decrease (in a 


limit - down to zero), but the thermal energy partially will be kept. 
Kinetic energy of the source material 


In any explosion the material of the source or combustion products will be involved 
in a movement. Kinetic energy of the source material finally will decrease down to zero, 


when the movement in the near field of explosion will be stopped. 
Thermal (potential) energy of the source 


The source originally contains all the energy of the explosion in the form of 
potential energy. During the explosion the portion of the energy of the source is transferred 
to the other parts of the system, and another portion remains in the source as thermal 
(potential) energy of the explosion products. This thermal energy dissipates by mixing. The 
mixing process is rather slow, and it is possible to accept with good accuracy, that the 


thermal energy of the source is the constant value at late time. 
Radiation 


The rate of energy radiation quickly decreases, so the losses by radiation reach 


constant value at quite early stage of explosion. 


Distribution of energy in the blast wave 


The redistribution of the energy in the blast wave with time is schematically shown 
in the fig. 1.1 according [1.3]. Note, that at the late stage of process, when the explosion is 
weak (far field wave), the total energy consists of thermal and kinetic energies of the wave, 
residual thermal energy, and thermal energy of fragments and products of explosion. 
Besides, some part of the energy has been lost as a result of radiation, however the losses 


by radiation are essential only for nuclear explosions [1.2]. 


Energy fraction , % 
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Fig. 1.1. Typical scheme of energy distribution in the blast wave with time: 1 and 2 - 
thermal and kinetics energies of a source; 3 and 4 - thermal and kinetics energies of 
fragments; 5 - residual energy, 6 and 7 - thermal and kinetics energies of the wave itself; 8 
- energy of radiation; 9 - energy of a wave in a zone of weak explosion 

Only portion of the released energy transforms into blast energy in the far field of 
explosion. The amount of this part must depend on the character of the explosion 
[1.2,1.4,1.5]. For accidental explosions, when the source volume is rather large, and energy 
release is, as a rule, rather slow, it is possible to expect the dependence of explosion 


efficiency on the character of heat release [1.3]. 


1.3. Propagation of explosion 


The blast waves are capable to propagate in all types of medium: air, water, ground 
and solids. The propagation of the air blast waves depends on the blast wave source type 
and the geometrical restrictions, i.e. from the presence of obstacles and friction losses 


[1.10]. 


For the blast wave propagation in gaseous or condensed medium it is possible to 


distinguish three characteristic zones (fig. 1.2): 


Zone of the initial (not burnout) mixture — U; 
Zone of reaction — R; 


Zone of combustion products — Pr. 


AP 


: x ¥ 
PED. ees 
Fig. 1.2. The schemes of pressure distribution for detonation and for deflagration: U/ — 
zone of initial products; R — zone of reaction; Pr — zone of combustion products; AP— 
overpressure. 
Two basic modes of explosive transformation are known: deflagration and 


detonation [1.10]. 


For deflagration mode the movement of the reaction zone on «fresh» (not burned) 
mixture is controlled by diffusion and heat conduction processes. The density, the pressure 


and the temperature vary continuously from zone U up to zone Pr. 


For detonation mode the movement of the reaction zone is controlled by the shock 
compression of «fresh» mixture behind the leading shock wave. The density, the pressure 


and the temperature in fresh mixture have a discontinuity at the shock front. 


In deflagration process, the pressure waves generated in the combustion zone 
propagate ahead of the front. In detonation process, the «fresh» mixture does not receive 
any information because of the supersonic speed of the shock front. The speed of the 
detonation wave is strongly connected to the speed of the combustion zone. The specific 


features of deflagration and detonation are listed in Table 1.3. 


Table 1.3 Properties of detonation and deflagration 


Property Deflagration Detonation 


Pressure Moderate High 


up to 10°... 10° kPa up to 10°MPa 


Subsonic Supersonic 


< 100 m/s > 1000 m/s 


Reaction zone — shock front Diverge interrelated 


Movement of explosion | From the combustion wave | Behind the shock wave front 


products 


The formal description of explosive process as isobaric one does not exclude 
possible effect on the environment. According to the Table 1.3 the speed of deflagration 
can change in a wide range: from several mm/s up to tens and several hundred meters per 


second. 


For the completeness of the description it is necessary to note, that in long channels 
at the presence of obstacles the speed of deflagration can exceed the speed of the sound in 


«fresh» mixture. Such processes were called as quasi-detonations [1.10]. 


The propagation of the explosion in elongated vessels, pipes, and galleries 
frequently is accompanied by pressure fluctuations with clearly distinguished phases of 
compression and rarefaction. Such a behavior is explained by the interaction of the 


combustion zone with the waves reflected from the end-wall of the channel. 


Non-uniformity of combustion speed along the cross-section of the tube and gas 
flow turbulence increase the combustion velocity. The acceleration of the combustion leads 
to the generation of the pressure waves according to the scheme given in the fig. 1.3. The 
compression waves ahead of the combustion zone overtake each other and form the blast 
wave with sharp pressure discontinuity. Under certain conditions such wave is capable to 


generate detonation [1.10]. 


AP 


x 
Fig. 1.3. The schemes of shock wave formation: AP is the overpressure; x is the distance 


The transition of deflagration to detonation is characterized by the jump of the 
speed, which is caused by “explosion into explosion” effect in space between combustion 
zone and the leading shock wave. The effect of “explosion into explosion” is also 


characterized by the retonation wave propagating against the flow [1.8,1.9]. 


Along with the smooth transition of accelerated flame into detonation wave, the 
initiation of detonation by a strong shock wave is possible, as it is realized at use of 


percussion caps in case of industrial HE. 


The classical Chapmen Jouguet model, developed more than 100 years ago, 
predicts the speed, the temperature, and the pressure of detonation products. The model of 
Zeldovich-von Neumann-Doring (ZND) describes the structure of the detonation wave. 


The detonation was observed in materials at various conditions: 


- The gas detonation was observed and described in 1881 by Berthelot, Vieille, 
Mallard and le Chatelier. They obtained the speeds of detonation up to 3500 m/s [1.10]; 


- The detonation in dust suspensions is rare phenomenon in comparison with 
gaseous detonation. The measure of dust to be exploded frequently described by an 
explosion index Ky =(dP/dt)max.V'’. Here: (dP/dt)nax is the highest rate of pressure rise during 
explosion in the vessel of volume V. Indexes Ky are used for the subdivision of dust 
systems according to explosion classes. A coal dust gets into the lowest, first explosion 
class, and the detonation in such system was detected only in tubes with the diameter more 


than 0,6 m [1.17]; 
- Black gunpowder or rocket propellants usually burn down in deflagration mode. 


The explosive transformation of condensed HE is realized in detonation mode [1.1]. 


In Table 1.4 the properties of some HE are indicated. 


Table 1.4 Properties of some condensed HE 
Explosive substances TNT equivalent Density, Detonation speed 


10° kg/m* 


Protective explosives LTS 


ANFO 0,80 
Water gels 1,20 


Emulsions 


Gelatins 


Nitroglycerine 


As an example we specify, that the detonation of | kg HE with specific heat of 
explosion of 5000 kJ/kg and packed in the tube of length | m results in power of 3010° 


kW. It is equivalent to the power of 30 electrical power stations with 1000 MW each. 


1.4 Exothermal reactions, thermal explosions and auto-accelerating processes 


in the condensed media 


The energy during combustion is released by means of the chemical reactions. 
Burning represents the self-sustaining exothermal process. The most important phenomena 
at combustion are the physical processes of heat and mass transfer. Heat of combustion, 
diffusion of chemical products and gas movement are determined by the transformation of 


chemical energy into thermal energy by exothermal reactions [1.10]. 


The explosions are frequently accompanied by the release of gaseous products. 
However it is possible to specify fast and power-intensive explosive processes without the 
formation of gaseous products [1.1]. For example the gasless explosion of thermit 


structures: 
Fe,03 +2Al > A1,03 +2Fe. 


This reaction releases about 4000 kJ/kg of heat, but is not accompanied by serious 


mechanical effects. 


Some exothermal reactions exhibit significant pre-warming-up, capable to convert 


into an open fire [1.3]. As an example we specify the reaction 
2FeS2 a 702 +16 H20 => 2H2SO4 + 2FeSO4 : 7H20 + 1321 kJ. 


With this reaction there is an accumulation of heat, and the self-ignition is possible 


if no cooling exists in a system. 


Other chemical substances are capable to decomposition in the absence of the 
oxidizer at increased temperature. The process can occur with the release (exo -) or the 


absorption (endo -) of heat. The thermal explosion is rather probable under conditions of 


exothermal decomposition and insufficient heat losses. The pressure starts to increase if the 
disintegration of substance occurs with gas production in closed volume. There is a real 
danger of vessel burst followed by the selfignition of disintegration products, if those are 


combustible substances [1 .3]. 


When solid combustible materials are burning, the heat conduction process 
dominates, and diffusion is negligible in comparison with gases. Therefore, the processes 


of heat transfer are basic at burning of solid combustible structures. 


The thermal explosions or auto-accelerating reactions are possible in vessels, 
reactors or power-intensive substances depots. In all cases the thermal explosions are 
connected to the accumulation of heat and the acceleration of chemical reaction. Dynamics 
of thermal explosion depends on the sizes and geometrical form of reacting substance, the 
temperature. If there is a gas release during thermal explosion, then the pressure growth is 
marked and the external explosive effects in the form of blast waves are possible. The 
reasons of auto acceleration of chemical reactions are various, and most principal of them 


are [1.1]: 


e Auto-catalytic processes with the formation of intermediate substances 
accelerating the basic reaction; 

e Sudden exothermal disintegration of half-products; 

e Rapid transition of explosive process into the mass of the not reacted 


«fresh» products. 


The uncontrollable course of chemical reactions at nitration, polymerization is not 
similar to the typical gaseous or dust explosion by a number of reasons. More often the 
auto-accelerating chemical processes occur in the substance with the high initial density at 
the large stored energy. The same processes differ by significant final pressures because of 
large gas generation in the confined space. The two-phase medium (gas bubbles - 
condensed phase) is formed during explosion. Emission of half-products from space, 
where there was a thermal explosion, into surrounding atmosphere quite often becomes the 
reason of the secondary fires and explosions. Here pertinently to note two dangerous 


scripts. 


The first script is realized, when the container (vessel, capacity) bursts under the 
influence of pressure from the heated combustible vapor or half-products of disintegration 
and the contents of the container, quickly evaporating, is thrown out into the atmosphere. 
Such phenomenon is accompanied by moderate blast effect [1.6]. More essential is the 
occurrence of the fireball making a powerful thermal impulse on the environment at 


explosion location. 


Besides, under the influence of thrust from two-phase jets, the parts of the reaction 
vessel are rejected to significant distances and make additional destructions. In the 
literature this type of explosion was called BLEVE (boiling liquid evaporation with the 
subsequent vapor explosion). Finally question on a demolition effect due to BLEVE is not 


solved yet [1.3]. 


The second script of explosive process develops, when the large mass of explosive 
material is thrown out from the confinement and is not ignited immediately. The fuel is 
mixed up with the air and forms a cloud of the mixture within explosion concentration 


limits. Three ways of development of a dangerous situation are possible: 


e The cloud does not find a source of ignition and safely dissipates in the 
atmosphere; 

e The cloud finds a source of ignition, but burns down in a safe, slow mode of 
burning without essential blast effects [1.3]; 

e After ignition the cloud burns down in fast (up to a detonation) mode of 


combustion with essential demolition effects. 


Last variant of explosion is called UVCE (unconfined vapor cloud explosion) [1.1]. 


1.5. Gas and dust explosions 


Typical examples of chemical explosions are gas and dust explosions. The 
necessary condition of their realization is the simultaneous presence of the ignition source 
and explosive atmosphere. Such atmosphere is formed when the air is mixing with the 


combustible gas or the combustible dust in the certain proportions. The upper and lower 


concentration limits outside which the explosive transformation is impossible [1.17] are 


known. 


The state of a gaseous mixture is supported quite long period of time due to the 
diffusion process. The dust accumulates due to the gravitational forces. The combustible 
gas mixture is formed easier, than the dust one. Necessary conditions for dust explosion are 


following: 


e The dust must be combustible, i.e. could burn in air or in oxygen; 

e The distribution of particles of a dust by size should provide the propagation 
of the flame; 

e The dust does not agglomerate at the emission in an atmosphere; 

e The concentration of a dust should be within explosion limits; 

e The concentration of gaseous oxidizer should be sufficient to support 
combustion; 


e The presence of the powerful ignition source. 


It is quite clear, that in real conditions the simultaneous performance of the listed 
conditions in all points of dusty space is not plausible. One of the reasons of formation and 
development of dust explosion is the occurrence of pressure waves at local energy release. 
These waves pass the dust layers causing the additional formation of the combustible 
mixture, either because of friction at the boundary of the dust layer - air or due to the 
vibration under the action of pressure waves. According to the described sequence of 
events the initial explosive process is called the primary explosion, and the subsequent 


stages — the secondary explosion [1.17]. 


1.6. Characteristics of gas or dust explosion 


The typical pattern of pressure history at gas or dust explosion is shown in the fig. 


1.4 [1.17]. 


AP 
A Pow 


Fig. 1.4. The basic parameters of a pressure wave at explosion: AP max — maximal 
explosion pressure increase; tga=(dP/dt) max — maximal rate of pressure increase; t; — 
time of pressure increase; t? — delay time 


Explosion can be characterized by the parameters: 


e Maximal pressure of explosion Pina 


e Maximal rate of pressure rise (dP/dt) max 


The listed parameters are determined experimentally by standard equipment at 
central ignition of the dust in a wide range of concentration. Sometimes in addition the 


time of pressure rise at burning 7, and the time of explosion delay 4 are measured. 


The value Pyyax at explosion in a spherical or cubic vessel does not depend on the 
size of a vessel. The rate of pressure rise (dP/dt)max decreases with the increase of the 


vessel size according to the formula (dP/dt)max V’** = Ke. 


Parameter K,, is the characteristic value of the dust for each substance. According to 
the value K.,, combustible dusts are referred to the various explosion classes, as it is shown 


in Table 1.5. 


Table 1.5 Explosion classes of the combustible dusts [1.17] 


Designations of classes Ky, MPamc’! 


Stl Ln,.20 


St2 20... 30 
St3 > 30 


It is necessary to remember that K,, depends on the conditions of explosion and use 
only the values measured by standard technique. Other factors (some of them are listed in 


tab. 1.6 [1.17]) may influence the value of explosion index. 


For example, the dependence Pax and (dP/dt)ngx on initial pressure for coal dust 


explosion in air is shown in fig. 1.5. 


Table 1.6 Factors, influencing the rate of pressure rise at explosion [1.17] 


Gases Additives of combustible gases and 


vapors 


DM vex > (dP/AD max. 


120 700 
100 600 
80 500 
60 400 
40 300 
20 200 

ie 0,5 1,0 14 _ 


Fig. 1.5. Effect of initial pressure Po on the explosion parameters: — Dependence of Pnax 
on Po; —— Dependence of (dP/dt) max on Po 

The practice shows, that any measures aimed to increase the dust particles size, dust 
accumulation, decrease of the turbulence, elimination of potential sources of ignition really 


reduce the danger of explosion [1.17]. 


1.7. «Ideal» explosions 


Some characteristics of explosions in air can essentially depend on the parameters 
of the source, which forms the blast wave. By definition, the process is considered as an 
explosion, if the energy is released rather quickly and inside the quite small volume [1.1]. 
The energy of explosion /, density of energy release //V and the rate of energy release, 1.e. 
power are the basic parameters of the source, which determine amplitude, duration and 
other characteristics of the blast wave. Is established, that the main properties of blast 
waves can be determined on the basis of unilateral parameter - total energy - without 
taking into account the density of energy release and the source power. The point sources, 
nuclear bombs, laser sparks and condensed explosive substances (HE) are the basic four 
sources of «ideal» explosions. The nuclear bombs and laser sparks are mentioned here for 


completeness of the review and will not be discussed further. 


The explosions generated by condensed HE are useful for investigating the 
accidental explosions, because the majority of experimental data on ideal explosive waves 
are obtained using detonation of condensed HE. These HE generate practically ideal 
explosive waves, because the density of HE is large in comparison with the density of air 
and, subsequently, the energy release per unit volume is significant. For the military 
purposes, the most important characteristic of HE efficiency is the energy content per unit 
mass (or specific heat of explosion), therefore usually in tables both this characteristic, and 


the density of HE are listed. 


The point explosion source is the model, which is convenient for the mathematical 
description of the real explosion process and, hence, has theoretical importance. Point 
explosion is considered as the explosive process, at which the blast wave is formed as a 
result of instant release of finite amount of energy at some point of the homogeneous 
atmosphere. Numerous studies of point explosions were performed for the cases of both 
the atmosphere «of real air», and the atmosphere «of ideal gas», where the ratio of specific 


heats at constant pressure and constant volume is v= 1,4. 


The blast waves are attenuating with distance from the source of explosion, 
therefore it is accepted to specify three basic zones. In the zone, nearest to a source, so- 
called zone of strong explosion, the pressure in the wave is so high, that the external 
pressure can be neglected. The flow is self-similar and is described by known analytical 
dependences [1.7, 1.8,1.9]. The intermediate zone represents significant practical interest 
because overpressure and impulse are rather high to produce significant destructions. 
However the description of the process in this zone has no analytical solution and should 
be done with the help of numerical methods [1.4, 1.6, 1.11]. The far field of explosion (or 
zone of weak explosion) is following the intermediate zone. The analytical solution in far 
field is admitted from the viewpoint that if at some distance the time history of 
overpressure is known, then it is easy to construct the appropriate dependence for the 


larger distances [1.15]. 


There are theoretical indications that in the zone of weak explosion the self- 
sustained N-wave can be observed, in which the process in a positive phase of 


overpressure does not depend on the flow in internal area [1.12]. 


1.8. Basic features of nonideal explosions 


The majority of sources of nonideal explosions has essentially smaller (than the 
ideal sources) density of energy release, final time energy release, and also requires some 
environment for realization of explosion. The explosions of chemical reactors, as a rule, 
are preceded by uncontrolled development of reaction in reactor, until the threshold 
pressure will be exceeded. The explosions of gases and dust in the open-ended volumes 
usually result into fires only. The damages of vessels under pressure with non reacting 
gases can serve as the sources of strong explosions, but such explosions will take place 


only at strength (before damage) construction. 


The researches specify the presence of two factors resulting in nonideal blast waves 
from the spherical source. These are the final volume of the source and the final time of 


energy release in it. 


By calculating the external work at slow adiabatic expansion of sphere to the initial 
pressure level, we shall obtain the formula for effective energy of exploding sphere at 


relative density of energy g=(7./To)-1 [1.13,1.14]: 


: =a [i+a)-+9) | 


Where Q is the total heat released by the source. The value EF; describes the 
maximal work, which can be made by any spherical explosion. In a limit at go for fixed 


Q (point source) £,/O—1, whereas at g—0 and for fixed QO: 


E/Q>(¥%1)/ ¥ 


The essential factor for nonideal explosion is the time of energy release in a source 
t-. Comparing this value to the time of sound propagation from the center to border of the 


source, we shall obtain dimensionless time of energy release 


T1.Co/l ro 


Here cp is the sound speed in the source before energy release, 7p is the radius of the 


source. 


It is necessary to indicate the other characteristics of nonideal explosions. At rather 
slow energy release (for example, at burning speeds smaller 45 m/s, or at dimensionless 
time of heat release r > 0,5) the formation of shock waves is not observed, though the 
positive impulse of pressure in quasi-acoustic compression wave is kept. Other specific 
feature of blast waves from the sources with small density of energy release is that the 
amplitude of pressure of a negative phase becomes comparable to the amplitude of the 
positive phase [1.15]. For this reason the blast waves from the specified sources essentially 
differ from the blast waves of ideal sources. In the blast wave from ideal source there is a 
negative phase, but its amplitude, as a rule, is small in comparison with amplitude of a 


positive phase, and the damage caused by the negative pulse, is insignificant. 
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CHAPTER 2. THE BASIC PARAMETERS OF BLAST WAVES FROM 
CONDENSED HE DETONATION IN AIR 


2.1. Initial status 


The application of condensed HE in practical tasks has resulted in the necessity of 
numerous measurements of the parameters of the blast waves generated by detonation 
products expanding into environment. Such an environment can be a solid body, liquid or 
multiphase medium. Below we consider basically the cases of HE charges detonation in a 


gas - usually in an air. 


What causes the necessity of new revision of a question about the parameters of the 
blast waves generated by the detonation of condensed HE in air? The existence of classical 
measurements of blast waves parameters in air [2.1], unfortunately, does not allow 
considering the case been resolved completely and finally. Moreover, the use of these data 
frequently results in mistakes connected to the restrictions, at which the measurements 


were carried out. 


It is interesting to note, that in the above mentioned reference [2.1] majorities of 
possible restrictions are considered and the appropriate precautions are given. However, 
these precautions were not analyzed subsequently, and, as a result, the most of the given 


expressions were used without any reservations. 


The dependence of blast waves parameters on the power and geometrical 
characteristics are experimentally determined. Form the other side, these dependencies, 
been applied as reference analytical expressions, often did not take into account the errors 
of measurements and areas of applicability, that was emphasized in the appropriate 


primary sources. 


There were no reliable measurements of parameters of air shock waves after 
publishing the works [2.1,2.2,2.3] in the domestic literature. At the same time, the detailed 


materials on the measurements of fields of explosive parameters under various conditions 


of HE charges detonation in air in view of the presence of limiting surfaces were published 


continuously in the world literature. 


Additional collection of actual measurements of air shock waves parameters is now 
gathered in [2.4... 2.25,2.39... 2.41]. On the basis of these measurements it is possible to 
add the results given in the mentioned above sources, and correct them in particular cases. 
Moreover, it is possible to consider additional formulas for calculating the parameters of 
the normal and oblique reflected waves at limiting surfaces. The generalization of well- 
known and new data allows more correctly to estimate the scale of errors arising from the 
use the so-called TNT equivalents for explosions of fuel-air mixtures, dust explosions etc., 


but for the condensed HE only. 
In the literature the various types of condensed HE detonations are analyzed: 


1. Free-air (or simply air) explosion — detonation of HE charge in a gas (more 
often in air) in absence of reflecting surfaces (fig. 2.1, a). 

2. Surface (ground) explosion — detonation of HE charge in a gas (air) in 
case, when the HE charge is placed on a limiting surface (ground) (fig. 2.1, 
b). 

3. Above-ground explosion — detonation of HE charge in air at certain 
distance from a limiting surface (ground) (fig. 2.1, c). 

4. Underground explosion — detonation of HE charge inside the ground. 


5. Underwater explosion — detonation of HE charge, inside the liquid. 


triple 
point 


Fig. 2.1. The scheme of HE charge detonation: a — Air explosion; b — Ground explosion; 
c — Above - ground explosion 


The wave overpressure is usually measured with satisfactory accuracy +10%. This 
fact must be taken into account when performing the analysis of the experimental 
dependences. Note, that other parameters of the waves are measured with greater error that 


must be also taken into account when they are used as the standards. 
2.2. Parameters of shock waves for free-air explosion 


The air explosion of HE is characterized by the minimal distortions of the process 
of explosion products expansion into the environment and, hence, is least complicated by 
various waves refractions and reflections at any obstacles nearby the charge (such, as the 


solid walls or the surface of a ground or liquid). 


In order to minimize the distortion of expansion process, the air explosions, as a 
rule, are carried out with the small or moderate (by weight) HE charges. In [2.1] most of 
the measurements were performed for ground explosion. The sufficient number of correct 
measurements of air explosion parameters is presented in [2.4... 2.9] at shock waves (SW) 
front overpressures AP = P; — Po from 5.5 up to 0.003 MPa. Here P; and Po are the 
pressures behind and ahead of the front accordingly. In the field of small overpressures AP 
< 0.003 MPa the data were obtained on the basis of the analysis [2.7] and some 
measurements in [2.8,2.9]. A lot of information on the parameters of air explosion may be 


gathered from [2.4,2.5]. 


Fig. 2.2. The scheme of measured parameters of the blast wave: t, is the time of arrival; t. 
is the duration of a positive phase, t. is the duration of rarefaction wave; AP nay is the 
amplitude of pressure wave; AP_ is the amplitude of rarefaction wave; I is the impulse of a 
compression phase; I. is the impulse of a rarefaction phase; AP» and Tt - secondary 
pressure rise and its duration 

The diagram of the pressure structure in the blast wave is given in the fig. 2.2 in 
order to obtain clear understanding of the parameters discussed further. The starting time is 
counted from the moment of the blast wave emerges from the center. The value 7, means 
the time of arrival of the wave to the position of measurement at the distance R from the 
center of explosion of the spherical charge of radius Ro. The impulse of the phase of 


compression is calculated using formula: 


t,t, 


I= [rar 


The impulse of the compression phase is the most important characteristic of an air 


wave together with value of overpressure AP . 


To characterize some consequences of explosion the magnitudes of amplitude of 


rarefaction waveAP_ and impulse of a negative phase of pressure 


T= [arar 


tatty 
can be indicated. 


The durations of compression wave Tt, and rarefaction wavet_ are important for the 
analysis of the resonant response of a target under the action of air shock wave. The 
parameters of the secondary pressure increase AP, and T2 are usually ignored, though for 


estimations it is possible to assume with high accuracy that APo = AP_. 
2.3. Shock waves overpressures at air explosion 


Various formulas for determination of overpressure at the shock front are presented 
in published works. Let's outline some of them. Brode [2.10] on the basis of data 
processing of the numerical calculations proposed two expressions for calculation of the 


overpressures: 


AP = 0,657R.* + 0,098 for AP > 0,1 MPa 
AP = 0,096R,! +0,143R, +0,574R,* — 0,002 for AP =< 0,1 MPa 


Here AP , MPa; R, = R/G°*, m/kg”*?; R, m, - distance; G, kg, - weight of HE 


charge. 


In [2.11] three expressions for calculation are given for different intensities of the 


wave (TNT explosion): 
= -1 2 3 
AP =1,38R,' +0,543R.° -0,035R-* For range 0,05 <R* < 0,3, m/ke”™: 


_ -1 ~2 3 
AP = 0,607R,' +0,032R,* + 0,209R, For range 0,3 <R* < 1,0, m/kg"™?, 


AP =0,065R;! +0,397R2 +0,322R° 


For range 1,0 < R* < 10, m/kg®*?. 


In [2.1,2.3] for TNT-RDX mixture (50 % TNT and 50 % RDX) the following 


expression is written: 


for air explosion, in the range 1> R* > 0,1, m/kg””*? 


AP =0,085R,' +0,3R,° +0,82R,° 


> 


or for ground explosion 


AP =0,106R,' +.0,45R,° +1L6RS 


(2.5) 


The above expression is more convenient, than three formulas from [2.11]. The 
coefficients in (2.1... 2.5) depend on the type of HE and must be corrected to take into 
account the heat of explosive transformation for particular HE. For example, the formula 


from [2.3] for TNT air explosion takes a form 


AP =0,081R,' +0,28R,° +0,707R. 


> 


(2.6) 
or for ground explosion 


AP =0,101R,' +0,42R,° + 141k 


(2.7) 


Generalization of the experimental data made in [2.5] in [2.12] for TNT air 


explosion yields 
AP =1,7-10" exp(-7,5R°* )+ 0,0156, 
and for ground TNT explosion 
AP =1,7-10° exp(-7,14R°”* )+0,0156 
At 0,1 <R* <8 m/kg*?. 


Ina far field of explosion, namely at R* > 8 m/kg”*’ for TNT air explosion 


AP =8-10*° exp(- LOTR ) 


> 


For ground TNT explosion 


AP =8-10" exp(-10,46R""), 
All presented dependences are written with the help of Sadovsky — Hopkinson 


variable R. = R/G°*’. Other Sadovsky — Hopkinson variables are r*=1/E°*, 


I, =I/E°* , where E is the energy of explosion. 


Sadovsky — Hopkinson variables are not dimensionless. Therefore all expressions 
for the blast waves parameters in Sadovsky — Hopkinson variables are valid only at that 
pressure of an environment and speed of a sound, at which the measurements are made. 
The dependences of the blast waves parameters on geometrical and power values are not 
universal, since the effects of initial pressure and sound speeds in the gas are not taken into 


account. Sachs variables are more universal for the description of blast waves parameters: 


P=AP/P, R= RPO? / EO? 


> 


33 


- = te = 4 \-0 
t=1c,PO3/E ‘T= Ie, (ER?) 


Here Po, co - pressure and sound speed in an environment. Sachs variables are 
dimensionless, and one can construct the dependences of the appropriate functions P = f 
(Rk); t= f(R); 1 =f (R) suitable for any initial conditions. Quite often the researches use the 
ratio of the distance from the center of explosion RX to a radius of the charge Ro. It is 


possible to show, that there are rather simple relation between three arguments: 


R=R(P/E,) 


-1 
— 4 
R= R/R, c (FapEs 


Here F;, is the specific heat of explosion, is the density of HE. It is possible to 


use the following formulas in calculations for different geometrical sizes: 
R/R,=K,R R,=K,R R/R, =(K,/K,)R. 


The values of the coefficients for some HE are given in table 2.1. Ktnr is the TNT 


energetic equivalent of HE other, than TNT. 


Table 2.1 Geometrical factors for the calculation of the blast fields and the TNT 
equivalents 


HE type p, kg/m? |Eya, J/kg Py =0.1 MPa Py) = 1MPa 


TNT/RDX 


6 
Sale 1660 5,107 - 10 67,37 3,66 31,5 1,72 |1,13 
TNT 1590 4,517 - 10° 64,29 351 30,2 1,65 {1,00 
RDX 1650 5,356 - 10° 68,76 3,71 32,3 1,74 1,18 
TNT/RDX ;: 
AOie 1690 518210 68,50 3,68 31,12 1,73 |1,15 
Nitromethane | 1000 4,517 - 10° 55,2 351 25,6 1,65 {1,00 
PETN 1770 5,800 - 10° 74,2 3,81 34.4 1,74 {1,28 
Nitroglycerine |1590 6,7 - 10° 77,1 41 35,7 188 {1,48 
HMX 1900 5,68 - 10° 753 3,79 34,9 1,78 {1,26 
Nitrogyanidine | 1620 3.00210° 58.8 311 273 1,44 |0,67 


Tetryl 1730 4,517 - 10° 67,7 351 31,4 1,65 {1,00 


HE type pkeg/m — |Eya, J/kg fu= 01M =0.1 — 2u= IMPs = MPa —__ 


Ammonium 


What is the difficulty of application of some of the given dependences for AP 


estimations? 


It was specified in [2.1], that when calculating the argument R* one must use not 
complete weight of the charge G, but only the weight of its so-called «active» part Gz. The 
value G, decreases with the decrease of G. There are also specified the values of @ = G/Ga 
for TNT. At G= 1 kg we have a= 30 %, accordingly at G=8 kg, a = 65 % and at G= 
kg, a = 80 %. In the case of amatol there is additionally dependence on the distance from 
the center of explosion (or from the wave overpressure). The corresponding data are listed 


in tab. 2.2. 


Table 2.2 Relative weight of an active part of a charge (amatol), % 


Weight of a Pressure in the wave, MPa 
charge, kg 2,0 
25 60 
73 
83 


As can be seen from above presented information, only at G> 50 kg it is possible to 
apply the formula from [2.1] to estimate AP without any exceptions. Even larger limit G > 


200 kg is indicated in [2.1]. 


Let’s analyze the type of curves AP/Po = f(R+), shown in the fig. 2.3 for various 


expressions, taking into account the above remarks. The line | is plotted using the 


Sadovsky formula, line 2 — Baker — Held dependence from [2,4], line 3 using expression 
(2.2) from [2.11], line 4 - Brode formula [2.10], and line 5 - are the data from [2.6]. As it is 
seen at AP/P) > 40 the formulas of Sadovsky and Brode are not suitable for estimations 
and deviate from the data of experiment. It was already pointed out in [2.3]. One must take 
into attention rather wide disorder of the AP measurements, which was also underlined in 
the review [2.4]. As an illustration, some more dependences are given in the fig. 2.4. Line 


1 is obtained in [2.7], line 2 - in [2.8], line 3 - in [2.9]. 


The analysis of the presented data shows, that the accuracy of the formulas from 
[2.1, 2.11] in the field of their applicability in definition of pressure does not exceed 15 %. 


A smaller accuracy of AP measurements (up to 30 %) was pointed out in [2.10]. 


The specified accuracy is satisfactory at the direct estimation of overpressure by the 
value of explosion energy in view of the remarks from [2.1]. However one can obtain 
significant error (up to 100%) in the backward expert estimations of the released energy by 


overpressure (as in cases of TNT equivalent use). 


The error is especially great in the practically important far field, at R* > 1 m/kg'?, 


where E ~ (APR). 


Finishing the consideration of expressions for overpressure in a wave from air 
explosion, we shall present the data for the amplitudes of rarefaction waves [2.6]. In table 
2.3 the amplitude of rarefaction wave in comparison with the amplitude of the shock wave 


are shown for various R*. 
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Fig. 2.3. Dependence of the shock wave intensity on the distance. I - by [2.1]; 2-by [2.4, 


-s 
2.12]; 3 - by [2.11]; 4-by [2.10]; 5-by [2.6] 
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Fig. 2.4. Dependence of the shock wave intensity on the distance: 1 - by [2.7], 2-by [2.8], 


3-by [2.9] 


Table 2.3 Parameters of rarefaction wave for air explosion 


R* mike? AP/Po AP/Po 


It is possible to derive the approximated formulas for the calculation of the 


rarefaction wave amplitude AP_/Py = 0,48 R.' at R, > 1,2 m/kg' 3 For R, <1 m/kg! 8 AP. 


/Po ~ const and AP. ~ 0,1 MPa at Po = 0,1 MPa. 


The important addition to the listed data on intensity of shock waves is the data 
from [2.17] for HE charges with various densities. It was shown in [2.17], that the energy 
similarity in the form [2.1] takes place only at R* > 0,8 m/kg"*, when the front pressure 


does not depend on HE properties. In the near zone 0,053 < R* < 0,8 m/kg"” 


the energy 
similarity is not valid, and the geometrical similarity of explosions is kept only at identical 
parameters of a detonation. For the description of parameters of shock waves in a near 
zone it is required to adopt additional parameters describing the rate of heat release, such 


as the pressure or detonation speed as offered in [2.17]. 


Note, that the measurements in the specified reference can be used for 


determination of the overpressures in the near zone of explosion at 4 < AP/P, < 200. 
2.4. Shock waves intensity for ground and above-ground explosions. 


The rule of doubling of explosion energy in any of the above presented formulas is 
often used to transfer the consideration from the air explosion case to the ground explosion 


case. The experiments show, that this rule is not universal and has a bad accuracy at the 


distances R < 25Ro from the center of explosion. Let's compare the energies for ground and 
air explosion by the impulse of compression phase and by overpressure using the data from 


(2.5, 2.6]. 


In the fig 2.5 the dependence of TNT equivalent of ground explosion in relation to 
air explosion (by impulse) as a function of dimensionless overpressure in a wave is given. 
It is seen, that even at AP/Py > 2 Ki > 2. At AP/Po = 10 according to the data [2.5] - (line 1) 
Kk; © 3,8, and according to the data [2.6] - (line 2) Ky ~ 3. The reason of deviation from 


equality Ay = 2 is, that at AP/Po ~ 10 and R/R, < 20 the behavior of the wave is 


influenced by the expansion of reaction products which, apparently, is unequal for ground 


and air explosions in the near zone. 


In the fig 2.6 the dependence of TNT equivalent by pressure for ground explosion 
in relation to air explosion [2.6] is presented. It is seen, that Kp is within the limits of 1,6 
... 1,8. Even more complex is the dependence of TNT equivalent by pressure and impulse 
for the above-ground explosion in comparison to air explosion. In the fig 2.7 the 
dependence of TNT equivalents for the above-ground explosion from overpressure in a 
wave is given. The lines 1, 2 correspond to TNT equivalent by pressure for A, = H/G°** = 
2 and 0,44 m/kg’** (#- is the height of explosion), line 3 and 4 — TNT equivalent by 


0,33 


impulse for heights 0,4 and 3 m/kg’~” . The line 5 indicates the change of TNT equivalent 


by impulse between the above-ground and ground explosions at \, =0,6 m/kr°. For the 
above-ground explosion at 0,1 < AP < 1 MPa and Po = 0,1 MPa its efficiency is higher, 


than in the case of ground explosion. 


AP/Py 


Fig. 2.5. Dependence of TNT equivalent by impulse for ground explosion in relation to air 
explosion on the intensity of the wave: 1 - by [2.5], 2-by [2.6] 


AP, ies 


Fig. 2.6. Dependence of TNT equivalent by pressure for ground explosion in relation to air 
explosion on the overpressure at the wave front [2.6] 
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Fig. 2.7. Dependence of INT equivalents by pressure (1, 2) and by impulse (3, 4, 5) for the 
above-ground explosion on the overpressure at the wave front (A, = 2 (1); 0,44 (2); 0,4 
(3); 3 (4); 0,6 m/kg”*? (5) 

In [2.16] the generalized dependence of TNT equivalent of any explosion on the 
height of the charge position above the surface of the earth is reported. In the fig. 2.8 these 
dependences are reproduced. Specific feature of the given dependence is the fact, that 
already at insignificant height of the charge H/G°** > 1,25 m/kg°** explosion cannot be 
considered as ideally ground. It is possible also to claim, that the ideal case of air explosion 
is realized at H/G°*? > 8 n/kg°*’, ie. the presence of limiting surfaces results in 


appreciable distortions to the given measurements. 


0 15 30 45 60 7,5  9,0A%,m/k 


Fig. 2.8. Dependence of TNT equivalent by pressure on An = H/G°* for the above-ground 
explosion (by data [2.16]) 

When performing the estimations of pressure fields the researchers frequently do 
not make distinctions between strong nuclear explosion and HE detonation. The existing 
data on pressure fields demonstrate the inconsistency of such approach. In the fig. 2.9 the 
dependence of TNT equivalent of strong nuclear explosion relative to HE detonation is 
shown. Depending on the intensity of the shock wave at AP/P» < 20 the TNT equivalent of 
strong nuclear explosion is reduced to 0,45... 0,60 [2.14, 2.15]. In [2.6] this circumstance is 
connected to the greater loss of energy by radiation at strong explosion. When explosion 
takes place nearby the limiting surface, the pressure in a wave normally reflected from a 


rigid wall worth to be measured. Using the data [2.5,2.6] for a rough estimation of pressure 


0,33 
g 


in the reflected wave it is possible to offer an expression P,/Py ~ (70 Ro/R)* at 8 < R/Ro < 


60. Here R is the distance to the reflecting surface. 


For the above-ground explosion the pressure P,, strictly speaking, is realized 
precisely under the center of the raised charge. At some distance from a vertical connecting 
the center of the charge to a ground, the pressure will vary. Therefore, for the above- 
ground explosion it is necessary to distinguish the zones of regular and irregular reflection. 
In the zone of regular reflection the gauges on a ground will record the pressure behind the 
oblique reflected shock wave. In the zone of irregular reflection the gauges will record the 
parameters of the wave corresponding to Mach stem. This wave is close by intensity to the 
wave from the air explosion with the corrected value of energy. In the zone of regular 
reflection and in the space above a reflecting surface, the gauges at first will record the 
parameters of the passing wave, and after that — the parameters of reflected wave. The data 
[2.5, 2.6] allow specifying the division of zones of regular and irregular reflection at the 
above-ground explosion. The trajectory of the triple point is shown in the fig.2.10 in 


Sadovsky — Hopkinson variables by height and radial distance. 
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Fig. 2.9. Dependence of TNT equivalent of nuclear explosion (by pressure) on the 


overpressure 
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Fig. 2.10. A trajectory of the triple point for the above-ground explosions at a various 


heights of charge position 


In the fig.2.10 H* = h/G°*? is the height of a triple point, and Hy/Ro is the same 
height relatively to the radius of the charge. The trajectories are given for two magnitudes 


of parameter A, = 0,5 and 0,8 m/kg°*? , Le. at the heights of explosion H/Ro = 9 and 15. 
2.5. Velocity of blast waves 


Using the pressure measurements one can determine the blast wave velocity. 
According to the data [2.5] at 1.5 < R/Ro < 70 the velocity of the shock wave can be 
estimated as M ~ (70 R/Ro)***. Here M = D/co is the Mach number, D is the speed of the 
wave, Co is the sound speed in unperturbed medium. There are also slightly different data 
in the literature. Thus, according [2.2] for the detonation of charges TNT/RDX 50/50 M = 
27,9 (Ro /R)°” at 1 <R/Ry < 10, for the explosion of charges PETN M= I7R IR! at 1 = 
R/Ro < 10. In the fig. 2.11 the corresponding dependences are compared. The line | is 
obtained for the speed of the wave [2.5], the line 3 for RDX [2.2], line 2 corresponds to the 
speed of the gas at the front of wave W/co [2.5]. 
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Fig. 2.11. Dependence of the shock wave Mach number on distance: 1-by [2.5], 2-by [2.5]; 
3-by [2.2] 


The magnitude of the dynamic head behind the wave is related to the speed and the 


intensity of a wave. The ratio of dynamic head to initial pressure O = g/Po can be estimated 


using the approximate expression from [2.5]: O = (35Ro/R)*® at 7 < R/Ro < 70. The 
dependence of Q on the R/Ro for TNT is given in the fig. 2.12. 
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Fig. 2.12. Dependence of the dynamic head and the time of wave arrival on distance: 1 — 
QO by [2.5], 2 — 1, by [2.5] 


2.6. Impulse characteristics of blast waves 


Let's analyze the relations for another important parameter of the air shock wave — 


the impulse of the compression phase. 


In case of impulses it is necessary to distinguish zones, where the movement of 
expanding products of explosion influences the pressure fall behind the wave front. 
According to the data [2.3] the width of a layer between the products of explosion and 


front of a wave and the radius of the pressure wave XR, will be: 
(SRS" )( Rs / Ro) ~ 0,045 (Rs / Ro )'* -1) 


Or at Rp/Rs << 1 5R,! = 0,045 (Ro/Rs) °. 


Even at the distance from the center of explosion Ro/Rs << 35, i.e. at R* < 1,8 
m/kg?** and AP/Py) > 2,6 SR;' ~ 0,25. The products of explosion occupy an essential part 


of volume involved into the movement behind the shock wave. Frequently, the calculations 
of the impulse of compression phase use the expression from [2.3] for air explosion of 


TNT/RDX 50/50 charge 
1, =20K; Je Great Rs RG 
And for detonation of TNT charge 
I, =18,4R;,' 


In these expressions / is measured in kgf s/m’, and G in kg. However, authors [2.1] 
present the similar formulas and make the caution about the low accuracy of the factor of 
proportionality. The further experiments have confirmed the insufficient accuracy of the 
given formulas. The results of measurements of impulses from different references are 
given in the fig.2.13 in coordinates /* and R* or AP/Po. The line | is plotted by expression 
I* = 20/R+, line 2 - by measurements of Adushkin and Korotkov [2.2], and line 3 - by the 
data [2.5, 2.6]. The experimental points (solid circles) correspond to the measurements of 
the pressure impulse in the rarefaction phase by the data [2.6]. As it is seen from the 
fig.2.13 the empirical formula /* = 20/R* for TNT/RDX 50/50 has the best accuracy. 
Therefore the formula /* = 20/R* is suitable only for the estimations of the impulse in 


rarefaction phase with accuracy ~ (20 ... 30) %. 


T 


Fig. 2.13. Dependence of the impulse of pressure wave on distance and intensity of a wave: 
1 - by [2.3]; 2 — by [2.2]; 3 — by [2.5, 2.6]; *— by [2.6] for rarefaction phase 

In the far field of explosion at AP/Py < 2.3 the impulses of compression / and 
rarefaction /. waves are close to each other, that for the first time was marked in [2.1]. At 
AP/P, > | the value of pressure impulse in the wave should be estimated by the corrected 
expressions. So at 10> AP/Py > 1 the factor of proportionality in the formula I* = A/R *, A 
< Aj, where A, from [2.1, 2.3]. For TNT/RDX 50/50 J, = (15 + 2)/R, is closer to 


experimental data, and for TNT /, = (13,8 + 2)/R,. Applying the energetic modeling for 


ground explosion relative to air explosion, the impulse of compression in a wave at ground 


explosion will be /, = (20 + 3)/R,... 


In the near zone AP/Pp > 10 the values of the measured impulse quickly reach a 
maximum and does not grow by approaching to the center of explosion. The influence of 
explosion products expansion in this zone is reduced, so a simple rule of power modeling - 
«doubling of energy of explosion» is not valid anymore. On a descending branch of 
dependence / = / (R, E), i.e. at AP/Po < 10 the dimensionless impulses can be calculated in 


Sachs form by the formula 


1 ~1,8R,/R 
Or dimensional impulse by the formula 


1=1,8(P,)” (£°")R,/eoR 


In a number of works the change of the pressure impulse in a wave at normal 
reflection of a wave from a rigid wall was studied. Korotkov and Tsykulin [2.13] have 
shown for the first time, that the increase of the pressure impulse in reflected wave 
corresponds to the change of pressure. In table 2.4 the ratios of impulses in reflected and 


incident waves [2.13], and overpressures in the waves are listed. 


Later, the measurements [2.13] were confirmed in [2.5]. It is possible to offer an 
expression J/,/7 ~ AP,/AP [2.5, 2.13] for the calculations of pressure impulse in the normally 
reflected wave. Here /,, AP; are the impulse and the amplitude of the reflected wave. The 
validity of equality /,/7 ~ AP,/AP means that the durations of phases of compression for 
incident t and reflected wave t, are close and with good accuracy it is possible to take t = 
t,. Because of the increase of the impulse at reflection it is possible to introduce the 
concept of TNT equivalent by impulse relative to the air explosion. In the fig.2.7 the lines 
3, 4 indicate TNT equivalent of the above-ground explosion for reference heights AH = 0,4 
and 3 m/kg’"*. The line 5 shows the change of TNT equivalent by impulse between the 
above-ground and ground explosion at Au = 0,6 m/kg””’. It is clear, that at such reference 


height the explosion is almost equivalent by impulse to ground one. 


Table 2.4 Parameters of the reflected and incident waves at normal reflection at the wall 


3,9 5,1 0,75 
3,4 42 0,36 
3,2 3,4 0,22 


3,0 mk 0,15 


27 27 0,096 


2,65 25 0,065 


The generalization of dependences for intensity of a shock wave and impulse of 
compression phase allows specifying universal link of the basic parameters of the blast 


wave and the energy of explosion. On the basis of the analysis of numerous measurements 


it is possible to find the following relation 


T= o(AP/P,)" 


Here / - dimensionless impulse in Sachs variables, q - factor of proportionality. 


For the data in [2.5] @ = 0,04. 


In dimensional form 


For a constant parameters of an environment, for example co = 340 m/s and Po = 
0,1 MPa 


J22 3107 APR 
Here /-Pas,AP-Pa,E-J. 


_ 7p-0,33 — 0,44 
In Sadovsky — Hopkinson variables [,=1E 0" = BAP” 


Dependence of the dimensionless impulse on the overpressure 1s shown in fig. 2.14 


at 


Fig. 2.14 Dependence of the dimensionless impulse on the overpressure 


The proposed formula T =0,04 (AP/Po) 44 describes a branch of a curve at AP/Po 


< 10. 
Dependence of a type 
ge n 770,33 
[=k AP"E™ where 0,4 <n < 0,6 


is suitable for other types of explosive transformation. So in [2.14] for strong 


2 * 4 p0.55 770.33 
explosion the validity of expression [=k,AP FE” was shown at AP < 400 MPa. For the 


gaseous detonation it is possible to specify the dependence 


T= EAP ES 


Let's analyze time characteristics of the blast wave. To calculate the duration of a 


phase of compression the expression t* = 1,2R«°° was offered in [2.1]. Here t - ms, 


= 0,33 
%=t/G™ Ge kg HE. 


In [2.12] another expression for r —0.06 exp(1.4R?*) was given at 0,1 < R* < 1,3 


ee 0,18 
m/kg"* and: <= 13,8—-14K, at R* > 1,3 m/kg!”. 


In the fig. 2.15 the dependences [2.1] is given by a line 1, and [2.12] — by line 3. 
The line 2 - is plotted by the results of measurements [2.3]. The greatest divergence of t* 


values in a zone R* < 1 m/kg'”® 


where, apparently, the data [2.3] and dependence from 
[2.12] are closer to the experimental data, than formula from [2.1], that was also specified 
by the measurements of t* in [2.5]. At R* > 2 m/kg"”, i.e. at AP < 0,2 MPa all formulas 


lead to the close values of t*. 
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Fig.2.15 Dependence of the duration of the compression phase on the distance: | - by 
[2.1], 2 - by [2.3], 3 - by [2.12] 


Using the approximate equality of pressure impulses in compression and rarefaction 


phases at AP/Po < 10 it is possible to specify the ratio of durations of these phases. The 


dependence of the ratio of durations of phases of rarefaction (t.) and compression (t.) on 
distance in the Sadovsky — Hopkinson form is presented in the fig.2.16. At R* <2 m/kg!”, 
i.e. at AP > 2,3 Po the ratio (t-/t+) is ina range 8 ... 14. In connection with elongation of 
the rarefaction wave in time and its small amplitude, the effect of the rarefaction wave on a 


target is small in comparison with the pressure wave [2.1]. 
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Fig. 2.16 Ratio of durations of rarefaction and compression phases in the blast wave 


Another parameter of explosive process is the time of arrival of a shock wave from 
explosion to the given location. As a first approximation for estimations it is possible to 


use an expression 


7, =10°(R/8R,)™ 


Here 7, =1,¢)(P,/ Ey - time of arrival of a wave for distance R in Sachs 


variables. In dimensional form 


0,33 1,74 
: col) \ 8Ry 


Dimensions f, - 8s; E — J; Po - Pa. 
The formula has satisfactory accuracy at R/Ro > 8 and is plotted in the fig.2.12 


2.7. Wave perturbations in the homogeneous atmosphere at the far distances 


from the center of explosion 


Due to the attenuation of the air shock wave at far distances from the source of 
perturbation, the speed of propagation smoothly reduces to the sound speed in an 
atmosphere. At quasi-acoustic and acoustic stages the pressure wave depends on the 
condition of the atmosphere and, first of all on the gradient of the temperature field and the 
speed of a wind. It is considered, that at the distance, where the wave overpressure reaches 
the magnitude AP ~ 2,5 kPa, i.e. at R/Ro > 500 it is possible to expect appreciable 
distortions of the wave front due to the effect of meteorological conditions. According to 
the correlation for acoustic pulses AP ~ po co u (where 9 = 1,29 kg/m’, co = 340 m/s are the 
density of air and the sound speed accordingly) the overpressure AP < 2,5 kPa corresponds 


to the speed of air movement behind the wave 6 m/s. 


Such speed of air movement is comparable to the speed of atmospheric movements, 
which explains the appreciable influence of an atmosphere on wave process. Depending on 
the condition of an atmosphere the movement of the wave front can deviate either into 
direction out or to the ground surface. For the specified reason the overpressure at the far 
distance from the center of explosion in a non-uniform atmosphere can change essentially. 
The expression for calculation of overpressure at the far distance from the center of 
explosion is given in [2.16] on the basis of the experimental data At ground explosion in a 


homogeneous atmosphere the overpressure versus the distance & will be [2.16]: 


AP = 668(2G)"" R™(P/B,)” 
Here AP - in kPa; G - weight of TNT charge, t; Po - atmospheric pressure at the sea 
level - 102 kPa; P - atmospheric pressure at the position of explosion, kPa, R - distance in 


km. 


The factor of doubling of the charge weight is considered only at ground explosion. 
According to the data measurement [2.16] the real air explosions without influence of a 


Ws 


surface are realized only at heights H/G°** > 6 m/kg 


Sometimes, it is useful to apply the other scale of measurements at the analysis of 
the magnitudes of quasi-acoustic perturbations, created by HE detonation. The detonation 
of HE charge at far distance is considered as a pulse of low-frequency sound signal. The 
overpressure AP in Pa corresponds to the intensity of the sound signal in dB by relation A 


= 20 log (Ap) + 94. For convenience of recalculation the diagram is presented in the fig 


2.07: 


| 


AP WiPa 


Fig.2.17 The diagram of recalculation of intensity of a sound signal by the overpressure. 


It is possible to specify the dependences for the amplitude and duration of a shock 
wave at the acoustic stage given in [2.24]. Thus, at explosion of HE charge on a soft 


ground the amplitude of the air wave will be 
41/3 14 
AP =4,4(G'?/R) 


At 1000> R/G'? > 10 m/kg"*; R - m, G - kg, AP - bar. 


At explosion of charges on a rocky ground 


1,5 


APES GTR) 


in the same range of lengths and weights of the charge. 


The duration of compression phase in ms is: 


r, =3,15G"" Ig(R/G"*) 


2.8. Wave perturbations in a non-uniform atmosphere 


The gradients of wind speed and the temperatures, which are characteristic features 
for an atmosphere, affect the propagation of quasi-acoustic pressure impulses. The 
reducing of air temperature with height, characteristic for daytime, promotes a deviation of 
wave front from the earth surface and attenuates the perturbations in comparison with 
predicted values of AP. In an inverse atmosphere, where temperature of air grows with 
height, the wave front turns to the surface of earth and the level of perturbations can be 


larger than predicted values of AP. 


The increase of the wind speed with height results in a deviation of the wave front 
from the ground during the propagation of the wave against a wind and the turn of the 
wave to ground during the propagation of a wave alongside with a wind. Thus the level of 
perturbations in the direction against a wind is less than predicted AP, and vice versa: they 
are higher in the opposite direction of wind propagation. The marked fluctuations of waves 
parameters are so essential, that frequently are taken into account at blasting operations, 


especially with the large HE weights [2.24]. 


It is well known, that atmospheric pressure and temperature change with from a sea 
level. Let's consider the dependence of the basic parameters of demolition action on 
atmospheric conditions. Let the amplitude of the blast wave from explosion under the 
normal conditions (Py = 10° Pa and 7» = 288) is APo, impulse Jo, time of arrival to the 


reference point fo, duration of phases of compression and rarefaction t,o and t.o. After 


rearranging dimensionless parameters of similarity at temperature 7,,/7) and pressure P,,/Po 
we obtain: 


AP, /AP, =(P,/P,) 


> 


n A 


Tdd (PP P.Y ORs Tey 


TL, [T,, = caer [t4 = Le Ita = ee /P, oe (288/7, ve 
All parameters with an index ” correspond to explosion at P,, and 7;,. As seen, the 
decrease of initial pressure leads to the decreasing of overpressure and impulse, transferred 


by the wave. 
2.9. Cup formation at explosion 


Underground, ground and above-ground the explosions of HE charge result in 
emission of the ground and the formation of a cup (crater). The sizes of a crater depend on 
the large number of the factors. The most essential ones are the energy of explosion and 
charge location relatively the surface of the ground. H = 0 - ground explosion, H < 0 - 
underground explosion and H > 0 - above-ground explosion. The type of the ground is 
essential, above or inside which the explosion takes place. The parameters of relative 


efficiency of ground emission for a number of HE are collected in table 2.5. 


Table 2.5 Relative efficiency of HE detonation on emission 


HE Type Weight factor of equivalence 


TNT 


Amatol 


Dynamite 


Pentolite 


Nitromethane 


RDX 1,34 


> 


The form of a charge and HE type in comparison with depth of the charge location 
and type of the ground are unimportant. The sizes of a crater can vary in ten times in case 
of transition from one type of ground to another. The generalization of the information has 


resulted in the empirical formula [2.16] 


V, IV, = Gexp|-5,2H(V,G)° | 


Where JV, - volume of a crater, m’: 

Vo - characteristic parameter of a ground type, m’/t: 

G - Weight of HE charge (TNT), t; 

H - Height (depth) of a charge location above (under) the ground, m. 


It is possible to calculate characteristic radius XR, and depths of a crater D,, using 


simple relations 


Resi 2 
D,=05V," 


The typical values of the parameter Vp are listed in table 2.6. 


Table 2.6 Characteristic parameters of crater formation efficiency for different soils 
Soil 


Damp clay 


Damp fine sand 


Damp large sand 


Damp shelf silt 113 


Dry clay 


Sand dry 


Sandstone, tufa 


Granite, basalt 20 


The accuracy is + 30%. 


After explosion the most part of the material, which has been thrown out from a 
crater, is concentrated at the distance R < (3 ... 5) R, from the center of explosion. The 
ground explosion of HE charge is accompanied by displacement of ground [2.20, 2.21]. 
The maximal speed of displacement for non-buried charge can be found by the formula 


from [2.16] 


Dyyx = 2700(R/ G28)" 


max 
> 


Where Dyna, in sm/s, X - in meters, G - weight of an equivalent TNT charge in tons. 


10*(R/Go@ 


For buried charge Dina ~ . The peak value of ground acceleration in cm/s” 


is about Gmax ~ 10°(R/G? i : 

It worth to note, that the formation of a crater at the ground and above-ground 
explosion is a characteristic feature of condensed HE detonation [2.42]. The critical value 
AP* ~ 2,0 MPa of the wave front overpressure required for the emission of the ground is 
given in [2.9]. At explosion of fuel - air mixtures (FAE) such level of pressure cannot be 
achieved anywhere, including inside area of the detonating cloud. For the specified reason 
the FAE above the surface of the ground is not accompanied by the cup formation, and the 


messages, known in the popular literature, do not correspond to the truth. 
2.10. Sizes of a cloud of explosion products 


The thermal energy released in HE explosion, results in emission of hot products of 
explosion (EP) into surrounding atmosphere. These products involve into the movement a 


part of the fragmented material of the soil. The EP cloud and dusty substance rises in an 


atmosphere and drifts with a wind, leaving a dusty trace on the earth ground. According to 
[2.16], in 120 s after explosion the top border of a cloud will reach the height H = 506G°”° 
m, here G is the weight of an equivalent TNT charge in tons. Finally, at the moment of 
cloud hanging, the top border will be at height H; = 670 G°’, and the bottom border of EP 
cloud will be at height H, = 335G°”’. Thus the diameter of a cloud is Dn = 335G°”°, and 
the volume Vn = 3° 10’ G?” m’ = 3 107 G°” km’. 


2.11. Afterburning of detonation products 


In spite of long — time investigation of the dangerous factors accompanying the 
free-air explosion of condensed explosive substances, a number of questions left which 
require additional consideration. One of this poorly known parameters of explosion effect 
concerns the amplitude - temporary characteristics of the thermal impulse due to the 
chemical interaction of expanding (after the detonation of initial charge) products of 


explosion with an environment (for example, with air). 


Conditionally it is possible to distinguish internal and external changes of 
temperature at HE detonation. The internal changes of temperature in the charge volume 
are caused by detonation complex passing through it and occur for a short-term interval of 
time with duration Ar < 10° s. To study these temperature changes the special methods of 
optical pyrometry [2.36] are offered. In HE detonation studies these methods are directed 
only to the obtaining the information about the nature of radiation from the detonation 
complex, temperature of detonation products, and the change of the diagram pressure — 


time in the charge itself [2.36]. 


The external changes of temperature are caused by the expansion of EP and their 
mixing with air. The afterburning of products of incomplete transformation is separated 
from the act of the detonation by an interval of time © >> 4? and occurs in an environment 
much later than its shock compression by the blast air wave responsible for a demolition 
defeat. In the literature there is no information on the amplitude - temporary parameters of 


external temperature fields or molecular or mixed HE of various structure. 


In recent time the complex investigation of the damage factors from the explosion 
of systems with time-delayed energy release due expanding EP afterburning becomes of 
growing interest. In this case (against to the standard schemes HE detonation study) the 
technique of registration of governing parameters, such, as EP temperature and intensity of 
radiation in the given ranges of wavelengths at all stages of explosion - from initial (with a 
time scale ten micro seconds) up to late, of ten milliseconds after the start of explosive 


process is necessary. 


The temperature and intensity of EP radiation for various stages expansion of 
detonation products for TNT and RDX based HE were registered in [2.36] for the first time 


with the help of original double-wavelength photoreception devices. 


As expected, the explosion of molecular HE is accompanied in general case by two 


temperature impulses (T-impulses):, 


The 1-st T-impulse has the duration up to 10“ s and is caused actually by the HE 
detonation itself; the danger of a defeat due to radiation from the first T-impulse can not be 


taken into account in view of its small duration; 


The 2-nd T-impulse has duration up to 10°' s and is caused by process of HE 
burning in surrounding air; the danger of radiation from the second T-impulse should be 
taken into account at the complex estimation of consequences of HE explosion because of 


baric-thermal nature of blast loading. 


The relation of amplitude - temporary parameters of the 2-nd T-impulse from the 
chemical structure of explosive substance is established for the first time in [2.36]. The 
explosion of molecular HE of RDX type (oxygen balance OB = - 21,6) differs 
insignificantly by secondary increase of temperature by amplitude. The consequences of 
explosion basically are connected with baric loading. The explosion of molecular HE of a 
TNT type (oxygen balance OB = - 74) is characterized by essential secondary increase of 
temperature. The estimation of consequences of such HE explosion should be complex and 


must take into account the thermobaric effect on the target. 


The explosives under study were TNT, RDX, mixture 50% TNT + 50% RDX (TG 
50/50) and mixture 95%RDX + 5% paraffin (RDX(p.)). Table 2.7 presents some properties 
of these high explosives: p — initial density, D and p — detonation velocity and pressure, 7, 
— temperature of detonation products and C — content of carbon (soot) in detonation 
products (calculated by [2.36]). The composition was calculated for the expansion level in 
2-3 times. All the charges of mass 100 g were cylindrically shaped with diameter 40 mm 
prepared by cold press technique. Initiation of detonation was performed by means of 


blasting cap and intermediate 8 g RDX charge (diameter 20 mm). 


Table 2.7 Properties of high explosives 


Figure 2.18 presents relative displacement of the charge and photodetectors (PD) in 
the explosion chamber. The distance between horizontally hanging charge and the floor 
was 0.94 m. One of the photodetectors was positioned along the axis of the charge, while 
another PD was focused on the sidewall of the charge. The main element of the 
photodetector is Si-Ge photodiode of sandwich structure containing the thin semi- 


transparent silicon layer that covers germanium substrate. 


The essential feature of Si-Ge photodiode is a very small intersection between the 
Si and Ge spectral responses. The advantage of the PD is the possibility of measuring 


radiation simultaneously in two portions of spectrum, namely 0,4 < A < 1,1 pum (Si- 


channel) w 1,0 < A < 1,8 um (Ge-channel). The photodetectors are equipped by specially 
designed preamplifiers. The preamplifier features low-noise high-speed operational 
amplifiers AD823. Since the absolute values of the expected radiation fluxes can change in 
an extremely wide range (up to five orders of magnitude variation), special attention was 
given to the possibility of changing the proper measuring range. The preamplifier is 
designed to operate at six sensitivity values. The faster time response achieves | ps. PC- 
based analog-digital converter T512 was used to register the outputs of the PD units. The 


photodetectors were equipped by objectives thus providing view angle 0.04 rad. 


Fig.2.18 Experimental arrangement. | — charge; 2 — intermediate charge; 3 — 
blasting cap; 4 — photodetectors; 5 — objectives. 


Fig.2.19 Outputs of Si channel (arbitrary units) of the PD placed perpendicular to 
the charge axis (sidewall view). a — initial stage of explosions; b — full-time registrations. 


Each photodetector unit (Si-Ge photodiode + preamplifier) was calibrated by the 


use of standard light source in accordance with commonly used procedure. 


Specific features of emission related to the HE explosions can be elucidated by the 
consideration of outputs of Si-channel of the photodetector placed across the axis of a 
charge (see fig.2.19). As seen all the emission traces exhibit regular trends, namely, the 
presence of two amplitude-time stages. The duration of initial stage (fig.2.19a) is about 
200 us after the initiation of a charge. During this stage the intensity of emission (radiation 
flux) changes significantly, 1.e. sharply increasing up to some maximum value and 
afterwards smoothly decreasing. The second (later) stage (fig.2.195) also can be 
characterized by the appearance of a maximum (at ¢= 5 — 8 ms). However, the value of this 
maximum seems to be one-two orders of magnitude lower than that at the initial stage of 


explosion. 


The intensity and dynamics of variation of the measured radiation flax depend 
heavily on the composition of a given HE as well as on the orientation of the 
photodetector. The most pronounced differences are observed between the cases of TNT 
and RDX explosions. The parameters of emission related to the burst of TG 50/50 and 
RDX(p.) exhibit intermediate behaviors. 


These peculiarities reflect both qualitative and quantitative variations of explosion 
emissivity properties in the set RDX > RDX(p.) ~ TG 50/50 > TNT. As seen from the 
table 2.7 this sequence is characterized by the monotonic increase of carbon content in 
detonation products. On the strength of these data one can interpret the observed 
photodetector’s outputs in the following way. Initial sharply increasing emission seems to 
be associated with the air heated by shock wave. However, since short time interval (of 
order 10 — 20 us) after initiation the optical thickness of the shock falls and expanded 
detonation products starts to be responsible for the luminosity. In this connection it should 
be noted that the intensity of radiation of carbon (soot) particles is substantially higher than 
that of gaseous components of detonation products. From this point of view at the initial 
stage of explosion (Fig.2.19a) the radiation flux in the case of RDX burst reflects the 
attenuation of shock wave and expansion of gaseous components of detonation products. 
Considerable increase of soot content, as in the case of TNT explosion, results in the fact 
that the contribution of soot particles into radiation intensity becomes dominant at least at 
40 — 50 us after initiation of the charge. Among other things the presence of carbon as a 
fuel is a likely reason for after-burning phenomena at the later stage of explosion. These 
phenomena result in increase of the observed radiation intensity at the later stage of 


explosions (Fig.2.195). 


"0,1 1,0 10 &ms 


Fig.2.20 Temperature-time histories in the course of longitudinal (a) and lateral (b) 
expansion of HE detonation products. 

The dynamics of after-burning processes can be clearly illustrated by the 
temperature-time histories calculated by the model of black body radiation. As it is seen 
from fig.2.20 at the first 800 — 900 us temperature of detonation products falls down to 
values about 7 = 1200 — 1500 K. Further one can easily observe relatively smooth 
temperature rise to be associated with the burning phenomenon. The intensity of the after- 
burning process (i.e. heating rate and maximum temperature value) depends significantly 
from the soot content. Maximum temperature levels are achieved in the case of TNT and 
TG 50/50 explosions, while in the case of RDX burst the after-burning effects seem to be 
negligibly small. An important thing is that the after-burning process is spatially non- 
uniformed. The longitudinal expansion (fig.2.20a@) results in more intensive combustion 
(i.e. combustion starts earlier and temperature seems to be higher) than that lateral 


expansion (fig.2.204). This difference can be explained in the frame of the model, in which 


the rate of the reaction of oxidizing of detonation products is controlled by the rate of 
turbulent mixing of EP with surrounding air. In the scheme of relative displacement of the 
charge and detonator the initial velocity of expanding products in longitudinal direction is 


2-2.5 times higher than in lateral direction. 


Note, that the up to date known spectrum of different HE by nature and 
composition make possible to investigate the processes of afterburning in mixtures with 
another than soot reactive components. An example is the work [2.37], where specific 
features of oxidizing of the secondary products of HMX with Mg and Polytetrafluorethylen 


additives were studied on the basis of qualitative pattern of emission. 


Performed tests show that the photodetector based on the double-wavelength Si-Ge 
photodiode presents a universal mobile tool for registration of emission and temperature 
measurements in expanding clouds of HE detonation products. Due to the simple design, it 
is especially suitable for multi-channel registration in both outdoors and laboratory 
conditions. The use of the PD system makes possible to compare the characteristics of 
emission in different kinds of explosives. Application of the double-wavelength 
photodetectors to the burst of TNT- and RDX-based explosives enables to reveal 
peculiarities of after-burning phenomena associated with turbulent mixing of detonation 


products with surrounding air. 


2.12. Shock waves from non-spherical charges 


The HE charges are the spheres or hemispheres only in specially prepared 
experiments. More often the form of the charge is the cylinder of length Z and diameter D, 
i.e. cylinder with ratio L/D > 1. Generally, the form of the shock wave generated by 
explosion of cylindrical HE charge is not spherical at rather far distance from the source of 
explosion. The nonspherical shape of a wave means the inconstancy of overpressure at the 
shock front. The nonspherical shape of a wave and the inconstancy of pressure in different 
directions is kept till the overpressure AP ~ 0,014 MPa. For the less intensive waves the 
nonspherical wave shape and inconstancy of overpressure in different directions are not so 


essential and can be disregarded. For the long cylindrical charges with L/D > 1 it is 


possible to use the following expression for calculating the orientation diagram of a 


pressure field 


y = [2,0467 — 0,1146 x + (0,1285 — 0,0342 x) cos (8) + (0,0621 — 0,3280 x) cos (20) 
+(-0,0029 + 0,0304 x) cos (30) + (-0,1532 — 0,998 x) cos (40)] + [-2,1617 + 0,1422 x + 
0,2079 + 0,1161 x) cos (0) + (-0,4178 + 0,3686 x) cos (20) + (-0,1372 + 0,0648 x) cos 


(3e) + (-0,3484+ 0,1191 x) cos (40)] A + [0,4366 + 0,0418 x) + (0,0138 + 0,0938 x) cos 
(8) + (0,1178 + 0,1451 x) cos (20) + (0,2556 — 0,043 x) cos (30) + (0,3121 + 0,161 x ) cos 
(48)] a7. 


For the short cylindrical charges at L/D < | the following expression for calculating 


the orientation diagram of a pressure field is offered 


y = [2,0467 — 0,1753 x + (0,1285 + 0,072 x) cos (8) + (0,0621 — 0,2503 x) cos (28) 
+ (—0,0029 + 0,0079 x) cos (38) — 0,1534 cos (40)] + [-2,1616 + 0,0464 x + (—0,2079 — 
0,2174 x) cos (8) + (-0,4178 + 0,3426 x) cos (20) + (0,1372 — 0,1171 x) cos (30) + (-— 
0,3484 — 0,3449 x) cos (40)] A + [0.4366 + 0,0059 x + (0,0138 — 0,0006 x) cos (8) + 
(0,1178 — 0,2695 x) cos (20) + (0,2556 + 0,2072 x) cos (38) + (0,3123 — 0,2160 x) cos 
(40)] 07. 


In these formulas 

x — In(Z/D); 

i. — In (0,0348 R/G°”?): 
P3=8- 10° e, bar; 


8 — azimuth (horizontal) angle. 


Calculation configuration is shown in the fig 2.21. 


Fig. 2.21. The elongated cylindrical charge 


In case of ground explosion it is necessary to double the weight of a charge to 


calculate A. 


For understanding the opportunities of modeling of various types of charges 
explosions of the nonspherical form and arbitrary oriented relatively reflecting surface (as 
a rule, the surfaces of ground) the data on the measurements which have been carried out 
in [2.19] are rather useful. For the cylindrical HE charge with the length to diameter ratios 
from 1:1 up to 5:1, fixed on a cylindrical surface on ground, in directions on an azimuth 0°, 
22,5", 45°, 67,5°, 90°, 112,5°, 135°, 157,5°, 180° the overpressures were determined. The 
azimuth 180° corresponds to a direction of the blast wave propagation, and the azimuth 0° 
is the opposite direction. The azimuth 90° corresponds the perpendicular direction relative 


the blast wave. 


Along with the basic charge position, when the cylinder lies on the ground, the 
additional measurements were carried out for the inclined charge in 45° and 60° relative to 
ground. Explosion of the cylinder vertically standing on the ground also is considered. The 
results of measurements demonstrated, that the distinctions in AP are kept for arbitrary - 


oriented charges only at AP > 0,1 MPa, ie. at distances R/G'? < 5 m/kg"” 


. Otherwise, the 
distinction of AP falls within the limits of the accuracy of measurements. It permits to not 
consider the azimuth direction of wave propagation and to not take into account the 


orientation of the charge axis relative the earth surface at R/G'? > 5 m/kg'”. 
2.13. Parameters of shock waves at underwater explosion 


The reference homogeneous systems for two-phase fluids depending on a mass 
(volumetric) content of the condensed phase are either gas, or degassed liquid. Therefore, it 


is worth to list the dependences of the basic parameters for underwater explosion. On the 


basis of generalization of the information from [2.25 ... 2.33] the amplitude of the blast 


wave P, duration of the phase of compression t, and the impulse of the phase of 


compression / at explosion of G kg of HE under water at the distance X m from the center 


of explosion will be 


P, =K,(@ ‘RY 
tT, =K,, (Gc ‘RY 
fake Rk) 


The empirical factors Ki1, Ki2, Ki3, a,8, and y are summarized in table 2.8. 


Table 2.8 The basic parameters for calculation of the blast wave at underwater explosion 


Reference |Bjorn [2.30, 2.33] |Cole [2.26] Poche Bebb [2.32] |Arons [2.29] 
[2.31] 

HE Tetryl RDX TNT TNT 

RIG”, 14... 60 3,2 ...210 |0,16...5,3 |0,46... 795 

m/kg"" 

Ky,Pa [506 - 10° 560-10? /411-10° [522 - 10° 

a 1,1 1312 1s 1,13 

Ki2, ms 87 — = 92,5 


Ky3, Pas. |5,9 5.15 6,27 — — — 


¥ 0,87 0,89 0,98 = = — 


For completeness we shall specify, that the measurements in [2.27, 2.28, 2.30, 2.31] 
were made for small charges, measurements in [2.26] - for the large charges, 
measurements in [2.29, 2.32] are executed at the far distances from the location of 


explosion. 
2.14. Estimation of the explosion parameters 


The preformed analysis together with the data [2.34,2.35], shows, that there is a 
significant scatter of experimental data for all blast wave parameters and for all types of 
HE explosion in gaseous atmosphere (air). The known formulas for calculating the 
overpressure have the accuracy not more than 10 % even for an ideal case of air explosion. 
Therefore the accuracy of calibrating measuring systems with the help of HE will never 
exceed of accuracy of measurements of blast waves parameters. In the near zone of 
explosion at AP/Pp > 2,3 and R/Ro < 35 expressions for the impulse of a wave in the form 
[2.1] can not be recommended for use. In connection with the marked discrepancy of 
measurements of the pressure and the impulse the backward recalculation of the energy of 


explosion by overpressure and impulse leads to the errors in 1,5... 2 times. 


Therefore, the use of the parameters of the blast waves generated by HE charges for 
power estimations of other types of explosions should be applied taking into account the 
possible errors and is suitable only for qualitative estimations. At the ground explosion the 
uncontrollable and significant loss of energy on crater formation can be observed, in this 
connection the energy equivalent of the ground explosion in relation to air one frequently 
falls below theoretically expected value [2.34]. The measurements of the ground explosion 


parameters potentially comprise an uncontrollable error because of the energy losses. 


For this reason the data on the blast waves parameters are less suitable for 


calibrating of the measuring equipment, than the data on the air explosion. It is necessary 


to prefer here the measurements of blast waves parameters for the above-ground explosion, 


because already at a specific reference heights A, = 0,44 m/kg”? : 


the energy losses due to 
cup formation are small, so they could be neglected. Use of the data on blast parameters 
performed in the frame of equivalent estimations leads to even larger mistakes, than the 
use of the data on explosion of condensed HE. It is predetermined by a wide scatter of the 
measured parameters of strong explosions, by their small accuracy, and by their small 


number [2.15]. 


In most publications only the measurements of blast parameters from detonation of 
HE charges in air is reported.. Actually, the dependence on surrounding air density and 
speed of a sound really exists, and the increase of speed of a sound in gas environment (at 
constant initial pressure), should be accompanied by the reduction of overpressure. In 
general, the expected wave attenuation is caused by the tendency to decrease the acoustic 


resistance by the growth of speed of a sound, under uniform pressure. 


The magnitude of acoustic resistance pcg = 2 depends on pressure by formula 2 = 
vPo/co Thus, having warmed up the gas, addition of lightweight gases (hydrogen, helium 
etc.) promotes the decrease of acoustic resistance of gaseous atmosphere. As a result the 
smaller pressure perturbations will penetrate into surrounding atmosphere than in the case 
of heavyweight gases. The further change of wave parameters will depend on speed of 
wave propagation and the time of arrival of a wave the fixed point. It is necessary to 
mention as practically important feature the expected distinction of pressure waves 
parameters in a hot and cold season. Really, the acoustic resistance of gas (air) depends on 
temperature as Q = Po(W/RT, yaa . Thus, as a result of the blast wave reflection from the 
boundary of a cold gas, the wave with larger amplitude penetrates into the cold gas. 
Nowadays on the basis of available measurements there is no opportunity to establish a 
degree of correction, which is must be implemented into the known expressions at the 
realization of expert estimations. Most likely, the expected scatter of blast waves 


parameters depending on the weather condition does not exceed 20 % error. 


Below we list a set of the formulas for estimations of HE explosion parameters: 


1. Air explosion: 


AP [MPa] = 1.710° exp{-7.5R°"*!+0.0156 for 0,1 m/kg’*’ < Rs = B/G? < 8 m/kg 
p g 


0,33, 


AP [MPa] = 810° exp{-10.7R2'} for Re > 8 mv/kg®; 
2 Ground explosion: 


AP [MPa] = 1.7L10° exp{—7.14R2?8$+.0.0156 for 0,1 m/kg®? <R« 8 m/kg?*®: 
p g g 


AP [MPa] = 8110° exp{-10.46R)"'} for Re > Snvikg””; 

3. Duration of the compression phase: 

t, /G°* [ms/kg°**] = 7; =0.06 exp {1 AR} \ for 0,1 m/kg**? < Re< 1,3 m/kg?”*; 
7, /G°? = 75=13.8— RO" for Re> 1.3 mkg”™, 


At transition from TNT to another HE it is necessary to take Gue = KrnrG instead 


of G, where the values Krnr are given in table 2.1. 
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CHAPTER 3 CRITERIA OF BLAST WAVES DAMAGE 


3.1 Interaction of the blast waves with targets 


When a blast wave meets to a target a partial reflection takes place (fig.3.1). The 
lateral and top surfaces of the target are subjected by the pressure behind the passing wave 
AP. Windward side of the target will be exposed by at least doubled pressure load due to 
the stop of the flow. Nearby the lateral and top surfaces of the target there will be a flow 
with significant pressure gradients due to the rarefaction waves. Because of these waves 
the pressure at the windward side will decrease down to the stagnation pressure of a 
gaseous flow. After attaining the backside the wave bends the target and generate the 
turbulent zones. After the blast wave moves away, the target is subjected only by the wind 


load due to the friction forces arising from the streaming flow around the body. 


a ee 


Fig.3.1 Scheme of blast wave interaction with a target. 1- target, ISWF — incident shock 
wave front, RSWF — reflected shock wave front. 
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Fig.3.2 Three cases of passing wave interaction with obstacle: a) wave propagation 
alongside the surface, b) normal reflection, c) streaming flow around the small body 


It is possible to distinguish three limiting cases of targets blast loading (fig. 3.2) 


1. Wave does not “feel” the target. So, the target is subjected to the effect of 
the pressure behind the incident wave. 

2. Wave reflects at the target. So, the target is subjected to the effect of the 
pressure behind the reflected wave. 

3. Size of a target is small in comparison with the size of a blast wave. This 
provides almost instantaneous release of the target from loading. The target 
is subjected only by wind loading, which is determined by gas flow 


breaking. 


The magnitude of wind loading QO depend on the gas density p, the gas velocity w, 


and the shape of a target: 
O=Cj0.5 pur 
The friction parameter Cz is found by experimental methods [3.4...3.7]. 
3.2 Dynamics response of the target to the blast loading. 


The simplest scheme of blast loading of the target of mass m is presented in the 
fig.3.3. The magnitude of the active load P is determined by the above described schemes. 


The response force R defines the resistance of the target to the displacement: 
P=R+mdx/dr’, 
where g= d’x/d?’ is the acceleration of the target. 
At the static loading and at low acceleration P=R. 


At minor acceleration P2R 


Fig.3.3 The scheme of the forces acting to a target 


The response force X depends on the displacement x and the possible variants are 
shown in the fig.3.4. The variant in the fig.3.4a corresponds to the ideal elastic case: the 
resistance is directly proportional to the displacement. The ideal plastic case is represented 
in fig.3.45: the resistance does not depend on the displacement. Usually, in practice, the 


mixed variant takes place (see fig.3.4c, d): elastic case at small shift and plastic one at 
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Fig.3.4 The main types of target response to the force. a) — ideal elastic, b) — ideal plastic, 
c, d) — mixed. 


large displacement. 


The notion of the dynamic loading factor DLF is often used to calculate the target 
response to the blast loading. Using the DLF the static loading can be found, which 
generate the same effect as the dynamic loading under study. Usually the question of 


destroying loads is discussed [3.1, 3.2]. 
Let the system behavior 1s described by equation: 
Mq+Kx=P(t),where q= d°x/dt and K — is elasticity coefficient. 


The profile of the loading is plotted in fig.3.5. 


Fig.3.5 Scheme of the short-time loading of the target. 


The value of DLF depends on the parameter /,/7 (fig.3.6). The period of natural 


oscillations of the target is 7=2(m/K)°”. 


Fig.3.6 Dependence of dynamic loading factor on the relative loading duration. 


Two limiting cases can be marked: 


e for low ¢,/T <<1 the DLF parameter is proportional to the loading duration, 
the so-called impulse loading; 


e for high 7/7’ > 1 the DLF parameter 2, the so-called step loading. 


The cases of impulse and stepped loading are illustrated in fig.3.7. 


Fig.3.7 Two types of loading: a — impulse, b — step 


The target response at the impulse loading is determined as } Pdt and weakly 


depends on the function P=P(t). In case of step loading the time of the action is not of 
importance and can be considered as infinite value. Of importance is the time of the 


loading growth ¢, —see fig 3.8. 


Fig.3.8 Scheme of loading with preliminary compression phase. 


The dependence of the dynamic loading factor DLF on the parameter /,/7 is shown 
in the fig.3.9. For 4,0 the DLF —2, as in the case of ideal stepped loading. In the case of 
long -time loading growth DLF 1, and one can say about residual type of influence at 


DLF=1. More detailed illustration of the loading type on DLF is presented in the fig.3.10. 
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Fig.3.9 Dependence of dynamic loading factor on relative duration of preliminary 
compression phase 
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Fig.3.10 Dependence of dynamic loading factor on relative duration of loading for 
different types of loading [3.6] 


Blast and shock waves generated by explosion processes produce the impulse and 
stepped loads. The compression waves with long phase of pressure rise are closer to the 


static load. 


The method of diagrams overpressure AP- impulse / is often used to estimate the 
effect of dynamic target response on its vulnerability [3.4, 3.19, 3.32]. The idea of such 
diagrams construction is based on the dependence of vulnerability on the duration of 
loading. The space below or above the line of defeat corresponds to the dangerous or not 


dangerous combination of pressure and impulse for the particular target (fig.3.11) 
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Fig.3.11 The alg — Impulse diagram for building damage. 1- threshold of the minor 
structural damage, 2 — threshold of the significant structural damage, 3 — the boundary of 
the full (total) damage. [3.6] ¢ - group 1, & - group 2, VW - group 3, Wi - group 4. 


The points of the group | are the two-story wood buildings. 
The points of the group 2 are the two-story stone buildings. 
The points of the group 3 are the one-story wood buildings. 
The points of the group 4 are the one-story brick buildings. 


The numbers near the points denote the relative costs (in %) of the damage. Thin 
lines correspond to the relation between the impulse and the amplitude of a blast wave for 


indicated TNT weight. 
3.3 Critical levels of overpressure. 


One of the most significant and prevail factor (at any explosion accident) of 
damage is the impulse of pressure (velocity) generated of the pressure wave of a shock 
wave. Blast waves are created in surrounding atmosphere by the fast energy release as a 


result of explosion of a mass of energy carrier of any type, high pressure vessels burst, the 


failure of the vessels containing overheated liquid, etc. The only parameter — overpressure 
AP is used at the moment almost in all guides, technical references and expert handbooks 
for assessing the damage. No attention is paid to the shape of the wave and the durations of 
compression and rarefaction phases. No attention is paid also to the impulse of pressure 


carried by the wave 
t 
I= [APar 
0 


Here AP=P));-Po=AP;, P); and Po are the pressures behind the wave and in the 
initial medium. Traditional approach means the correlation between expected damage to 
the blast wave amplitude. Thus, the particular level of damage (often subjectively) is 
attributed to the characteristic value of AP; . The values AP; are presented in table 3.1 by 
phenomenological observations [3.1, 3.2, 3.3]. Note satisfactory agreement between 


domestic and foreign data. 


Table 3.1 Critical overpressures 


Limit of sound pressure for 
industrial personnel 
Rare damages of large\l window broken per . LkPa 
windows people 
(137dB) 
ringing in the ears 


ringing in the ears Pus2khP A 
(160 dB) 
1 window broken per 1 man 


Dangerous __ for light 
se rreees 


Injuries window 
fragments 
Light damage of roof and 

walls 


10kPa Typical pressure for S4.0kfa 
glass window failure 


20kPa 5% probability of severe| Rooting out of small trees 
damage 

22kPa More than 10% 
windows are broken 


22k e Damage of 
boards 

28kPa Lower limit 
structural damage 


40...70kPa All large windows are|Rooting out of trees, death) 30... 50kPa 
broken, rare damage oflof flying birds 
window frame 
20...35kPa Contusion of uncovered 
animals 


45... 85 kPa Window frames — are|Big and small animals lung|40...60kPa 
broken, laceration 
partial damage o 
buildings 


70...140kPa Asbestos covers are 
broken, bending and 
blown off of aluminum 
covers. Failure of wood 
panels. 


S5kPa i. 4 Lower lethal limit 


100 kPa Damage of steel frames 
140kPa Partial demolition o 
walls 


140......210kPa_ |\Demolition of _ slag- 
concrete walls 


SOkPa Light damage 7) 
buildings 


Lower limit of severe 
damage 


Demolition of — brick 
walls 


Light damage 0 
equipment 


60... 50 kPa 
70...60kPa 


50... L00kPa 


100... 200 kPa 


630 kPa Deformation of loaded 
rail cars 


The discrepancy of quantitative characteristics for the same level of damage is 


350... ..490 kPa |\Complete demolition o 
buildings, turning over 
of rail cars, crushing o 


Demolition of _ steel 
panels 


approximately 2 times, which is not acceptable for expert assessment of the real accident. 


In fact, the discrepancy in explosion energy estimations can reach 5-10 times. This 
can prejudice the determination of the energetic potential of the accident. The difference in 
AP; values is conditioned by the subjective nature of the damage evaluation. Below, as an 
example, we consider the critical values of AP; corresponding to the damage of window 


glasses, which are basically used in expert evaluations. 


The existing data on the critical overpressures able to destroy the window glasses 
are contradictory and therefore must be used very carefully. The critical value of 
overpressure depends on the width of the glass, window area, type of glass fixing, 
character of blast loading (duration, front steepness, amplitude modulation), and many 
others no controllable of badly controllable parameters. Thus, the indication of any exact 
value of the critical damage level of AP; without the fixing of the discrepancy 5( AP; ) is 
incorrect. American Bureau of Standards proposes the modern method on _ the 
consequences of single HE detonation in air. According to the standard S-2-54 (1976) [3.2] 
for every level of the overpressure the number of window glasses per 1000 habitants is 


indicated. For the urban areas the reference data are following: 


Overpressure AP, MPa Number of broken windows per 1000 


people 


The plot of the probability P of window failure on the overpressure in the wave AP 
is presented in fig. 3.12. The additional ordinate axis indicates the number of window 
glasses per 1000 habitants. One can expect the failure of at least | window per 1000 people 
at the pressure perturbations much less than usually reported in typical handbooks. For 
example, most of the domestic documents declare, that there is no window demolition at 
AP < 0.001MPa. According to USA standard at AP = 0.001MPa 1% of all windows (or 
200 windows per 1000 people) will be broken. The complete 100% window demolition 


occurs for the waves with AP > 0.01MPa. 
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Fig.3.12 Dependence of the probability of window failure on the overpressure in the wave 


Presented above discussion demonstrates, that the knowledge of the critical 
overpressures (shown in the table 3.1 or referred in [3.1, 3.2, 3.3] is necessary, but not 


enough instrument for expert assessment. 


3.4 Method of Pressure-Impulse diagrams 


More exact representation about an arbitrary system response to the blast loading 
can be obtained from the analysis, which takes into consideration the finite time of 


pressure wave action [3.4...3.7]. 


The time-space scale of the action is described in such a manner that both the 


amplitude of the wave and static /,, and dynamic J,;, impulses are taking into account. 
t 
I, = | APat 
0 
t 


i [ ewrdt: 
0 


Here pu’ is the doubled dynamic head, y and w are the air density and the velocity 
accordingly. The pressure amplitude and the static impulse determine the expected 
mechanical damage of the target, and the dynamic impulse and the amplitude of velocity 
head set the conditions of the possibility to blow off the object from initial location. The 
use of the impulse and the amplitude of the pressure wave allows to deeper understanding 
level of conditions and criteria of damages. The so-called danger number DN can be 
introduced for the quantitative hazard evaluation. For simple systems the danger number is 


evaluated as [3.4, 3.12] 
DN =(AP-AP*)(I-J;) 


As earlier AP; is the critical pressure amplitude necessary for demolition. 
Additionally the notion of the critical pressure impulse J; is introduced. In simplest case 


there exists hyperbolic dependence of DN in coordinates pressure-impulse. 


Let us consider a set of typical examples. Fig 3.13 illustrates the data on the human 


body damage by air blast waves [3.4, 3.5, 3.8, 3.11]. It is seen, that for all types of damage 


the decrease of wave duration (decrease of impulse) leads to a situation when any increase 
of the amplitude does not give rise to irreversible consequences. Qualitatively the same is 


the response of some mechanical targets to the blast wave. 


The typical pressure-impulse diagram for window breaking is shown in the fig. 3.14 


[3.14]. Another example of pressure-impulse diagram is presented in the fig. 3.15. 


AP,kPa 


0 07 14 21 28 3,5 42 4MPa-ms 


Fig.3.13 Diagram of human body damage by shock waves. 1- eardrums burst, 2- heavy 
contusion. 
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Fig 3.14 Pressure-impulse diagram for window breaking (line 1). Size of the window 
45x35 cm, width 3 mm. I-AP dependence for lkg of TNT (line2) 
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Fig. 3.15 Pressure impulse diagram for some constructions. 1 — turnover of truck, 2 — car 
cabine deformation, 3 — demolition panel constructions (antenna, panels, arcs) 


The problem of attaining the particular damage with the help of pressure — impulse 
diagrams is solved by the plotting corresponding functional relationship of the main 
explosion parameters in any of the diagrams in the fig. 3.13, 3.14, 3.15. Mutual position of 
the damage curve and the curve of wave characteristics determines the possibility and the 


level of expected effect. 


The experience of damage assessment caused by HE explosion nearby typical 
industrial or civilian objects showed, that the radius of the particular damage zone can be 


calculated using expression [3.14]: 


2 —0,17 
n,n (289) 


Here XK; — is the coefficient determining the level of damage i = I, H, UI, IV, G - 
TNT charge weight in kg. The value K; =3.8 corresponds to the zone of complete 
demolition. At Ky =9.6 the buildings are destroyed, roofs are blow off, the walls are 
damaged significantly. Ky =28 corresponds to the moderate damages. Finally, Ky =55 


indicate the zone of light (reparable) damages. The fraction of the broken windows in this 


zone is about 10%. Using the Rs; dependencies the corresponding pressure — impulse 
diagrams can be plotted by the following procedure: find the R using the known mass G 
and recalculate the coordinates at AP-/ diagram using the known dependencies from [3.4, 


3.5]. 


Fig. 3.16 Diagram for different levels of danger I-IV in combination with maximal levels of 
pressure-impulse for HE (line 1) 


Above-mentioned diagram is represented in the fig.3.16. The damage evaluation by 
overpressure is possible just for />>J;, where J; is the characteristic pressure impulse at 
the compression phase for the danger level under consideration. Lines I, II, HI, [IV 
characterize the response of the industrial and civil constructions to the blast loading. To 
evaluate the possibility to achieve the marked danger levels the functional relationship of 
pressure-impulse in the blast wave must be plotted at the same diagram. The use of the 
magnitude of critical overpressure (as usually done in practice) “silently” assumes that the 
pressure wave has significant duration, which is incorrect for the explosions of the limited 
fuels masses. It was shown in [3.4, 3.10] that the method of pressure —impulse diagrams 


gives the possibility to diminish the magnitudes of the safe distances without decrease of 


the level of the safety for the given objects. The expected profit in the safety distances 


estimation can be represented by the dynamic coefficient of explosion loads [3.2, 3.4]: 


: 3180 2 0.17 
aoe 


The dependence of the coefficient fon the weight of the charge (line 1) is presented 
in the fig.3.17. As seen from the figure, the profit in the safety distances due to the 
dynamic character of the blast loading can reach 1.2...1.4 times in the range of G under 


consideration. 
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Fig.3.17 The expected profit by the safety distances depending on HE weight. 


Recently, on the basis of analogous approach in France, the magnitudes of the 
safety distances for the small (up to 300-500 kg) HE charges were reconsidered. According 
to the new rules of HE storage in France the safe distance R’ =f R’, where R’ is the safe 


distance for the wave generated by strong explosion. According to the standards in France: 


G 
=1-0.5exp| -—— 
fi es =| 


The dependence of f (line 2 in the fig.3.17) on the charge weight additionally takes 


into account the effects of the fragments from explosion at small distances. 


It must be emphasized, that the real subdivision of the pressure-impulse plane into 
two parts (the internal — demolition zone, and the outer — stability zone) has no sharp 
boundary. With approaching to the danger threshold the probability of the particular 
danger increases. Finally the 100% probability of danger for the known pressure and 
impulse can be attained. This typical specific feature of the diagrams of danger can be 


represented by introducing the probit — function Pr; [3.8, 3.13]. 


The probability of hardly reparable damages of the industrial buildings (line IT in 
the fig.3.16) can be estimated by the equation: 


Pr, =5—-0.26InV, 


The factor V; is calculated by the overpressure AP; (in Pa) and the static impulse /; 


(in Pa’s) in the pressure wave. 


8.4 93 
17500 290 
y= Le 
AP. i 


The probability of non — reparable damages of buildings (line I in fig.3.16), which 


is characterized by their destruction, can be estimated by the equation: 
Pr, =5—0.22InV, 


In this case the factor /> is calculated as follows: 


7.4 11,3 
40000 460 
V, = +|— 
AP. I, 


3 3 


There also exist the formulas for calculating the probability of bioorganizms 


(including human body) blast damages. The probability of lungs contusion is estimated by: 


Pr. =5=5,74108,, 


where 
42 13 
S; = bes : 
; P i 


In last expression the non-dimensional pressure and impulse are determined by: 


here m — is the mass of bioobject in kg, pressure in Pa and impulse in Pa’s. The 
transition zone is shown in the fig.3.18. Some references report about the dependence of 
probability of human eardrums burst on the overpressure in the air wave: 


Pr, =-12.6+1.524In AP. 


Usually, the pressure wave acts on the live organisms not only by the shock 
increase of static pressure, but also by the dynamic impulse, which can blow off the human 
body with high velocity into the direction of wave propagation [3.34...3.36]. The 


probability of human body blow off is estimated using the value 


Pr, =5.0—2.44In S; , where 
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Fig.3.18 Diagram of human damage by blast waves. I — safe zone, II — damage, III - 
transition zone, | - lower threshold, 2— 1% lethality, 3 — 10% lethality, 4 — 50% lethality, 
5 — 90% lethality, 6 — (AP-I) for 1 kg TNT. 


Additional evaluation of the damage level at the displacement of live organism with 
finite velocity can be performed with the help of pressure-impulse diagram from [3.6, 3.8] 
reproduced in the fig. 3.19. Lines 1,2,3,4 in the fig. 3.19 corresponds to the following 


damage levels: 


e safe impact 

¢ critical (impact threshold — lethal case is possible) 
e 50% lethality 

e¢ 100% lethality 
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Fig.3.19 Diagram of human damage due to dynamic head. 


Analysis of expressions for functions Pr; (excluding Pr,) shows their close 
connection with the impulse characteristics of the pressure wave. Any of the damage levels 
(typical for Pr;, Pro, Pr3, Prs) can be attained only for the finite value of compression 


phase impulse, i.e. for /5~0. 


3.5 Secondary effects 


In preceding chapters only two direct factors of danger owing to the blast waves 
were considered: static pressure jump, and sudden air movement by the compression wave. 


However, of importance are also the secondary effects accompanying the blast loading. 


One of the most significant secondary effects is the fragments flight from destroyed 
structural elements. Above analysis shows, that the window glasses have the lowest 
threshold of strength. The important practical information on the velocity of window 
fragments for different window areas and glass widths is contained in the reference [3.9]. 


The fig. 3.20 represents the data from [3.9]. 


The dependence of the maximal range (flying distance) of the glass fragments on 
the blast wave pressure is an important characteristic for hazard evaluation and risk 
assessment. The plot in bottom part of the fig.3.20 determines the expected range of 
fragment flight. It is seen, that for the blast overpressure AP = (2...20) kPa the fragment 
flying distance is about 10...30 m, and the velocity of the fragments is about 10...30 m/s. 
The mass of the fragments of glass does not exceed 100 g in accordance with the data from 


[3.9]. The probability of light injuries is evaluated by the probit function [3.6,3.8]: 


Pr, =-29.15+2.10InS,, where 


Here mo is the mass of the fragment in kg, and v is the velocity of the fragment in 


m/s. 


0 5 10 is 20 AP, kPa 


Fig 3.20 The velocity (a) and the maximal distance (b) for the flight of the fragments of 
glass windows caused by blast wave. The shaded area corresponds to the scatter of data 


from [3.9]. 


The function probit Pr; is related to the probability p of particular level of danger by 
the table 3.2. 


Table 3.2 Probit vs. probability [3.8, 3.13] 


‘Muon HuoHuMMNO 


0 267 2.93) 312.- 3;25° 338° SAD 352°. 359° -3.66 


10 3.72 3.77 3.82 3.86 3.92 3.96 4.01 4.05 4.08 4.12 


20 416 419 423 426 429 433 436 439 442 445 


30 448 450 453 456 459 461 464 467 469 4.72 


40 475 4.77 480 482. 485 4.87 4.90 492 4.95 4.97 


50 200° °3.03 -5:05" 508: 310: S12 SAS: SES) -3200 “323 


60 5.25. 5.28 531 533 5.36 539 541 544 547 5.50 


70 Doe. SOD SS 5.6L $64 3.67 S71 . D742 SAT S81 


80 5.84 588 592 595 5.99 604 608 613 618 623 


90 628 634 641 648 655 664 675 688 7.05 7.33 
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3.6 Effect of bio-object position on the damage caused by blast wave. 


When the blast wave propagates around the human body, the loading significantly 


depends on the mutual orientation of the wave front and body position [3.15]: 


e The blast wave propagation along the lying body into direction legs — head. 

e The facing of the blast wave with the standing body, or its propagation 
perpendicular lying body. 

e The facing of the blast wave with the standing or lying body nearby 


reflecting surface. 


The damage diagrams for the above noted cases are presented in the fig. 3.21a,5 in 
coordinates compression phase duration 7, - blast wave overpressure AP by [3.4, 3.5, 


3.15]. 


Fig 3.21 Critical levels of blast amplitude and duration for different human body 
orientations relative to the blast front. a)The blast wave propagation along the lying body 
into direction legs — head, b) The facing of the blast wave with the standing body, or its 
propagation perpendicular lying body, c) The facing of the blast wave with the standing or 
lying body nearby reflecting surface. Line | — for 1 kg charge, line 2 — for 1000 kg charge, 
lines 3...7 — survival levels, line 8 — thorax damage threshold. 


The cases where minimal blast pressure and impulse are required for the damage 
correspond to the positions of human target lying or standing nearby the reflected surface 
(fig.3.21c). However, the consideration of the shape and type of reflecting surface, the 
incidence angle and the distance of human target from reflecting surface is of great 
difficulty. 


Moreover, in the real conditions, the reflecting surfaces nearby the human target 
could be absent, which restricts the application the lethality (or survival) curves 
applications for the case presented in the fig. 3.21c. In case of blast wave propagation in 
free space the maximal effect on human body occurs at perpendicular facing of the incident 
shock (fig.3.214). This position is rather probable and is considered as the basic one at 


primary damage assessment at accidental explosions [3.4]. 


3.7 Blast wave effect for people 


Blast damage of human body. 


The blast damage of human body is conditioned by the existence of various internal 
organs with different density and the gas-filled cavities inside the corp. First of all of 
importance are the lungs. Every lung contains about 300 billion of air cavities — aveols. 


The thickness of their walls is 5=0.2 um, that is 370 time narrower than hair width [3.4]. 


The blast wave, faced to the human body, partially penetrates inside and propagates 
through viscous multiphase medium. The pressure wave meets the aveols and compress 
their gaseous filler. This results in their destruction and the burst of neighbor capillary 
[3.37]. Finally internal blood flow starts which causes lethality. Also dangerous is the gas 


inflow into blood at rarefaction phase (caisson disease analog). 


For the detonation of HE charges up to 20 kg the severe damage threshold takes 
place at blast wave amplitude AP>270 kPa. Such a wave is expected at 250 g HE 


detonation at the distance of outstretched arm. 


In recent years the attention of physicians are attracted to the blast damage of such 
internal organs as stomach and intestinal tract, resulting in internal blood flow, burst and 
perforation of walls. The listed injuries are latent and are difficult diagnosed, which causes 


heavy long-time physiological consequences. 
Blast damage of human extremities. 


The severe damage of extremities (arms, legs) occurs at the extremely powerful 
blast load exceeding AP>1.5 MPa. The extremities contain small amount of gas, so can be 
considered as incompressible. The damage takes place as traumatic amputation, therefore 


requiring larger blast loads than in the case of breathing organs. 


At the same time, during the expertise one must take into account the fast 


attenuation of the blast wave with distance. 
Blast damage of human ears 


The ear, a sensitive organ system that converts sound waves into nerve impulses, 
responds to a band of frequencies ranging from 20 Hz to 20,000 Hz. This remarkable organ 
can respond to energy levels as low as 10 “'’watt/m’, which causes the eardrum to deflect 
less than the diameter of a single hydrogen molecule [3.33]. Not being able to respond 
faithfully to pulses having periods less than 0.3 ms, it attempts to do so by making a single 


large excursion. It is this motion, which can cause injury to the ear. 


The human ear is divided into the external, middle, and inner ear. The external ear 
amplifies the overpressure of the sound wave by approximately 20 percent and detects the 
location of the source of sound. Rupture of the eardrum, or tympanic membrane, which 
separates the external ear from the middle ear, has captured most of the attention of 


clinicians although it is not the most severe type of ear injury. The eardrum and ossicles of 


the middle ear transfer acoustical energy from the external ear to the inner ear where 
mechanical energy is finally converted into the electrical energy of the nerve impulse. The 
middle ear is an impedance matching device as well as an amplification stage. The middle 
ear contains two dampers: the stapedus muscle and associated ligaments, which limit the 
vibration of the stapes when subjected to intense signals, and the tensor tympani muscle 
and its adjoining ligaments, which limit the vibration of the eardrum. The first damper is 
the most important. These dampers have a reflex time of approximately 0.005 to 0.01 s, 
which is longer than "fast" rising air blasts. The manner in which the malleus and incus are 
linked allows far more resistance to inward displacement than to outward displacement. 
However, if the eardrum ruptures after inward displacement during positive phase of 
loading of the blast wave, the malleus and incus are less likely to displace as far outward 
during the negative phase of loading of the blast wave as they would if the eardrum 
remained intact. In this case, eardrum rupture could be beneficial. The maximum 
overpressure and its rise time, however, control the characteristics of the negative phase 
and is, therefore, of prime importance. Rupture of the eardrum, thus, becomes a good 


measure of serious ear damage [3.4]. 


3.8 Expert assessment of blast damage. 


To perform the fast (urgent) expert estimations of damages from accidental 
explosions and terroristic acts it is necessary to have the simple and convenient tools based 


on the validated experimental and theoretical observations [3.15, 3.34...3.37]. 


One of the useful and simplest instruments is the nomogramm presented in the fig. 
3.22. This nomogramm in the simplest form allows determining the consequences of blast 
loading on unprotected human body at the explosion of up to 22 kg TNT charge in near 


zone. 


The nomogramm can be used in the following manner. 


Let we have the explosive object of 4 kg of TNT. The distance to the object is 5 m. 
According to the nomogramm, the expected blast wave will have the amplitude AP=75 


kPa. One can expect the eardrums burst with the probability 25% [3.34...3.37]. 


The same HE charge placed at the distance 1.5m will generate the wave with the 


amplitude AP=1 MPa, that guarantees 15% lethality, complete deafness, and the total burst 
of eardrums [3.37]. 


Fig. 3.22 Nomogramm for determining the blast damage of bioobject depending on the 
distance and HE weight. Horizontal lines: 1 — lower threshold of blast damage for ears, 2 
— 50 % damage of ears, 3 - lower threshold of blast damage for lungs, 4 — 50 % damage of 
lungs, 5 — 1% lethality, 6 — 10% lethality, 7 — 50% lethality, 8 — 90% lethality, 9 — 99% 
lethality. 


3.9 Assessment of shock wave effect on the bio-objects at blasting operations. 


The principles of ecological safety [3.16] require the development of the methods 
for the prediction of the shock waves (SW) effect on the bio-objects at technological 


blasting operations. The safety distances for blasting operations are determined by 


empirical formulas [3.17, 3.18]: 


1 =kG* (3.1) 


where 7, is the safety distance in m, G is the energetic characteristic of explosion in 
kg of TNT equivalent, v and & are empirical coefficients, characterizing the sensitivity of 


the object to the blast loading. 


The formula (3.1) must be used carefully, because the safety distance has no 
agreed-upon notion. According to the safety rules [3.17] at vy =1/3 and & = 15 in (3.1) the 
safety distance for the technical staff at the ground explosions can be obtained. However, 
these distances are differing from the safety distances accepted in other countries. For 
example, the requirements of USA Ministry of Defense (DOD 5154.45) of June 23 1980 
[3.19] concerning the technical operations with explosives declare, that at the distances, 
which can be calculated by (3.1) with vy =1/3 and & = 16...20 the expected blast wave 
effect with high confidence level guarantees the absence of lethal cases, however the 
injuries connected with walls and windows demolition are possible. Ambiguous definitions 
of the safety criteria give rise to the necessity of serious analysis of the existing norms of 
safety and the development of the methods for safety distances assessment at blasting 


operations in different conditions [3.20]. 
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Fig.3.23 Nomogramm [3.16] for determining the SW effect on human body from 
condensed HE explosion. 1 — 10% of people are temporary deaf , 2- safe distance 
according norms [3.22], 3 —norms [3.17]. G HE weight in kg, R distance from the center 
of explosion in m. 


Analogous problems arise in the determining the safety distances for the 
underwater blasting operations [3.21]. The expected effect of underwater shock waves on 


the fish is evaluated using the formulas similar to (3.1). 


For underwater explosion and for the most sensitive kind of fish [3.18] 7. is 
determined from (3.1) at vy =1/2 and k = 45. The dependencies of (3.1) type are valid for 
the point charges in the narrow range of the explosion energy and represent the particular 
cases of the bioobject response to the blast loading. To develop the general methods of the 
damage predicting for bioobjects the data on the blast loading in pressure — impulse 
coordinates AP- / are required. The curve | in fig. 3.23 corresponds to the probability of 


the 10% of people become temporary deaf at normally reflected SW. 


The nomogramm of damage of sensitive kinds of fish by blast waves is presented in 
the fig. 3.24. This nomogramm is taken from [3.16, 3.21] and was obtained on the basis of 
the experimental data on safe and the lethal distances at the explosion. According to the 


definition [3.21] 7, denotes the boundary of the zone, inside which the organism is killed, 


and 7 denotes the boundary of the zone, outside which the organism is free of damages. 
The curve | in the fig.3.24 corresponds to the 50% lethality threshold, and the curve 2 
characterizes the situation, where the SW effect results into the temporal lost of spatial 
orientation for the fish and the burst of blood capillary in 50% cases. The comparison 
between the damage curves for the fish and for the human eardrums judges about the 
generality of the damage laws for the bioobjects. In each case the damage curves are 
hyperbolas with two asymptotes corresponding to the critical pressure P., at quasi-static 


loading regime and to the critical impulse /,,. at impulse loading regime. 


The curves | and 2 in the fig.3.23 and in the fig 3.24 clearly demonstrate the 
existence of impulse, transition and quasi-static regions of loading. Also emphasized is the 
probabilistic character of the bioobject damage, which is conditioned by the individual 


characteristics of the organism affected by SW. 


Fig.3.24 Nomogramm of damage of sensitive kinds of fish by blast waves [3.16]. 1 - 50% 
lethality threshold, 2 - temporal lost of spatial orientation for the fish and the burst of 
blood capillary in 50% cases, 3 - measured overpressure and impulse after explosion of 
cylindrical charge (v=2) with linear density 16 g/m in the bubbly medium with volume void 
fraction B= 2%. For spherical charges (v=3). 


The damage diagrams in coordinates AP- / are the characteristics of an object and 


does not depend on the type of condensed HE or another energy sources able to generate 


SW. Moreover, the shape of SW (plane, cylindrical, spherical) is unimportant for 


determination of the damage zone. 


Mathematically, the hyperbola equation for the curves in the fig.3.23 and 3.24 can 
be written as [3.23,3.24]: 


DN =(AP-P,,)(I-I,,) = const (3.2) 


where AP and / are the overpressure and impulse of dynamic load. 


The equation (3.2) determines the third critical parameter DN , which characterizes 


the object sensitivity to dynamics loading. 


Extrapolating the experimental dependences (fig.3.23 and 3.24) into the region of 
impulse and quasi-static loading, the critical damage parameters were obtained in [3.16] 


and represented in table 3.3. 


Sensitive types | Lethal 
of fish Lost of spatial 


orientation for 
the fish and the 
burst of blood 
capilla 
Temporal deaf 
Level of 
accepted 
influence 


0,12 2,3:10° 
800 0,07 75 


The analysis of the data obtained in [3.20, 3.25] as well as the AP- / diagrams 
presented in [3.16] gives rise to the conclusion, that the critical parameters of any two lines 


of bioobject damage in the first approximation will be related as follows: 


Fon fon PM, _ Fon fon (3:3) 


Let us use the AP- / diagrams in order to determine the threshold curve of accepted 
shock-wave influence on the human body at blasting operations in urban areas. According 
to the requirements [3.17, 3.22] the level of sound wave pressure, measured at frequencies 
sensible for man ear, should not exceed 125 dB. The experiments with G=0.01...1.00 kg 
HE charges show, that level of sound pressure is 27 dB lower (in average) than the 


maximal overpressure at the incident SW front. 


Therefore, for the blast loading the accepted level of sound wave pressure 125 dB 
corresponds to the SW overpressure AP =800 Pa. Using the latter value as quasi-static 
asymptote of blast loading P.,. and substituting into equation (3.3), the curve 2 in the 
fig.3.23 was obtained in [3.16]. This curve defines the region of SW parameters, in which 


the effect on the human body is under the requirements of technical norms [3.22]. 


The critical parameters, characterizing the curve 2 (table 3.3), can be used for 
determination of the safety distances at technological blasting operations in industrial, 


urban and other conditions. 


The prediction of SW effect on the bioobjects with known DN , P-, Jc, require the 
knowledge of the dependencies AP(R:+) and /(R:). 


The dependencies of blast waves parameters and the distances to the HE charge are 
plotted in fig. 3.23 as a grid r-G. The intersection of the gridlines AP=/() with curves 1...3 
define the distances corresponding to the different SW damages of human ear organs. 
Thus, for G=lkg the magnitude of r, = 77m by [3.22] and r, = 15m by [3.17]. At the 
distance 47 m one can expect that 10 % of people will be temporary deaf. The comparison 
of above mentioned data judges the significant underestimation of 7. in [3.17] and 
inapplicability of these norms at technological blasting operations in industrial, urban and 


other conditions. 


Fig. 3.24 presents the graphics solution [3.16] of the problem of fish damage 
regions at underwater explosion. The dependencies AP=//J) for various G are 
approximated by the straight lines, which permits to determine the necessary distances 
using the coordinates r-G for different symmetry. An example of determining 7, and rx is 
demonstrated in the fig 3.24 for the explosion of 1 kg TNT charge. The demolition cords 
are often used as intermediate detonator at underwater blasting operations. For the 
cylindrical charge with linear density G=16g/m the distance r, will be 49 m, which is close 
to the r, corresponding to Ikg HE charge. The doubling of the specific energy of 
cylindrical charge will increase the 7, value up to 86 m. At the same time, the magnitude r, 


=94 m can be attained just for 5 kg HE charge. 


The effectiveness of using gas-contained media for explosion effect localization is 
clearly demonstrated by the curve 3, calculated by the wave parameters at explosion of 
cylindrical charge with linear density G=16g/m in water containing 2% of the air (bubbly 
medium). The intersection of lines 3 and 2 corresponds to 7, =0.75 m. the addition of 2 % 


vol. of air bubbles permits to decrease the fish damage radius in 50 times. 


Specific feature of gas-containing media, used for shock waves damping, is the 
shift of AP=f() curves into direction of quasi-static loading, that can be explained by 
lower SW impulse suppression in comparison with pressure. Note, that due to probabilistic 
character of SW damage for bio-object, the notion of safety itself becomes relative one. It 


is more correctly to say about permission distances but for safety distances. 


The permission risk level for particular case must be determined taking into account 
the economical and social factors, providing the environmental protection and, at the same 
time, allowing the progressive technologies. Based on above mentioned facts, the 
conclusion is made in [3.16] that the present requirements [3.17,3.18,3.21] of ecological 
safety in technologies, using pulse energy sources in offshore water, correspond to the risk 
level, represented by curve 2 in the fig. 3.24. The use of pulse energy sources in special 
conditions, when it is required to exclude additional effect of dynamic loading on the fish, 
the critical pressure (curve 2) must be diminished in 2...3 times, changing all other critical 


parameters according formula (3.3). 


3.10 Estimation of the wave critical parameters for window glasses. 


Important part of the requirements for explosions use the predicting of the shock 
waves effect on the objects located nearby the place of blasting operations [3.38]. The 
complex character of the object response to the dynamic loading does not allow to connect 
unilaterally the demolition to the either overpressure or the impulse of the blast wave 
[3.26...3.28]. As already mentioned above, one of the methods suitable for dynamic 
loading is the construction of the AP-/ diagrams. The curve of the object damage (or 


demolition) under dynamic loading has the hyperbolic dependence (3.2). 


The AP-/ diagrams, corresponding to the window destruction, are taken from [3.28] 
and presented in the fig.3.25. The curves of damage probability (1, 50, and 100%) for 
different window glass sizes are plotted on the basis of processing of large number of 
experimental data [3.14, 3.23, 3.27, 3.38, 3.39]. The curves in diagram demonstrate the 
existence of two regions, corresponding to the quasi-static and impulse loadings, as well as 
the probabilistic character of glasses destruction related to the individual properties and 


ways of fixing of the windows. 
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Fig.3.25 Damage diagram for the window demolition and the curves of pressure and 
impulse change for HE detonations. Destruction: 1-1%; 2-50%; 3-100%, [3.23]; 4-50% 
[3.39] for glass size 450x550x2.1 mm. 


The curve 4 corresponds to 50% destruction probability of glass [3.23,3.39] with 
fixed size. The critical parameters of damage for this curve are following: P.,= 5kPa, [= 
17.9 Pa s, DN = 5.68 10° Pa’ s. The knowledge of listed critical parameters for the glass of 


specified size is sufficient to draw the AP-/ damage diagram. 


To predict the damage of glass with arbitrary geometrical sizes it is necessary to 
find the effect of the width a, the length a, and the thickness 6 on the critical parameters 
P.,and I, . 


The analysis of presented AP-/ diagrams for the glass and other objects [3.4] shows 
that the relationship between the critical parameters in first approach can be written as 
[3.38]: 


L,.P..DN* =const (3.4) 


cr” CF 


The equation (3.4) allows constructing the AP-/J diagram and determining the 
critical damage parameters for any glass, if the critical parameters for particular glass are 
known. Moreover, the processing of the diagrams from [3.4] made in [3.38] yields the 


proportionality law for the critical values of the glass damage: 


ch = Ae (3.5) 


Using the (3.5) one can determine all set of the critical parameters (including DN 


by (3.4)) if the change of one particular parameter is known. 


The critical blast wave pulse loadings for glass destruction are not realized in 
practice. Therefore, the existing data on dynamic loads are, in fact, referred to quasi-static 
loading. The data on critical impulses presented in [3.29] are not actually critical ones in 


accordance to the definition (3.2). 


The critical and permission loads at quasi-static loading are analyzed in [3.38]. 


Experimental data on the glass destruction at gaseous explosion permit to obtain the 


minimal value of P.,, for 50% damage for types of glasses corresponding to the State 
Standard. Thus, for the glass thickness 6=3mm we have Pe mnin(50%)=2.2 kPa, and for 
O=4mm the Pe min(50%)=4.1 kPa [3.38]. Moreover, the State Standard declares, that the 
permission pressure should not exceed 850 Pa. According to building rules, the maximal 
wind pressure on the glass (which corresponds to the quasi-static loading) varies depending 


on the region from 270 Pa (Moscow region) up to 1 kPa (Arctic Ocean), for Kiev — 350 Pa. 


Permission overpressure at special regions was 500 Pa [3.27]. Thus, at blasting 
operations in urban areas, when the types of window glasses are unknown, on must refer to 
permission pressures Pe presented in technical documentation. For special regions one 
can recommend the critical values Prem = SOO0Pa, or Poem = 850Pa if the glasses are 


constructed according the State Standard. 


If the dependences AP/(R:) and /(R:) for energy source, which generates the blast 
wave, are known, then using the relationship between geometrical glass characteristics and 
critical damage parameters [3.38], one can determine the regions of different damage of 
windows as well as the safe distances. Due to the complex wave diversion in case of 
industrial plants the dependencies AP(R:) and /(R+) differ from the analogous dependencies 
obtained analytically for one-dimensional wave propagation. In this case it is more 
convenient to use directly AP-/ damage diagrams together with the experimental 


dependencies AP(R«) and /(R:). 


The dependencies AP(R+) and /(R«) obtained for charge explosion in air or in other 
media, used for explosion localization (dusty gases, foams, etc.) at any one-dimensional 
symmetry of space, permit to determine the effect of the blast wave on the window glasses 
at known parameters P-,, J--, DN. The curves AP(R«) and /(Rx) calculated by [3.30] for 
spherical charge detonation are plotted in the fig. 3.25. The dashed lines correspond to the 
equal distances from energy source to the object. The intersection of the dashed lines with 
curves corresponding to the different damage probabilities determines the distances where 


the indicated damage level is realized. 
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Fig. 3.26 Dependence of the safe distance rg (1) and factor x (2) on the energetic 
characteristic of explosion for glass (450x350x2.1) at shock wave parameters by [3.14]. 


For weak shock waves in air (AP< 10 kPa) the empirical dependencies from [3.30] 


are approximated by in [3.38] by relationships: 
AP=AR™", I=BG°?R", R=rG?*, 


where 1.13, B=200, A=109. Here G is the TNT equivalent in kg, 7 is the distance 
to the charge in m, AP is the overpressure in kPa, / is the pressure impulse in Pa s. Then 


the safe distance 7, can be expressed as a function of the charge weight G, namely: 
ro=k(G) G". 


Here k(G) is the weak polynomial function on G. 


Fig.3.27 Dependence of the standard glass (6=3mm) destruction on the weight of HE 
charge: 1 — 50% damage, 2— 1% damage, 3 — safe distance for passing wave from 
spherical charge, 4 — safe distance at ground explosion taking into account shock wave 
reflection. 


As example, the dependencies 7(50%) and 7(G) for glass and shock wave source 
are plotted in the fig. 3.26 [3.14], where good agreement between calculated and measured 
values is pointed out. Dependence 7(G) has the limits: 7=2/3 at G0 and 7=1/3 at Go. 
Fig. 3.26 demonstrates, that at 7=0.5 a transition zone is realized for the charges weights 
1...10 kg, which are commonly used in metal-working by explosions. It explains the fact 


that often used empirical formula r.=k G'” [3.17] predicts the safe distances rather good. 


The dependencies 7(G) for standard glasses (thickness 6=3mm) and for various 
damage levels are plotted in the fig. 3.27. The calculation parameters for damage diagrams 


are listed in the table 3.4. 


The choice of critical parameters for reflected semispherical wave was done by 
taking into account the doubling of shock wave overpressure and wave amplification in the 


course of ground explosion. 


Table 3.4 Critical parameters of damage diagram 


Risk level 


2200 5,68'10" 


1% 
io passing 


spherical shock 
wave 


fe reflected | 337 1,21 259 151 
semi-spherical 
shock wave 


3.11 Accident prevention at blasting operations 


The information presented in chapters 2 and 3 shows that the dangerous 
consequences of energy release by explosion sources of various types can be parried just 
by passive and active action on the blast wave with the aim to decrease the overpressure 
and impulse. This difficult task can be resolved either by direct action on the shock wave 


or by effective absorption of the explosion products energy. 


To develop the regulations of ecological safety at technological blasting operations 


it is necessary to know: 


e the critical parameters of the dynamic loads resulting in damages of the 
objects to be protected; 

e the relationship between the shock waves parameters and energetic and 
geometrical characteristics of the HE charge in the particular conditions of 
the blasting operations; 

e the correlation dependence between the parameters of shock wave loading 
at explosion and the existing regulations of practical pulse sources of 


energy. 


The promising way of the solution of the explosion technology safety problem is 
the development of the explosion localization methods on the basis of gas-contained media 


(foams, bubble structures, dusty gases etc.). 


In accordance with the concept [3.16] the approximate consequence of explosion 


technology realization can be presented as a scheme in the fig.3.28, suggested by BI. 


Palamarchuk. The decision — making algorithm of accident prevention at blasting 


operations is shown in the fig. 3.29. 


Main operations Main content 


Task project description 

Goals, motivation, solutions 
Specific features planning of neighbor objects 
of blasting operations location and safety vulnerable objects 

size of damage zones 

choice of detonation location 


Object object plan 
disposition of bore-pits 
disposition of pre-weaken places 


Explosion technology method of destruction 
type of material 
type of bore-pits 
explosion circuit 
protection 
special methods of protection 
protection by local measures 
Scenario of operations coordination scheme and responsibility 
complete timing scheme of explosion 
sequence 


Fig. 3.28 Scheme of the explosion technology realization. 


Shock waves parameters Critical parameters of 
from HE detonation blast loads 


Possible methods of Comparison of the main 
blast loads suppression methods of suppression 


Fig.3.29 Block-scheme of decision — making algorithm 
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CHAPTER 4. ELEMENTARY WAYS OF EXPLOSIVE LOADINGS 
REDUCTION 


At early stages of the explosive technologies implementation it seemed, that the 
simple placement (contact or non-contact) of HE charges in various protecting envelopes 
will result in lower levels of explosive loadings of an acceptable level [4.1]. However 
further practical work has revealed hopelessness and significant inconveniences of the 


given technical decisions. Let's outline some examples. 
4.1. Contact placement of charges in a protective envelops 


The numerous measurements, made for an estimation of ammunition efficiency, 
establish the dependence of air shock wave parameters on relative weight of HE charge (G) 
and weight of closely placed armoured envelope (W). Dependences AP =f (R, W/G), 1=R 
(W/G) were experimentally constructed. It occurred, that the detonation of HE charge is 


equivalent to explosion of the non-protected charge with weight 
Gz =|0.2+0.8/(1+W/G)|G 


Thus by [4.2] for W/G < 1 the efficiency of demolition action of the armoured 
charges is higher, than non-armoured charges, because of the absence of dispersion of an 


external HE layer (see chapter 2). Measurements in [4.1, 4.2] show, that the contact 


placement of HE charge in an armoured protective envelope should not be considered at all 
as the method of shock wave attenuating, except for the obvious case, when the 


environment does not collapse under the action of explosion. 


4.2. Non-contact placement of charges in a protective envelops 


More effective, than just above considered, is the way of protection from the shock 
wave with the help of the continuous screen which is not placed into the contact with HE 
charge [4.3]. The typical scheme of installation of the non-contact screen and HE charge is 


shown in the fig. 4.1. 


Fig. 4.1. Installation of the non-contact screen 


The screen | represents an envelope as a lateral surface of the cylinder or prism, 
inside which the HE charge 2 is placed at its longitudinal axis. The envelope has a height 7 
an internal diameter D. The problem of explosive wave propagation around of the 
protective envelope with the height H = 100 cm (or 400 cm) and the radius & = 20, 25, 35, 
50, 75, 150 and 250 cm was considered in [4.3]. The explosion of spherical TNT charges 
with the radius Ryz = | cm (weight 6,83 g) and the radius Ry¢ = 8 cm (weight 3497 g) was 
studied (fig. 4.2). 
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Fig. 4.2. Possible variants of the screens placement 


In preliminary calculations the optimum height of the screen by the magnitude of 
energy carried away by the shock wave was determined. In a fig. 4.3 it is possible to see 
the reduction of the energy carried away by the shock wave, with the increase of the height 
of the limiting cylindrical envelope. The time /o corresponds to the moment of wave arrival 
to the open section of an envelope. The obtained results for the screens of 120 and 200 cm 
height are close to each other, and the optimum height of the envelope for <1 kg TNT 
charges is 100 cm. The shorter screens actually do not reduce the level of energy of the 


wave and do not suppress the explosive loading. 


The levels of the maximal explosive pressure P,, = AP + Po and the blast wave 
amplitude AP at initial pressure Py = 10° Pa for various variants are shown in the fig. 4.4 
and 4.5. The reference lines 1 in both diagrams correspond to the ground explosion of TNT 


charges with the radii Rug = 1 cm and Ryg = 8 cm accordingly. 


All variety of the variants considered does not confirm the principal difference in 


expected level of the blast loading suppression. Therefore, the dependence of factor K of 


the blast wave amplitude attenuation on the relative distance R/Ryg is presented in the fig. 
4.6. It is seen, that non-contact protective screen cannot essentially suppress the blast 
loading. Thus the way of suppression of blast loading by the non-contact screen cannot be 
recognized effective on the basis of the presented results. The conclusions from [4.3] 


completely correlate with [4.2], where the other types of non-contact screens are discussed. 
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Fig. 4.3. Change of the energy carried away by the shock wave, at increase of the height of 
limiting cylindrical envelope 
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Fig. 4.4. Levels of the maximal blast pressure P, and the amplitude of blast wave AP at 
initial pressure Py = 10° Pa for various radii of envelope (Ryz = 1 cm): 1 - R = «0; 2-R= 
25 cm; 3 - R = 35 cm; 4- R = 735 cm; 5- R= 50cm; 6- R= 20cm; 7-R=150cm 


P kPa 
soo fmkP 


AP. 
sD soo 


Fig. 4.5. Levels of the maximal blast pressure P,, and the amplitude of blast wave AP at 

initial pressure Po = 10° Pa for various radii of envelope (Rug = 8 cm): 1 - R = c; 2-R= 
75 cm; 3- R= 150 cm, 4- R = 100 cm; 5- R = 35 cm; 6-R = 50 cm; 7-R = 25 cm; 8-R = 
20 cm 
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Fig. 4.6. Dependence of the factor of reduction of the explosive wave amplitude K on the 
relative distance R/Rur 


Fig. 4.7. The diagram for a choice of the effective thickness 6 of protective envelope 
ensuring the catching of fragments: I — the fragments penetrate the envelope, II — the 
fragments do not penetrate the envelope, III The zone of uncertainty [4.4] 


The scheme represented in the fig. 4.1, is rather effective for catching the fragments 
at the known thickness of cylindrical cowling. Using the information given in [4.4], it is 
possible to construct the diagram (fig. 4.7) for a choice of the effective thickness 6 of 
protective envelope ensuring the catching of the fragments. The diagram is constructed in 
Sadovsky — Hopkison coordinates. The abscissa is the reference distance R/G°**, and the 
ordinate is the specific thickness of the concrete envelope 6 /G°**. In the field of 
parameters I the fragments penetrate the envelope, in zone II — the fragments do not 
penetrate the envelope, and III is the zone of uncertainty. The diagram is correct if internal 


space of a protective envelope is filled by air under normal atmospheric conditions. 
4.3. | Vacuumization of gaseous medium around of a charge 


From the dependences of the blast wave overpressure AP and impulse of 


compression phase / follows, that AP ~ Py’**, and J = Py”. Thus, the trivial way of 


reduction of the blast loadings is the decrease of initial pressure Po in the gas layer 
separating HE charge from protected object. However, the vacuumization of cavity around 


HE charge requires the placing of the charge in special tight containers and actually 


excludes the application of the combined ways of protection. Nevertheless in special cases 


the given way of explosive loadings regulation can be considered as rather effective. 
4.4. Shock waves attenuation by perforated plates 


One of the effective ways of blast waves amplitude decreasing is the placing of thin 
perforated plates (PP), grids and other permeable barriers across the gas flow behind the 
wave. The shock wave propagation with the constant parameters behind the front through 
PP was studied in [4.5 ...4. 12]. However, in practice it is necessary more often to deal with 
SW with non-stationary pressure behind the front. The force influence of such waves 


depends on the overpressure at the front AP and the impulse of the compression phase: 
I= } AP(t)dt 
0 


Here +r is the duration of compression phase, AP(t) is the pressure — time 
dependence. It is impossible to make any practical conclusions by the results from [4.5... 4. 
12] concerning the change of parameters / and rt at SW propagation of final duration 


through PP. 


In [4.13] the results of experiments on the shock wave attenuation by the perforated 
screen are given. The shock wave was formed at destruction the gas-filled cylindrical 
envelope. It was shown, that the duration of SW compression phase does not change at the 


passage through the screen. 


In [4. 1, 4. 15] the attenuation of waves formed in environmental space by the grids 
and PP after HE detonation was investigated. In [4. 1] the data on measurements of the SW 


impulse behind a barrier are given. 


However at the conditions of the experiments in [4.1,4.15] the two-dimensional 
effects could essentially influence on the values AP an J behind PP, as the barrier were 
placed in the open space. Thus, for practically important one-dimensional processes in the 


pipes and the channels the use of results from these two references is wrongful. 


In [4. 18] the laws of interaction between the plane SW of final duration formed in 
a shock tube, with PP, positioned in one of the sections of the tube are found. The 
parameters of the incident SW (AP.,7.) and passed (transmitted) SW (AP.,7) were 
determined by the piezoelectric pressure gauges. The impulses of a compression phase of 
the incident SW (/.) and passed SW (/) were measured with the help of the electronic 
integrator of a signal from pressure transducer. Parameters of registration system provided 


the accuracy of AP measurement no more than 10 %, and /— 15%. 


The range AP = 0,05 ...0,25 MPa was investigated. Characteristic values of rt. = 
0,5... 2 ms. The permeability of partitions, i.e. the ratio of the area of apertures to the total 
area of PP, in experiments was @ = 0,009... 0,25. Thickness of partitions 3... 10 mm. The 
change of SW parameters at passage through PP can be described by coefficients 6 = 7. /r, 
, 9 = AP/AP, , and €=1/T1.. With the help of the gauge located directly inside the 
partition, the parameters of reflected SW AP« and rt» were determined. In a range under 
consideration was obtained, that the duration of SW compression phase at passage through 


PP does not change, i.e. + 7, = t and d= 1. 


In the fig. 4.8a the dependences of SW attenuating factor by pressure 77 on PP 
permeability are given. The curve | is plotted by the data from [4.5... 4.11] for the plane 
SW with the constant parameters behind front AP = 0,001... 0,2 MPa. The shaded area 2 
and points 3 are the data [4.1] and [4.15] accordingly for the spherical blast waves at AP; 
= 0,08... | MPa. The points 4 are the results of experiments for plane SW of a triangular 


pressure structure, which at 0,009 < @< 0,3 can be described by the empirical relation: 
n=1.140” (4.1) 


In the fig. 4.84, the dependences of factor € on PP permeability are shown. Area I 
corresponds to the data from [4.1], points 2 are the results of the experiments [4.18]. The 
curve 3 is constructed by (4.1) in coordinates é-@. Thus, within the experimental error for 
the plane waves of triangular pressure structure €=77. In the case of spherical blast waves 
[4.1] we have another pattern. If the pressure at the front considerably decreases behind the 


barrier, then the impulse of compression phase is close to initial one, and in the field of 


permeability a> 0,04 even exceeds it. It is probable, that in conditions [4.1] the increase of 
the impulse is connected to the outside flow around the barrier, as well as the reflected 
wave. It is shown in [4.12,4.13], that at a < 0,1... 0,2 the pressure in the wave reflected 
from PP only slightly differs from the normal reflection pressure. In the fig. 4.9 the results 
of y=AP-+/AP» measurements (AP» is the overpressure at reflection at impenetrable barrier) 


for the incident SW with a triangular structure of pressure are presented. 


Fig. 4.8. Influence of PP permeability on the parameters of shock waves: a — the factor of 
SW attenuation by pressure n; b — the factor of SW attenuation by impulse & 


It is seen, that at ~< 1, the pressure AP: differs from AP» less than in 30 %. With 
the increase of AP, this distinction decreases. The measurement of AP, by the pressure 


gauge located directly in PP, allows to estimate the effective SW influence on the partition. 


The detailed study of the process of the formation of the SW reflected from the wall, 


equipped with an aperture, is presented in [4. 17]. 


Fig. 4.9. Results of y7=AP «/AP 9 measurements. Points 1 correspond to AP. =0.2 MPa. 2 - 
0,08 MPa. Curves 3,4 are numerical calculations [4.7] for the plane SW with constant 
parameters behind the front accordingly at AP, =0.2 and 0,08 MPa. 


For the practical applications it is essential to study an opportunity of using several 
plates (so-called «package») for SW suppression. Consider a case of two consequently 
located PP with the permabilities a and a. If the distance between the partitions / 1s large 
enough in comparison with the distance from the last PP to the place of measurement, then 
it is possible to consider the plates «as independent» and effective factor of the SW 
attenuation by amplitude 7). will be determined by the values 7, 72 for each PP. Then 


taking into account (4.1) we have: 


m= m Mm 1.3 (a a). (4.2) 


In the case, when the measurements are carried out at sufficient distance from the 
«package», the repeated SW reflections in a gap between partitions will result into the 
generation of a stronger wave. For the effective permeability of two plates 77:2 in the case 
of «mutual influence» of PP the semi-empirical formula is offered in [4. 1] at 0,01 <q , 


a <0, 13): 
Oh. = Oh Oy (a +a)" (4.3) 


The «package» of two PP with the permeability 0,07 and 0,12 was investigated at / 
= 25 mm. The partitions were positioned in various combinations. The distance to the place 
of measurement from the last PP was 13 cm. In the fig. 4.10 the dependences of the 
attenuating factor 77,2 on the permeability of the second partition are given at the fixed a. 
Here, the curves 1, 2 are calculated by (4.2); 3,4 — experiment; 5, 6 are calculated by 
(4.1), (4.3); a = 0,07 for 1, 3, 5 and 0,12 for 2, 4, 6. It is seen, that the value 77;2 is higher, 
than in the case of «independent» partitions, but is not described by the formula (4.3). It is 
necessary to note, that the duration of the SW compression phase passed through the 


«package», is increased and, hence, the factor € in this case can differ from 77. 
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Fig. 4.10. Dependences of the attenuating factor 412 on the permeability of the second 
partition for the fixed a. 


Alongside with triangular compression waves the problem of interaction of short 
rarefaction perturbations with PP was investigated. In the experiments with the partitions 
of a= 0,14; 0,07; 0,012 the measured attenuation factors by pressure were 7 = 0,54; 0,34; 
0,17 accordingly. From the fig. 4.8a it follows, that 77 for the rarefaction waves are 1.5...2 
times larger than the corresponding values for the compression waves. The duration of the 
perturbation at the passage through PP does not change, and, hence, the impulse in the 
rarefaction wave decreases proportionally to the amplitude, as well as for the plane 


triangular compression waves. 


Thus, the perforated plates are the effective way of attenuating of the non-stationary 
pressure waves in the tubes and the channels not only by he amplitude, but also by the 
impulse. By the appropriate selection of a, and using of several partitions it is possible to 


vary the parameters of the transmitted SW and the reflected wave in a wide range. 


4.5. Shock waves attenuation by screens from granular materials 


In works [4.14,4.19,4.20,4.24,4.25] the propagation of SW in bulk and porous 
media is investigated. The obtained fast attenuation of waves allows assuming, that these 


media can serve as an effective method of SW suppression in channels. 


For the determination of attenuating effect generally it is necessary to consider a 
complete pattern of the wave propagation inside a porous sample. The experiments [4.19] 
testify that behind the plane SW, penetrating in volume with regularly loaded particles, 
there is an ordered wave structure, and the forward front has the poorly expressed breaks. 
The system of the reflected waves, which generates non-uniformity of thermodynamic 
parameters and speed across the section of the channel, is considerably weakened behind 
the front at the distances about several diameters of granular particles. The marked features 
allow assuming, that for the description of dynamics of SW attenuation in considered 


systems the one-dimensional (hydraulic) approach is applicable. 


The interaction of a flow with particles results in losses of the impulse of gas and 
redistribution of energy because of formation of the secondary jumps and the flow 
turbulisation. Factor of resistance of particles can be presented as the sum cy-cg+cy. The 
surface friction, as a rule, is small in comparison with the resistance of pressure (shape) 
(1.e. Cy” << cy ), and they can be neglected. To calculate cy in motionless layer of particles 
with a diameter d in conditions of the steady gaseous flow with density p the Ergan 


equation can be used: 
Cy * Cjp = 75/ Re + 0,875, (4.4) 


Where Ke = epdu / (1 — &)/u at the porosity 0,4 < e< 0,65. From (4.4) follows, 
that at Re> 10° the values cy tend to the constant value of 0,875, i.e. the medium obeys the 
square-law law of resistance. In the shock waves propagating in porous media of bulk 
density, such values of Re number are reached at the overpressures larger than 5 kPa (at d 


~ | mm). 


From the analysis of the equations of the perfect gas flowfield in a tube with 
friction losses follows, that the SW intensity at the exit from the screen with thickness / is 
defined by the incident SW Mach number and dimensionless parameter @ = 1.75 (1 — 


é)h/ed, describing the properties of the granular partition [4.22]. 


The results of the experiments on shock waves attenuation by the screens made 
from various granular materials (see tab. 4.1) with parameters close to that, used in 


practice, are given in [4.22]. 


Table 4.1 Screens, used in experiments [4.22] 


Size of 
Material Form of | the particles d, | ¢ Thickness of the screen /, 
particles mm mm 
Polythene Round shape | 3,9 0,38 10, 25, 60 
Steel Cylindrical O10 x 12 0,41 35, 110, 160 
35, 60, 110, 160 


60, 160, 250 


Calculated and measured ratios of AP/P, (AP is the overpressure in SW transmitted 
through the screen) versus the Mach number M, and the parameters of granular screens are 
presented in the fig. 4.11. From a fig. 4.11 it follows, that at O= 30 the amplitude of the 
passed SW practically does not depend on the initial intensity of waves in a range of My 
numbers under consideration. Hence, for the shock waves suppression the parameters of 


screen should satisfy to the condition O> 30. 


At h/d + | the parameter @, and, hence, and amplitude of the wave depend 
basically on the values Mo and «. It is necessary to note, that the ratio for pressure losses 
used in calculations, is valid under the condition //d >> 1. Nevertheless, similar result was 
observed in the studies on the shock wave suppression by perforated plates: the SW 


attenuation is defined only by permeability of the plates and the incident Mach number. 


Fig. 4.11. Change of amplitude of the shock wave: 1 — 3 — calculations [4.21], 4 — 6 — 
experimental data (granular screens), 7 — 9 — experimental data (perforated partitions), 
Mg = 147 (1,4, 7), 1,28 (2,5, 8), L135 (3,6.9) 


Points 7...9 in the fig.4.11 represent the experimental data for the perforated plates 
with permeability 0,44. Comparison with the experimental results for the screens made 
from granular material (¢ = 0,38... 0,50) shows, that in the limiting case at /d > 1 the bulk 
layer is simulated by the perforated partition with permeability equal to é (if to neglect the 
possible change of ¢ at h/d > 1 in comparison with h/d >> 1). In this connection, it worth 
to compare the action of SW on the perforated partition and the granular medium by other 
parameter, for example by the overpressure AP» in the wave reflected from the bulk 


medium. 
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Fig. 4.12. Change of amplitude of the reflected shock wave versus the porosity of a 
material: 1... 3- experimental data with the perforated partitions, points — experimental 
data with granular screens, 4... 6 — pressure of reflection from rigid wall, Mo = 1,47 
(1,4), 1,28 (2,5), 1,15 (3,60) 


Dependences of AP+/AP» on the parameter ¢ for the screens with minimal length 
(h/d = 2.2 ...3.2) are presented in the fig. 4.12. Also shown are the results for the SW, 
reflected from the perforated partition, and the ¢ value corresponds to the permeability of 
the partition. The satisfactory agreement between the screens and the perforated plates is 
observed. Slightly higher pressures of reflection for the screens is connected to the absence 


of the gaps even in the layers h/d = 2... 3. The pressure of reflection from the screen 


represents also the practical interest, as it is the value that defines a choice of strength 
properties of the screen. As well as in the case of the perforated partitions, the value AP: is 
lower than the pressure of reflection from a rigid impenetrable wall, decreases with the 
increase of ¢ and weakly grows with the increase of the screen length. Many of the 


conclusions are confirmed recently in [4.23]. 
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CHAPTER 5. APPLICATION OF POROUS SCREENS AND FILLERS FOR 
PROTECTION AGAINST BLAST LOADINGS 


5.1. Effect of the porous screen on blast loadings 


The interaction of plane shock waves with porous compressed materials made from 
polyurethane foam (PUF) and plastic foam was investigated in [5.1...5.5]. It was shown, 
that the maximal pressure amplitude at the wall under the PUF layer can considerably 
exceed the shock wave pressure at normal reflection from the rigid wall in the absence of a 
porous covering. For the theoretical description of the effects observed in [5.4, 5.5] the 
method based on representation of two-phase medium (PUF + gas) by the equivalent gas 
with changed in comparison with pure gas thermo-physical properties is offered. This 
method allows calculating the maximal loading and understanding the dynamics of the 
reflected wave. However, within the framework of equivalent gas model some specific 
features of the phenomenon revealed in [5.1... 5.4] are not described. The reasons of 
dependence of the pressure peak on the height (thickness) of a layer and the absence of the 
«shelf with constant parameters in pressure profile at the wall (at a level of peak pressure) 
remain unclear. Alternative model of the phenomenon based on consideration of the 


sample movement at sudden loading is described below. 


The essence of nonlinear process of the porous layer response to the blast loading 
can be presented as follows. At the moment of air shock wave reflection from the border 
gas - PUF a porous layer at the end plate of the shock tube is instantly loaded by the gas 
piston with the pressure close to the pressure of normal reflection. If the pressure 
considerably exceeds the characteristic value of compression stress of the porous material, 
then, at the initial stage, the force, counteracting to movement, appears negligibly small. 
During the development of deformation the resistance of a sample grows under the certain 


law. 


Let's present (fig. 5.1) the column of the compressed porous material of the unit 


cross-section area of mass m and height / by the equivalent mechanical system with one 


degree of freedom consisting of a point mass m and a combination of ideally plastic 
Coulomb element with a zero restoring force (region 1) and elastic elements with elasticity 


factor A and damping factor a (region 2). 
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Fig. 5.1. Calculated dependence of the relative loading coefficient 6 on time tat h = 50 
mm and P = 1.4 MPa. Dashed line corresponds to experiment [5.5] 


Let's consider the movement of mass m under the action of the suddenly enclosed 
constant loading P , corresponding to the overpressure in the reflected shock wave. In the 


region | the body gets the speed 
vo=(Pp')” 


where p is the density of the material (for PUF p = 33 kg/m’). The further motion 


of the mass is governed by the equation 


mx" + ax'+ kx =P 


(5.1) 


with the initial conditions 


x (0) = 0, x' (0) = vo. 
(5.2) 


"? is the critical damping 


At various relationships between a and do (ao = 2 (km) 
coefficient) three types of the solution of the equation (5.1) can be obtained. For the 
problem under consideration the case a<dag is the characteristic one. Pressure loading of a 
rigid support at an arbitrary time />0 is R=Ax. The relative loading coefficient 6 =R/po is 


found from the solution of the equations (5.1) and (5.2) to be: 


6 = exp(—atr) L(y —a) B' sin Bt—cos Bt] | 
(5.3) 


where a =a/2m;B=(k/m-a?) sy =hy/P 


The maximal value of the relative loading coefficient 6, is attained at the moment 


of time f,): 


wp eel e-ae-o) 


(5.4) 


Let's estimate the speed vo, with the assumption, that the mass ™m is initially located 
at the center of mass of PUF sample and the compression stress of a sample is much less 
than applied loading P. If to take into account, that at 40 % deformations the characteristic 
magnitude of compression stress of PUF is 10° Pa, then the last assumption is valid for 
P>>10° Pa. Then vo = (0.4 Pp”)’”, and the speed of the free border of the sample is 2vo 
and at P = 0.15...0.22 MPa is 85...103 m/s. The speed of the border, measured at the same 
conditions, is 95 m/s [5.5]. The calculated speed of the border PUF - gas does not depend 
on the height of the sample, that prove to be true by the measurements in [5.4, 5.5]. The 
magnitude of & for PUF is determined on the basis of one of the experiments (/ = 50 mm, 
P= 1.4 MPa) and from the condition of equivalence between the values of 6,, in the model 


and in the experiment. The / value is (4...4.2) 10° N/m (assuming @/do = 0.5...0.6). The 


dependence of loading 0=6/1) at h = 50 mm and P = 1.4 MPa is indicated by continuous 
line in the fig. 5.1. The dashed line corresponds to the experimental data [5.5]. The 
presence of a smooth curved part at the measured profile d/t) is connected, probably, to 
nonzero restoring force in the region | (fig. 5.1), and it is possible to explain the 


deceleration of the pressure decrease by the hysteresis phenomena at PUF loading. 
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Fig. 5.2. Dependence of the maximal relative loading coefficient 6, on the height of the 
PUF sample at different pressures P = 0.32 (1); 1.4 (2) and 5.0 MPa (3). k = 4.2 10° N/m’ 
and a=0.6 ao. The points 4 and 5 correspond to the experimental data [5.5] for P = 0.32 
and 1.4 MPa. 


The dependences of the maximal values of 6,, on the height of the PUF sample at 
different pressures P = 0.32 (1); 1.4 (2) and 5.0 MPa (3) are presented in the fig. 5.2. The 
specified dependences were calculated for the same magnitudes of k = 4.2 10’ N/m? and 
a=0.6 ao. The points 4 and 5 correspond to the experimental data [5.5] for P = 0.32 and 1.4 
MPa. In the pressure range under consideration the maximal relative loading coefficient 


Om can be approximated by the formula: 


6, =~ 140.3(khPt) 


As it is seen, the maximal loading increases with the decreasing of P (i.e. with 
reduction of the incident shock wave intensity). Note, that in the case of extremely small P, 


comparable to the PUF compression stress, the value of 6, expected to be decreasing. 


Thus, the character of loading changes at the shock-wave action on the structure 
with porous compressible covering. In the absence of the covering the structure is loaded 
quasi-statically. The presence of the covering results in the occurrence of the pulse 
component, which have to be taken into account in strength calculations. The suggested 
model correctly describes the phenomenon of amplification of shock-wave loading on a 
rigid support covered by the porous compressible material, reflects the dependence of the 
loading on geometrical sizes of the covering and on the pressure at the surface. On the 
basis of the proposed model the method of determining of shock waves loadings on 
structures faced by porous compressible materials can be developed. The solution of the 
inverse problem is possible, when the elastic properties of the covering are determining by 


known loading. 


5.2. Protection against explosive loadings by porous screens 


The expansion of technological processes involving explosion demands the further 
development of the reliable ways of explosion localization (damping, suppression). 
Widespread method of explosion localization is the use of explosive chambers (EC) [5.6, 
5.13]. An issue of the day in technology of material processing by explosion is the increase 
of maximal weight of HE charge to be detonated safely in the chamber. One of the 
solutions is the decrease of stress forces in the walls of the chamber by filling the cavity by 
various damping media and protection screens [5.7]. The reason of the decrease of stress 
forces, apparently, is that the damping material gives rise to an essential stretching of the 
pressure impulse in time and, subsequently, to the quasi-static loading of EC, at which the 


main factor determining the stress in the walls, is the lowered pressure, instead of impulse. 


At the same time, as shown experimentally in [5.8], the porous material placed at 


the rigid wall (as it was done in [5.2...5.4]) can increase the pressure amplitude of the 


reflected stationary wave. It causes the necessity of more careful research of the 
opportunity to use the given materials for loadings suppression and EC protection, that was 


not considered in [5.13]. 


In reference [5.8] the screens of various structures containing porous media are 
considered and the results of research of damping properties of porous media filling the 
internal volume of the chamber are offered. The screens were placed at the EC wall to 
decrease the impulse loading. The experiments were performed for the purposes of 
optimizing the parameters of porous fillers able to damp the loading. The screens made 
from foam-rubber (density 40 kg/m’) of different thickness, as well as alternating layers of 
steel and foam-rubber (stratified screens) were investigated. The loading was carried out 
by a spherical RDX charge (weight 0.15 kg) located at distance 0.5m from the wall of EC. 
The screens were mounted at EC wall, and the pressure gauge was set into the wall. The 


results of experiments are presented in table 5.1. 


Table 5.1 Parameters of loadings on walls depending on the barrier 


structure and thickness [5.8] 


Number — of | Thickness of | Maximal, Duration of the | Type of the barrier 
experiment impulse, Ls 


foam-rubber 


foam-rubber 


foam-rubber 


stratified screens 


stratified screens 
stratified screens 


Stratified screens were used in the experiments 116,117 and 124. Structure of the 


screens and the thickness of layers were following: 


e Experiment 116: steel - 2,5 mm + foam-rubber - 30 mm + steel - 2,5 mm + 
foam-rubber - 30 mm (2 - tier screen). 

e Experiment 117: steel - 5 mm + foam-rubber - 60 mm (1 -tier screen). 

e Experiment 124: 3 layers, each consisting of a steel plate (thickness 2.3 


mm) and foam-rubber layer (thickness 30 mm) (3-tier screen). 


As evident from the table, the use of foam-rubber screens with layer thickness < 80 
mm, as well as in [5.1,5.4], results in the increase of the SW amplitude at the wall behind 
the screen. Thereof the SW amplification at reflection could occur, as can be seen from 
experiments 111 and 105. Screens with layer thickness 100... 150 mm result in the 
decrease of the pressure amplitude (in comparison with the conditions of EC walls loading 
by explosion without the screen, experiment 102). Typical pressure records at walls are 


shown in the fig. 5.3. 


Fig. 5.3. Typical records of pressure at the walls of the explosive chamber [5.7,5.8] 


The comparison of the stratified screens of equal thickness and mass (differ only by 
the amount and the distribution of the layers, the experiments 116 and 117), shows, that 
two-layer | —tier screen is the most effective for the damping. The two-layer screen in view 
of the stretching of the impulse and the reduction of the pressure is better, than the screen 
of the larger thickness and mass consisting of three pairs of layers steel — foam-rubber (the 
experiment 124, 3-tier screen). It is connected with the fact, that the SW splitting in the 
screen with large number of layers (in acoustic approach leading to the attenuation of the 
head wave) is suppressed by the nonlinear effects arising at imposing of the reflected 
waves on the head wave. The impulse, received by the plate due to the interaction with 
SW, does not depend on structure of the screen and under constant conditions of loading 


remains constant. 


Kinetics energy F of the plate is larger for the lower mass © of the plate at the 
constant acting impulse 7. From this point of view it is favorable to have the massive 


elements of the screen on a surface perceiving the impulse. 


It is revealed experimentally, that at the constant parameters of the blast loading on 
the screen the most effective among the investigated structures are the two-layer screens 
consisting of the layer of steel (from the side, perceiving the SW impulse) and the porous 
material. The efficiency of such structure proves to be true by the studies of the blast effect 
on the animals protected by body armour (flack jacket) (a steel plate + amortization layer 
from a porous material) [5.12, 5.14]. Of interest is to determine the optimum porous 
material for the maximal pressure decrease and find the relationship between the degree of 
damping, amplitude and duration of the blast impulse acting on the screen. For this purpose 
the RDX charges of the spherical form and weight 0.034 and 0.215 kg were used. 
Changing the weight of the charge and the distance X to the screen, the amplitude AP and 


blast impulse duration 7 was determined in each series of experiments. 


Table 5.2 Parameters of the loads at EC walls at the loading by SW generated by the 
detonation of the charge of weight 0.215 kg through the screen: steel - 2 mm, porous 
material - 40 mm by [5.8] 


Structure 
of the 
screen 


Without 
screen 


Pearlite 


Rubber 
crumb 


Table 5.3 Parameters of the loads at EC walls at loading by SW generated by the 


detonation of the charge of weight 0.034 kg through the screen: steel - 1 mm, porous 


material - 22 mm by [5.8] 
Distance to the screen R and (2) 


Structure 
of the 
screen 


Without 
screen 


0,76(35) 


Sand 
Foam- 
rubber 


Pearlite 


Rubber 
crumb 


Steel shot 


The results of experiments [5.7, 5.8] with two-layer screens of identical thickness 


and the structure differed only by the type of a porous material are presented in tables 5.2 


and 5.3. The amplitudes and the characteristic duration times of the loadings on EC walls 


behind the screens are listed. The characteristics of the blast impulse at the EC wall 


without the screen are presented for the comparison. The relative distances to the screen 


R=R/R,, where Ro is the radius of HE charge are given in the tables. 


As seen from the tables 5.1, 5.2 and 5.3 the preferable porous materials for damping 
screens are foam - rubber, rubber crumb and pearlite, which allow essentially reducing the 
amplitude of the pressure at EC walls with simultaneous increase of the impulse duration 
in a wide range of governing parameters. The screens containing sand and steel shot, are 
less effective for the damping. In a number of cases such screens result in the increase of 
the pressure amplitude at EC walls. It is important to note the dependence of damping 
properties of screens containing different porous materials, on the duration and the 
amplitude of the blast pulse. So, for the blast pulse generated by the charge of weight 0.215 
kg at the distance 1.4 m(R = 35) the degrees of pressure amplitude reduction differed a 
little for various screens. On the contrary, for R = 13 the essential difference in damping 
properties was observed. The research of the loadings at EC walls filled by porous 


materials is of interest for the localization of explosion. 


Fig. 5.4. Scheme of determination of blast loadings: 1 - HE charge; 2 and 4 - pressure 
gauges; 3- deformation gauges; 5 - wood chips; 6 - air [5.8] 


The reaction of a thin-walled spherical envelope with a diameter 0.45 m and 
thickness of walls 0.007 m, loaded by explosion of the concentrated HE charge 1, located 
in the center of symmetry of the envelope was investigated in the experiments [5.8] (fig. 


5.4). The charge was initiated by the detonator. The pressure (gauges 2,4) at the internal 


wall of the chamber together with deformations (gauges 3) was measured. The experiments 
were carried out for the cases, when the envelope was filled by the wood chips 5 with 


initial density 100 kg/m’ or air 6 at normal conditions. 


In the experiments with sphere filled by air at weight of charges 0,03 and 0,06 kg 
the peak stresses at walls were accordingly 240 and 370 MPa. The maximal pressure was 
17,6 and 15 MPa, and time of loading at the wall was 46 and 50 us accordingly. When the 
envelope was filled by the wood chips at explosion of charges of weight 0.03 and 0.06 kg, 
the parameters of loading and deformation are equal accordingly: pressure — 1,3 MPa, 


time of action 1800...1230 us, maximal stress 45 and 150 MPa. 


Clearly seen, that for the charge of weight 0.05 kg the filling of EC volume by the 
wood chips results in decrease of the maximal stress more than in 5 times, and for the 


charge of weight 0.06 kg of the deformation forces decreased more than in 2 times. 


The important result obtained in [5.7... 5.9] is the dependence of the maximal 
pressure at the wall of the explosive chamber on the degree of compression of porous 
material in SW. The pressure records at the wall of the explosive chamber filled by 
wooden chips are presented in the fig. 5.5. Weight of the charge in the experiments is the 


same. 


Introducing the designation y=o/p, where is the density of the wood chips 
before explosion; is its density after explosion, it is possible to obtain, that the pressure 
decreases with increasing compression degree of the chips (or with decreasing 7). So, the 


curve | corresponds to y= 0,9, the curve 2 — y = 0,86, and the curve 3 - y= 0,76. 


Fig. 5.5. Influence of preliminary compaction of porous media on blast loading [5.8] 


Thus as was established in [5.8], the most effective among the investigated types of 
stratified screens are two-layer ones consisting from a porous and dense (for example, 


steel) material, located at the side perceiving the impulse. 


Among the investigated porous materials in the two-layer screens the most effective 
are foam-rubber, rubber crumb and pearlite, allowing to reduce the pressure at EC walls 


with simultaneous increase of loading duration. 


The damping effect essentially depends on the amplitude and duration of loading, 


acting on the screen. 


Experimentally shown, that the maximal pressure at the EC wall depends on the 


degree of compression of porous material. 


The filling of EC volume by porous material results in the reduction of the maximal 


stress in its walls, at least for the investigated relative distances. 


5.3. Optimization of protective screens. Problem formulation 


The processes of SW propagation and reflection in porous and stratified screens at 
the loading by the pressure pulse are investigated in [5.9], where the problem of 
optimization was solves. The stratified screens made from steel plates and polyurethane 


foam (foam-rubber) were studied. 


The optimization of protective screens has the practical importance for designing of 
the explosion protection equipment with the purpose of reduction of its dimensions, weight 


and cost. 


The behavior of foam rubber is modeled in [5.9] by equilibrium mixture of ideal 
caloric gas (the first phase) and incompressible particles (the second phase). The governing 


equations in Lagrangian form in case of one-dimensional flow at plane symmetry are: 


vot Orv at or 
Be __povdr 
Ot ot dt 


L=1,=T;R=B =P3u, =u, =u; 
m +m, =1,p,=m,p,5e, = ¢,1;; 


FL = PART; Py = const; p, = ap; 


e=¢cT C=C. 4+@6,:0, +0; = 1 


where v=1/o 1s the specific volume of the mixture; # is the mass velocity; P is the 
pressure in the mixture; e is the specific internal energy of the mixture; c; are the specific 
heat capacities of the phases; »;, are the true densities of the phases; 9; are the average 
densities of the phases; m; are the volume concentrations of the phases; a are the mass 


concentrations of the phases. 


In case of the mixture with a small volume concentration of condensed phase the 


model take a form: 


Lav _ au Lau__aP de__ pov. 
vot Orv at Or’ Ot Or 
P= pRT;R=a,h,e=cT;3c=a,¢,+ a,c, 


Such a mixture can be considered as ideally caloric gas (R=const, c=const) with the 


specific heats capacity ratio y=(c+R)c. 


The loading of a flat layer A;A2 of rubber-foam with thickness d by a triangular 
pressure pulse P(/) was considered in [5.9] (fig. 5.6). The border Az is an absolute rigid 
wall. The mathematical formulation of the problem is following: it is required to find 
functions uv, P, p, v, e, T, Py, Px» Pin, P22 M1, M2 in the area z = {A)<rs<Az, O<f<tp}, 
satisfying to above presented set of the equations, with corresponding initial and boundary 


conditions. 


Fig. 5.6. The scheme of the screen loading by the blast wave [5.9] 


The initial conditions: at the moment / = 0 the layer A; Az is not loaded, therefore 


at (= 0: 0 = 0, P=0, pri=pri" , P2=P22" , P=Po , (A1SrSAd). 


Boundary conditions: at A; border, at ‘> 0 the dependence P(t) is given, at the 


absolute rigid wall A> the condition z(7) = 0 is given. 


P(t) has a form: 


P(i-t/t,), forO<t<t 
P(t)= o(l-t/tp), for D (5.5) 
0 , for? >f, 


Density of the foam-rubber at the initial moment f = 0: pp = 40 kg/m’, true 
(absolute) density of the condensed phase >>” = 1.210° kg/m’, density of the gaseous 


phase (air) 97)" = 1,29 kg/m’. 
5.4. Shock waves propagation in porous screens 


The parameters of the loading pulse are obtained at Py = 5.8 MPa, fp = 200 us in 


expression (5.5). 


The calculated and experimental maximal pressures at the absolute rigid wall for 


different thickness of the sample are presented in the fig.5.7. 


The analysis of the calculation results shows, that the model 2, in which the 
volumetric concentration of a gas phase (m);<</) was not taken into account, describes the 
experimental data unsatisfactorily, because during the interstice reduction the volumetric 


concentration of the condensed phase can become about unity. 
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Fig. 5.7. The maximal loading at the wall for screens of different thickness at constant 
parameters of the blast wave [5.9] 

The comparison of the numerical results obtained in [5.9] with the use of the model 
1, taking into account the volumetric phase, with experimental data shows their good 
correspondence, at least up to the thickness of the samples 70... 80 mm. The little bit 
smaller magnitudes of the experimental data at the thickness of the samples more than 100 
mm are explained, apparently, by multi-dimensional character of a physical task, and non- 
equilibrium of the phases in real conditions as well. Results on the basic model are 


discussed further in view of the volumetric concentration of the condensed phase. 


The calculated and experimental dependences of pressure at the rigid wall on time 


are given in the fig. 5.8 for the thickness of the sample 60 mm. 
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Fig. 5.8. The diagrams of change of pressure at the rigid wall [5.9] 


As seen from the numerical and experimental data presented in the fig. 5.7, for the 
triangular loading pulse with parameters Po = 5.8 MPa, fp = 200 us there is a critical 
thickness of a porous layer d*, below which the amplification (but for the attenuation) of 
the pressure amplitude in the reflected wave occurs. At the same time, as shown in [5.1], 
the amplification of stationary SW amplitude, i.e. at fp >, takes place for all thickness of 


the layer. 


In this connection the dependence of critical thickness of a layer on the parameters 
of the loading pulse is of interest. In table 5.4 the critical thickness for four-types of 


loading are presented. 


Table 5.4 Critical thickness of protective layers [5.9] 


No. Parameters of the pulse d*, mm 
Pressure Po, MPa tp, US 

l 1,0 130 18 

2 6,6 130 87 

3 10 130 112 

4 6,6 1300 850 
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Fig. 5.9. Critical thickness of the porous screen for constant blast 


loading 


It is seen, that the critical thickness of rubber-foam layer depends both on the 


duration of the pulse, and on its amplitude. 


The characteristic dependence (Po = 10. MPa, tp = 130 us) of the pressure 
amplitude at the absolute rigid wall on the thickness of the porous layer is shown in the fig. 


5,9, 
5.5. Efficiency of blast loading suppression by porous screens 


The loadings by a triangular pressure pulse of the stratified screen consisting of 
one, two or three identical tiers were considered in [5.9]. Each tier contains a steel plate 
and rubber-foam layer (fig. 5.10). The stratified screen is contiguous with the rigid wall. 


The metal plates were considered incompressible. 


Fig. 5.10. The combined protective screen: a layer of steel (1) and a layer of rubber-foam 
(2) 

At t= 0 the pressure P(/) is given at the border of the stratified screen from the side 
perceiving the loading, and 2?) = 0 at the absolute rigid wall. At the contact borders of 


various materials the condition of speed and pressure continuity is given. 


The calculated (curve 1) and experimental (curve 2) time dependences of the 
pressure at the rigid wall behind the stratified screen consisting of one tier are shown in the 
fig. 5.1la by the data from [5.9]. Thickness of the steel plate is 4 mm, rubber-foam layer 


— 60 mm. 


Fig. 5.11. Pressure at the rigid wall behind the stratified screen: a - steel (4 mm) — foam 
rubber (60 mm); b — 2 tiers: steel (2,5 mm) — rubber-foam (30 mm); c — 3 tiers: steel (2,5 
mm) — foam rubber (30 mm). 1 — numerical calculations; 2 — experiment 


The time dependences of the pressure at the absolute rigid wall obtained in the 
numerical calculations (the curve 1) and in the experiments (the curves 2) are compared in 
the fig 5.11 by [5.9] for the stratified screens, consisting from two tiers, where each tier 
contains a steel plate by thickness of 2,5 mm and foam rubber layer by thickness of 30 
mm. The similar comparison of numerical and experimental dependence of pressure at the 
rigid wall is given in the fig. 5.11 for triple-tier barrier, where each tier consists of a steel 


plate by thickness of 2,5 mm and foam rubber layer of 30 mm. 


From the comparison of dependences of the pressures for three investigated types 
of barrier (single, double and triple) it follows, that the most effective for the shock wave 
loadings damping is the single-tire barrier. This barrier has the least dimensions, weight 
and is most simple for the practical use. The characteristic dependences P(/) at the rigid 
wall for three types of single-tire screens are presented in the fig. 5.12. The thickness of the 
steel plate was 4 mm in all experiments, and the thickness of the foam rubber was 40 mm 


(the curve 1), and 60 mm (the curve 2) accordingly. 
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Fig. 5.12. Change of P(t) at the rigid wall for single-tier screens 


All types of the investigated screens effectively transform the shockwave impulse 
of pressure into direction of the reducing of pressure amplitude with the simultaneous 


increasing of duration. 


5.6. Optimization of stratified porous screens 


The problem of optimization of stratified screens is considered in [5.9]. The single 
— tier double-layer screens were investigated as the most effective (system a plate + foam 


rubber layers ). 


Dependences of the maximal pressure at the absolute rigid wall on the thickness of 
foam rubber layer are shown at the fixed thickness of steel plates (fig. 5.13) and on 
thickness of the steel plate at the fixed thickness of foam rubber layers (fig. 5.14). The 


calculations were carried out at parameters of pulse loading: Po = 5.8 MPa, fp = 200 us. 
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Fig.5.13 Change of the maximal pressure at the rigid wall at various thickness of the foam 
rubber layer d, 
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Fig.5.14 Change of the maximal pressure at the rigid wall at various thickness of the metal 
Plate dn 


When there is no restriction on the thickness of the stratified screen, the decrease of 
the maximal pressure at a rigid wall can be achieved by the simple increase of the 


thickness of layers. Thus the more thickness of layers, the less maximal pressure at a wall. 


Let the permission pressure P at the protected wall is given. As shown in [5.4], the 
given pressure can be achieved at different thickness of layers, increasing the thickness of 
that or other layer. The choice of the optimum screen is reduced to selection from all 
stratified systems by the criteria of minimal weight or cost. This mathematical task in [5.9] 
is solved as the search of the thickness of layers d,, (metal) and d, (foam rubber), which 


minimize the function 
Pin Am + Py d, > min 


restricting the maximal pressure at the wall Ping < P*. Here p, and , are the 


density of steel and foam rubber accordingly. 


The analysis of dependences shown in the fig.5.13 and 5.14, has allowed making a 
conclusion that the value P is the function of the product d,,d,. with sufficient degree of 


accuracy. This conclusion proves to be true by the data from table 5.5. 


From the condition Pia=P* we obtain P (dnd) = P*, and therefore d,,d,, =const. 


Denote d,,d,. =a. Then d,=a/d,, and the required function will take the form: 
Pin Am + Py A/a > min 


the minimum of the given function was obtained in [5.9] at 


dn=(Or aipy” d,=(Pm wp)? 
Thus, knowing P*, we can find a from the diagrams (fig.5.14 and 5.13), and finally 


— the necessary thickness of layers. 


Of interest is the maximal degree of decreasing the pressure at the rigid wall at the 
given thickness of the single-tier stratified system d/,, + d, =const. Besides, it is important 
to elucidate the influence of characteristic parameters of triangular loading Po and fp on the 
damping properties of the screens. The ratio of layers thickness d,,, /d,, has to be determined, 
at which the minimum of pressure at the absolute rigid wall is reached, under condition 


that thickness of stratified system is constant, i.e. d,, + d, =const. 


Table 5.5 Maximal pressures at the wall for different thickness of layers of steel and foam 
rubber in the screen 


No Prnax, 10° Pa Thickness of layers, mm And, 
din | d,. 


Fig. 5.15, Attenuation of the long-wave pressure pulse: Po = 1 MPa, tp = 1000 pus, dm + di, 
=50 mm. 
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Fig.5.16 Attenuation of the short-wave pressure pulse: Po = 10 MPa, tp = 100 us, dy + d, 
=50 mm. 


Dependences of the maximal pressure at the rigid wall P,,,. on the ratios d,, /d, are 
shown below, when d,, + d, =50 mm, and the loadings were characterized by the 
parameters: Po = 1 MPa, fp = 1000 us (fig. 5.15) and : Po = 10 MPa, tp = 100 us (fig. 5.16). 
The parameters Po and fp in both cases are chosen from the reasons of preservation of the 
constant impulse Pofp /2, acting on the screen. It is seen, that the effective decrease Pinay 1S 
observed up to d,, /d,, =0.5, and the minimum is reached at d,,, /d, ranging from 0.3 to 0.9. 


Further, there is an insignificant growth of loading with the at increase of dy /d, . 


The position of the minimum of pressure weakly depends on the characteristics of 


the applied pulse. 


The results of similar calculations for the case d,, + d, = 80 mm and parameters of 


the loading pulse Po = 10 MPa, ¢p = 100 us are presented in the fig.5.17. 


The characteristic pressure profiles at the rigid wall for the two-layer screen with 
the thickness of 80 mm and different ratios of thickness of layers of steel and rubber-foam 


are shown in the fig. 5.18. Such screen results in strong transformation of the initial pulse. 


The amplitude of pressure decreases in 100 times (curves 4), and the duration is increased 


in 300 times in comparison with parameters of the initial pulse. 
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Fig. 5.17. Dependence of the maximal pressure at the rigid wall on the ratio d,/d, : pulse 
Po = 10 MPa, tp = 100 pus; dy + d, = 80 mm. 
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Fig. 5.18. Characteristic pressure profiles at the rigid wall for the two-layer screen with 
the thickness of 80 mm: 1 —dy, = 10mm, d,. = 70mm, 2- dm = 20mm, d, = 60mm; 3- din = 
30mm, d, = 50mm; 4- dy = 40mm, d, = 40mm by [35.9] 


Let's describe the effect of volumetric concentration of the condensed phase (foam 
rubber) on the maximal pressure at the rigid wall. Consider the two-layer screen consisting 
of the single tier. The tier contains a steel plate with thickness d,, and foam rubber layer 
with thickness d,. The pressure P(/) is given by the dependence (5.5) with the parameters 
Po = 5.8 MPa, fp = 200 us. Average density of foam rubber varied from 0.04 up to 0.12 


g/em’. 


The maximal pressures at the rigid wall for stratified screen with the thickness of 
foam rubber layer 80 mm, and the thickness of a steel plate from 2 to 8 mm are shown in 


the fig.5.19. 


The analysis of the dependences obtained showed, that the increase of volumetric 
concentration of the condensed phase results in the increase of the maximal pressure at the 
rigid wall. However this effect is insignificant and shown to be revealed at the small 


thickness of the steel plate. 
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Fig. 5.19. Dependence of the maximal pressure at the rigid wall on the density of a porous 
layer in the screen, consisting from foam rubber with the thickness 80 mm and the steel 
plate with the thickness from 2 to 8 mm 


The dependences of the maximal pressure at the rigid wall on the thickness of foam 
rubber layer are presented in the fig.5.20 for different volumetric concentrations of the 
condensed phase in foam rubber. The effect of volumetric concentration of the condensed 
phase on the parameters of the shock wave propagation is insignificant for the system 


under investigation. 
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Fig. 5.20. Dependences of the maximal pressure at the rigid wall on the thickness of a 
foam rubber layer for different volumetric concentrations of the condensed phase in foam 
rubber. 


5.7. Recommendations for the selection of the stratified screen 
Thus, the observations and calculations in [5.7... 5.9] show, that: 


if For the loadings with triangular profile, there exists the critical thickness of 
a porous layer d*, below which the amplitude in reflected SW from the rigid wall in a 
porous layer increases in comparison with the amplitude of applied loading Po. At d>d* 
the pressure amplitude at the wall behind the screen will be less, than in initial loading. 


This paradox is confirmed in [5.11]. 


The critical thickness ¢d* depends on the amplitude Py and duration fp of the blast 


loading. 


2 The most effective among the stratified screens is two-layer screen. The 
denser layer should be located from the side perceiving loading. Moreover, the two-layer 
screens are easier, have smaller weight and dimensions for the given degree of pressure 
amplitude attenuation. The degree of the pressure amplitude attenuation is proportional to 


the layer thickness. 


The maximal pressure attenuation is reached at the ratio of thickness of the dense 
layer to the thickness of porous layer — 0.8 and practically does not depend on the 


amplitude P, and duration fp of the triangular pulse. 


3. The increase of the porous layer density from 40 to 120 kg/m’ practically has no 


effect on the damping degree both for two-layer screens, and for porous layers. 


5.8. Protection of the human body against explosive loadings by special 


equipment with use of porous materials 


The scientific and technical level of special protective equipment (flack jackets, 
helmets) permits to decrease the effect of loadings in the area of thorax and stomach in 


40... 50 % and in the area of head in 90 % [5.12... 5.15]. 


The reduction of the damage level for bioobject due to protective clothes is 


illustrated in the fig. 5.21. 


Nomogramm in the fig. 5.21 allows comparing the damage levels the unprotected 


body (a) and the man protected (4) by special equipment (miner) [5.15]. 


The direct studies of the flack jackets influence on the heaviness of the blast trauma 
at animals (pigs) have shown the protective effect of the «rigid» structures (steel armour 
plate and damping support from gas-filled, cellular polymeric material) [5.12,5.14]. As a 
sample, the flack-jacket of the second class of protection “Module 3M» was chosen. The 
protective elements in this jacket are executed as figured orthopedic plates. As a source of 
blast wave the HE charges of two types were used: TNT and LDE of complex structure. 
The experiments were carried out in open field (HE weight 5 kg) and in enclosure (1 kg 
HE). In each experiment two comparable pairs of animals in flack jackets and without 
protection were settled down equally spaced from the center of explosion, left side to 
explosion. 15 comparable pairs of animals were obtained under different experimental 
conditions (HE type, detonation conditions, distance from explosion). The results of the 
blast effect were determined as a degree of the heaviness of the blast trauma on the basis of 
the morphological analysis of animals at macro- and microlevel. Besides, the biochemical 
studies of the blood were performed. The analysis of the data obtained shows, that the total 


degree of the heaviness of the blast trauma for the animals in jackets is 25 % lower, than 


for the unprotected animals. This parameter practically did not depend on the conditions of 


experiment and characteristics of HE charge. 
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Fig. 5.21. Nomogramm for definition of the damage level of bioobject: unprotected (a) and 
protected (b) by special equipment, depending on distance and weight of HE charge. 
Horizontal lines: 1 — lower threshold of eardrums damage; 2 — 50 % of eardrums 
damage; 3 — lower threshold of lungs damage; 4-50 % of lungs damage; 5-99 % of 
survival; 6-90 % of survival; 7 — 50 % of survival; 8 — 10 % of survival; 9-1 % of 
survival 


The soft protective clothes are inefficient for counteraction to the blast loading. 
Moreover, there are direct indications to the amplification of the blast loadings by tight 
compressed materials that is harmful for the state of gas-filled organs of the human body 


[5.12,5.14]. 
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CHAPTER 6. EFFECTS OF SHOCK WAVES SUPPRESSION BY TWO-PHASE 
MEDIA 


When choosing the methods of influence on the parameters of blast waves, it is 
pertinently to use as the basic ones the information on the features of the waves 
propagation after the detonation of the concentrated HE charges. First of all, let's make 


some assumptions and remarks. 


As the energy transmitted from the gas phase into the liquid phase, practically does 
not cause the pressure increase and, thus, cannot maintain the shock wave propagation, it 1s 
necessary to expect essential increase of the effect of the absorption of the shock wave 


energy by two-phase medium when explosive source is placed into it. 


The experiments show, that at the detonation of HE charges placed directly in two- 
phase medium or at the formation of short perturbations in such medium by other sources, 
the intensity of wave disturbances is considerably reduced in comparison with the case of 
the detonation of HE charge placed in continuous medium. Depending on the volumetric 
contents of bubbles in the liquid the blast waves attenuation can be in 10... 50 times 
effective in comparison with a pure liquid at volumetric concentration of a gas within the 
limits of 1 % < a < 50 %. The difference of the considered scheme for the blast waves 
parameters decreasing from the scheme used earlier is, that the source of a wave is located 
directly in the two-phase mixture and the products of explosion or compressed gas at 
expansion contact to the two-phase medium without intermediate layers of continuous 


medium. 


6.1. Relationship of the blast waves parameters and initial conditions of an 


environment. 


The representation of the dependences of the basic parameters of spherical blast 


waves on the basis of the theory of similarity with the application of dimensionless values 


Ic > R 
nn Aaa TS —__; k= ——_ 
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permits to establish the influence of initial conditions of the medium (pressure Po, 
temperature 7) and sound speed co) at their arbitrary change. Let the parameters of 
explosion in medium with Po, 7o, co are known. It is necessary to understand the 
transformation of the intensity, the duration and the impulse of the wave for explosion in 
the medium with P); 4 Po, 71; # To, or ci: # Co. The sound speed in the atmosphere is 
changed due to the fluctuations of the temperature and the humidity of the air. The 


relationship reflecting these effects (m/s) is known: 


0 


5 ie 
¢,, =340(T, /273)” frost 
(6.1) 

Here 7, K is the temperature of the air, P,, is the water vapor pressure. 


The initial pressure of an environment varies depending on the height above the sea 
level. Let's consider the basic consequences from the theory of similarity. The adjustment 
for the SW overpressure will be a multiplier P,\/Po, ie. P21/P, ~ (Pi1/Po)??. At the 
assumption, that the pressure in the explosion products is much greater than the pressure in 
the environment, the explosion at higher initial pressure will produce the blast wave with 
the parameters increased proportionally to the initial pressure. The value of the dynamic 


head behind the wave also increases in the same proportion. 


More complicated is the adjustment multiplier for the pressure impulse in the 


passing or reflected wave: 


0,66 
fu [4h [=] 
rf Fy Ci 
It is possible to see, that the impulse increases with pressure, but decreases in 


inverse proportion to the sound speed in environment. 


Let's proceed to the relationships for adjustment multipliers for values describing 
the time of wave arrival, duration of compression and rarefaction phases. It appears, that in 
the environment with the initial speed of the sound cj;, and the pressure Pj,, all time 


intervals change in a proportion: 


Together with the pressure rise, not only the time intervals of the shock-wave 


process are shortened, but also the spatial scales of system change, where the impulse 


moving are essential. It is caused by the change of the characteristic length of explosion Ro, 


which in spherical case is 


Note, that the theory of similarity does not give the information on the dependence 
of the overpressure at the blast wave front on the relationship of the sound speeds in the 
explosion products c4 and in the environment exposed to the shock compression co and on 
the pressure in explosion products P4. Actually such relationship exists and in the first 
approximation can be represented by the basic expression from the theory of shock tubes 


for the pressure P» at the front of the blast wave: 
Pye (Bat) 
fy (Ly Mo 
(6.2) 


Here: 7% is the specific heats ratio in explosion products; % is the specific heats ratio 


in the surroundings. 
At P4 >> Po the following approximation is valid 
2 
(7s =) eer: (Pee = 1) | = Dy. | 27. +(%+ eee - 1) | 
(6.3) 


It is seen, that any event promoting the increase of the sound speed in the explosion 
products or the reduction of the sound speed in the environment, should be completed by 


some amplification of the blast explosive wave. 
6.2. Change of blast loadings parameters by two-phase medium 


On the basis of the above-mentioned information on the blast waves parameters in 


air together with the hypothesis about an opportunity of the two-phase systems 


representation as an equivalent gas with corrected thermo-physical properties, we shall 
describe some important ways of regulation of the blast waves parameters. We present the 
expected ways of the blast loadings regulation accompanying the carrying of the explosion 


source from the continuous to the two-phase medium [6.1,6.3,6.4,6.48... 6.50]. 


Basing on the concept of equivalent gas, we shall write the known expressions 


[6.36,6.37] for the speed of the blast wave from the strong spherical explosion as 
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Here w* is the energetic parameter [6.38, 6.39], o is the gas density, o is the 


concentration of particles of continuous medium, / is the time. 


For the constant energy of explosion in the same moment of time the speed of the 


blast wave in two-phase medium D* and in a gas D are related as: 


o.| ee)" | 140 ‘ 
D |a'p.| | (1+ 7)(1+7)(1+ 776) 


(6.5) 


Here o is the ratio of the specific heats for the continuous medium and the gas, 
n=0/p is the loading parameter (the ratio of the mass of continuous medium to the mass of 
a gas), 7=(1+70)(1+y770). With the velocity of the BW the distance passed by the wave 
also changes. If the distance passed by the blast wave in a gas is Rs; then the distance 


passed by the blast wave in two -phase mixture Rs* will be defined as: 


0,2 
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(6.6) 


The velocity of the blast wave at the distance RX for the strong explosion is 


D = oS = RY 
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(6.7) 


At the same distance from the center of explosion the blast wave velocity in two- 


phase mixture D* is lower than it's velocity in a gas D: 


0,5 
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(6.8) 


If the temperature of the particles is constant then 


D*/D = (1+ny°° 
(6.9) 


It is seen, that due to the heating of the particles (at yO >1) the effective wave 
attenuation is realized. The inertia of the medium is characterized by the value 77 and the 
ability to extract the heat from a gas - by the value yo. For the large kinematics inertia 


(77>>1) and high ability to the fast heating (v7d>>1) we have 


The reduction of the blast wave speed in water-air foams marked in the experiments 
(fig.6.1) have appeared closer to the case of isothermal particles, rather than to an 
equilibrium case. When the wave generated by the strong explosion propagates in two- 


phase medium, the mass of the substance involved into the movement increases. In the 


same moment of time after the explosion the mass of two-phase medium m* exceeds the 


mass of the gas m: 


14. 


m 


=(1+n)" [eerie ney 


(6.10) 


Fig. 6.1 Dependence of the SW speed on the relative (by the weight of the charge) distance 
for air (1) and foam (2) by [6.48... 6.50] 


The difference in the masses and the propagation speeds of the substance at the 
explosion in two-phase medium and in a gas leads to the change of the explosion pressure 
at the same distance. For the case of complete equilibrium the explosion pressure in two- 


phase medium will be P;*, and in a gas P) that gives: 


P1*/P, x+nyo)' (6. 
11) 


i.e. at R/Ro = 14 it is possible to note the decrease of overpressure in a wave, as 


shown in the fig. 6.2. 


The ratio of impulses of pressure in blast waves in foam /* and in a gas / is 


t.| (14+T)(1+7) } 
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(6.12) 


In two-phase medium with the frozen temperature of the particles, for example at 6 
= 0, the impulse in the wave exceeds the impulse in the wave propagating in the gas. At 
strong explosion in two-phase medium there are the zones, where /*//<1 and /*/7 > 1. The 


first situation is more probable far from the charge, and second one - nearby the charge. 


Not only in the strong explosion case the placing of the pressure waves source into 
two-phase medium will result in the decreasing of demolition action. The same effects will 


take place in the case of placing the high-pressure vessel in a dusty system. 


The filling of the high-pressure vessel by two-phase mixture will result in reduction 


of the part of the energy transferred to the shock wave at explosion products expansion. 
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Fig. 6.2. Change of pressure with distance for foam (1, 2) and air (3) by [6.48... 6.50] 


The effective suppression of shock waves also can be reached at detonation of HE 
charges in the tubes filled by the water-air foam or by the liquid with suspended gas 
bubbles. As shown in [6.1... 6.4,6.50], the physical model for the description of the blast 


waves attenuation process in the foamy structured medium can be offered. 


According to the model [6.1...6.4, 6.50] there is an opportunity of the approached 
representation of two-phase medium as a gas with an effective parameter of the adiabatic 
coefficient 7: Within the framework of such model the change of the shock waves 
parameters in two-phase medium can be investigated at their propagation into the direction 
of the acoustic resistance decrease. The cases of transition of the pressure wave from a 
continuous liquid in gas-liquid medium, from gas-liquid bubbly or foamy structure into the 
air [6.6, 6.16, 6.17, 6.22, 6.26] are considered. It was shown, that at these transitions the 
intensity of the shock wave decreases in 4... 6 times at the volumetric contents of the gas in 


the liquid about 5 %. 


The reduction of explosive loadings practically does not depend on the amplitude of 
the incident wave. For the gas-liquid systems of bubbly structure it was found, that the 
attenuation of the shock waves, transmitted from a liquid to a gas-liquid medium, is 
effective at the contents of a gas in a liquid (by volume) @ = 2... 4 %. The processes of the 
waves attenuation at crossing interphase boundary into direction of the acoustic resistance 


decrease can be described on a basis of a simple hydrodynamic model [6.51]. 


The practical importance of studying the change of the pressure waves parameters 
at transition through the interphase boundaries is, that it is possible to create a simple 
calculation model of the interaction of the pressure waves with bubbly screens in a liquid. 
Such model is constructed for the cases of normal and oblique incidence on a filter. It was 
shown that bubbly screens become ineffective for the attenuating of extended waves (for 
which the medium trends to be continuous by compressibility). The efficiency of the 
screens is reduced with the increase of the initial pressure in the medium. The larger 
efficiency of the shock waves attenuation was obtained at their oblique incidence on the 
screen. Direct experiments showed, that the greatest efficiency of the oblique screen is 
achieved if the screen has a convexity directed towards to the movement of a wave. One 
more practically important result is marked in the experiments: the bubbly screen placed 
directly nearby a rigid wall, is not equivalent by efficiency (i.e. does not result in the 
decrease of the wave amplitude) to the non-contact screen. The efficiency of the screen 


rises with its removal from a rigid wall. 


6.3. Regulation of the blast waves parameters at the charges detonation in 


various media. 


The calculated and experimental data on the intensity of blast waves generated by 
HE charges in homogeneous and heterogeneous media, allow to specify a range of 
expected overpressures at the transition from one medium to another [6.72]. The 
overpressures at the shock waves front from the HE charge (TNT) detonation in dense soil 
with the minimal number of gaseous inclusions (line 1) [6.76], in the gasless water (line 2) 


[6.81], and in the air (line 3) [6.39] are indicated in the summary diagram (fig. 6.3). All 


lines are plotted for the initial pressure Po=0,1 MPa. The addition of the gaseous gaps (a < 


10% vol.) leads to the sharp decrease of the blast waves amplitude. 


The shaded area of the overpressures (4) shows, that the level of the pressure can be 
lowered more than in 10 times when replacing of a dense ground by the porous filler. 
Transition from a continuous liquid to the system gas bubbles + liquid also results in the 
decrease of AP to a level typical for the pure gas (it is observed at the bubbles volume 
fraction @ > | %). The greater inclination of the dependences AP = AP (R«) in the gas 
bubbles + liquid system (also getting in the shaded zone (4)) testifies to the amplification 


of the dissipation processes. 


However, the dependence AP = AP (R«) for a gas (line 3) is not limiting. The 
smaller overpressures AP can be obtained at HE detonation in fine water foam with the 
concentration of the liquid o = (5...50) kg/m? (area 5). Let's pay attention to a number of 
specific features. The comparison of the dependence AP = AP (R«) for the air explosion (3) 
and the explosion in a porous ground (4) shows, that at R« = 4... 6 m/kg*? (R/Ro = 70... 
100) the amplitudes of quasi-acoustic pressure perturbations in the air and in the porous 
ground are almost identical. At R/Ro = 150... 200 the amplitudes of quasi-acoustic pressure 
perturbations in the air are higher, than in the porous ground. The final time required for 
completing the relaxation processes in two-phase medium, means, that in a direct vicinity 
of HE charge the impulse, as well as the overpressure in two-phase medium is higher, than 


in a gas. 


In the summary diagram the band of experimental values AP at the explosion in the 
foam is complemented by the calculated dependence (dashed line) AP = AP (R«) (taken 
from [6. 17]) for the water foam with the density o = 3...6 kg/m? in the far zone of the 
explosion at R/Rp > 15. The inclination of AP = AP (R«) for the explosion in the foam and 
in the air at R/Ro > 150 is approximately the same, i.e. the achieved level of the reduction 
of the blast wave amplitude is kept constant and is close to the limiting achievable one. 
The obtained feature of the blast waves attenuation in the foam and in the air specifies the 
relationship between the achievable level of attenuation and the ratio of the sizes of HE 


charge Ry and the foamy envelope A;. This relationship can be written as follows 


AP R 
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This dependence is valid for the envelopes with the sizes 10 <R/Ro < 26 at any 


distance from HE charge. 
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Fig. 6.3. The summary diagram of dependence of blast wave overpressure for various 
media environments: | - dense ground; 2 -gasless water; 3 — air; 4 — porous ground; 5 
— water-air foam (water-mechanical foam) 

The value of y** in dB does not depend on the distance from the center of explosion 
in a range 180 < R/Ry < 1800 or 8... 10 < R« < 10° m/kg®”’. The intersection of the line 3 
with the area 4 for the bulk systems shows, that at the distances from the center of 
explosion A/Ro > 150 the levels of blast loading can be diminished not only by placing the 
HE charge in a gas-liquid foamy envelope, but also in a protective envelope made from 
any granular or porous material. On the basis of generalization of the large amount of the 


experimental data B.I. Palamarchuk has derived the dependence [6.1]: 


y* = (32-4)42 for R/Ry > 150 (6.13) 
Here z = (oV/G) °°? or z=(m/G)°**, 

V — Volume of a protective envelope; 

G — Weight of HE charge; 

m — mass of a protective envelope. 


The given formula is suitable for the estimations of protective properties of 
envelopes made from the water foam, the snow pillow, protection made from the layers of 
the sand, the glass fiber, the metal wire. Let's emphasize, that y** > 0 at z > 4/3, i.e. at m/G 
> 2,4. We present some information from [6.1, 6.3, 6.4, 6.16, 6.39... 6.41, 6.48] about the 
parameters of the blast waves in the foam envelope itself, i.e. in the near zone of explosion. 
The dependence of the relative impulse of the reflected blat wave on relative distance in 
the air (1) and in the foam (2) is represented in the fig. 6.4. The impulse of the reflected 
wave in the foam at R < 0,64 m/kg?** is more, and at R> 0,64 m/kg?** is less, than in the 


air. 
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Fig. 6.4. Dependence of the scaled impulse of the reflected wave on relative distance: 1- 
air; 2 - foams at o = 5 keg/” 
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Fig. 6.5. Dependence of the blast wave overpressure on the concentration of the liquid at 
Rx = 0,85 m/kg?? 

For determining the influence of mass concentration of the liquid on the blast 
waves parameters in the foams the experiments with RDX charges [6.40] were carried out. 
The dependence of the blast wave overpressure on the mass concentration of the liquid is 
presented in the fig. 6.5 at Rx = 0,85 m/kg°*’. As it is seen, the absolute value of the 
overpressure at the front of the blast wave in the foam is very sensitive to the concentration 
of the liquid up to o = 4...5 kg/m’. Only slow reduce of the overpressure is marked with 


the further increase of the foam density. 


The relative decrease of the blast wave overpressure in the air and in the foam at the 
different scaled distances is shown in the fig. 6.6. The relationship between the attenuating 
factors for the blast wave in foam with the concentration of the liquid inside foamy layer 


can be approximated by the linear dependence: 


Ale = 140,750 


| foam 


(6.14) 


The formula (6.14) is valid at 2 kg/m’ < o< 20 kg/m’, 0,8 m/kg®*’ < Rs < 2,5 m/kg’’. The 


0,33 


attenuation by AP and / will occur at Rx > 1 m/kg’~”. It means, that for a charge of mass G 


= 1 kga foam layer of radius X = | m is required. 


All above conclusions are confirmed by the calculations and the experiments in 
[6.17]. The scheme of HE charge placement represented in the fig. 6.7 was studied: 
spherical HE charge (TNT) | was surrounded by the spherical foamy envelope 2. 
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Fig. 6.6. Dependence of the attenuating factor for the blast wave in the foam on the 
concentration of the liquid at different distances from the center of explosion 


Fig. 6.7. The scheme of a configuration explosive substance — gas-filled envelope: | — 
HE charge; 2 — water foam; 3 — microphone 
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Fig. 6.8. Dependence of overpressure in a blast wave on the distance for air and foams 
with various multifold: | a= 60; 2 a= 200; 3 a= 360 


Radius of the charge R = Ro, radius of the foamy envelope Ry Parameters of the 
shock waves were calculated for R/Ro = 40 at the foam multifold a= 60, 200, and 360 (that 
corresponds to the concentration of water o = 16,6; 5 and 2,78 kg/m’ accordingly). The 
results of calculations are presented in the fig. 6.8. At R/Ro < 10 the shock wave pressure in 
the foam exceeds the shock wave pressure in the gas. However at R/Ro > 10 the effective 
decrease of the overpressure is achieved, which increases with o changing from 2,78 kg/m* 


up to 16,6 kg/m’. 


In the fig. 6.9 the HE charge is separated from the air by a water layer 2. Water is 
kept in a thin aluminum envelope (3). In [6.55] the calculations of air shock wave 
amplitude were performed for the detonation of spherical TNT charge of 0,45 kg and the 


radius Ryo = 4 cm. 


The variants are investigated, when the aluminum envelope contacts the surface of 
HE charge or is separated from it by a water layer with the radii Ry/Ro equal to 1,5; 2; 3; 


and 4 accordingly. Thus the ratio of the mass of the water m, to the mass of the charge G 


was m,/G = 2,1; 5; 16,8; and 40 accordingly. The results of the comparison are presented 


in the fig. 6.10 in coordinates wave amplitude— relative distance. 


Fig. 6.9. The scheme of a charge placement in the water and in the envelope by [6.55]: 1 
— HE charge; 2 — water; 3 — envelope 
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Fig. 6.10. Dependence of overpressure on distance by [6.55]: 1 — in air; 2 — in thin 
aluminum envelope; 3 — in water + Al envelope 


As the scale for the comparison the radius of a charge Ro = 4 cm is taken. The line | 
is obtained at HE charges detonation without an envelope in the air. The line 2 corresponds 
to the calculation of AP = f{(R/Ro) at HE charges detonation in a thin aluminum envelope. 


The differences in the dependences | and 2 can be not considered in view of the 


measurement errors and accepted assumptions. The band of the values 3 corresponds to the 
detonation of a charge, immersed into the water. The bottom boundary of the area 3 
corresponds to the case m,,/G = 40, and the top boundary is responsible for the case m,/G 
= 2. At the distance R/Ry = 40... 50, i.e. at the distance 160... 200 cm, the distinction of all 
considered variants is negligible, i.e. the effect of attenuation of loading by the amplitude 
of a pressure wave is insignificant. Thus, the scheme of a charge placement, shown in the 
fig. 6.9, does not represent practical interest for the reduction of the blast loadings at HE 


detonation in the air. 


6.4. Attenuation of weak shock waves generated at the detonation of HE 


charges, placed in gas-contained envelopes 


From the ecological safety point of view [6.1,6.3] even the weak quasi-acoustic 
shock waves are the physical pollution of an environment [6.2,6.5]. In this connection the 
problem of creation of reliable and effective ways of SW attenuation [6.1, 6.2,6.3, 6.7] is 


urgent. 


As shown in research works of the Institute of Electric Welding IEW (Ukraine), at 
the welding and the cutting by the explosion in industrial conditions (and also near to 
inhabited buildings and structures) the methods of protection on the basis of multiphase 
media [6.4,6.40... 6.42] can be successively applied with the purpose of the ecological 
safety maintenance. Thus for the solution of a lot of practical tasks it is important to know 
not only the law of SW attenuation in two-phase mixture, but also the behavior of SW 
being weakened in the protective screens after their exit into the environment. Of interest is 
how the protective envelopes made from the gas-filled material decrease the SW 
overpressure, the impulse and the flow of energy [6.5,6.6,6.17, 6.25, 6.26,6.27,6.29... 
6.35]. 


Let's present the results of experimental studies in IEW, during which the laws of 
weak SW attenuation were determined at detonation of superimposed HE condensed 
charges. Also studied are the effects of gas-filled envelopes on the SW parameters. Water- 


mechanical foam and the snow also were used as protective tools [6.40... 6.42]. 


The efficiency of SW parameters reduction was determined from the comparison of 
the parameters of the blast waves generated at detonation of superimposed HE condensed 
charges and at detonation of HE charges in protective envelopes made from gas-contained 
materials. The SW parameters were measured with the help of the impulse noise-meters 
and the pressure gauges having resonant frequency 200 kHz, with the subsequent record on 
the registration system. Use of the specified methods of measurement has allowed 
expanding the frequency and the dynamic range of the registered overpressure and to 
receive the information on SW structure, transmitted through a protective screen 


[6.4,6.6,6.7,6.25... 6.27]. 


The experimental dependences from [6.40] for the SW overpressure and the scaled 
impulse of the compression phase on the scaled distance Rs=R/G°** are presented in the 
fig. 6.11 for the detonation of HE charges with the mass from 20 g up to 100 kg, where G 
is the TNT equivalent of a charge expressed in kg. The relationship between the 


overpressure at the front of SW and the scaled distance R« has a form 


AP=AR;” 
(6.15) 


Where A and m are the constants dependent on the properties of a surface, on which 


there is an explosion; Rs — in m/kg’**, AP — in Pa. 


At explosion of the superimposed (on a ground) charges on the sandy ground (A = 
5,510°, m = -1,6) the formula (6.15) gives rise to the AP values with an error 15... 20 %. 
For the SW of such intensity at the detonation of the freely suspended HE charges m ~ - 
1,13 [6.4,6.5]. The effect of the increase of the parameter m at detonation of the 
superimposed charges is explained in [6.4,6.5] by the energy losses as a result of SW 
interaction with the surface. Indirect confirmation of this conclusion was obtained in [6.4, 
6.5], namely the influence of surface properties on the rate of SW attenuation. So, at rather 
large distances from the charge (R> 150K o) the presence of a snow-cover results in the 


increase of m. 
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Fig. 6.11. Dependence of the parameters of the quasi-acoustic SW on the distance: 1 — 
peak pressure in SW after passage of snow-cover; 2 — the same with protection by 
multifold foam; 3 — impulse of a positive phase without protection; 4 — impulse in case of 
protective envelope made from multifold foam; 5 — peak pressure in SW at a snow 
surface; 6 — the same at a rigid ground [6.40] 


The performed studies [6.40] have shown, that the SW, generated in the air after the 
passage through the protective environment, are similar to the waves generated by the 
detonation of HE charges with smaller energy. The typical records of the quasi-acoustic 
SW are shown in the fig. 6.12. As it is seen, the pressure waves profiles slightly differ from 
each other, and the difference between their macroparameters is no more than 10... 15 %. 
So, the pressure impulse of a positive phase for an open charge is equal 2,8 Pa’s, AP = 1 


kPa; for the protected charge / = 3,0 Pas, AP = | kPa. 


Fig. 6.12. Profiles of pressure for quasi-acoustic shock waves at the distance 30 m from 
RDX charges with mass 0,15 kg in air (a) and I kg - in protection from snow (b) by the 
data of measurements from [6.40] 


Following conclusions of the Institute of Electric Welding [6.40, 6.42], the 
macroparameters of the SW generated by HE detonation in the gas-filled envelopes can be 
determined by TNT equivalent G* of the explosion. The value of G* is equal to the energy 
of the unprotected HE charge, at which detonation at the fixed distance from the place of 
explosion the parameters of SW coincide with the parameters of the blast wave generated 
by the charge in protective envelope. The efficiency of the protection in this case can be 


defined by the factor of reduction of the energy of explosion 77: 


G*=nG 
(6.16) 


In the case presented in the fig. 6.12 7 = 0,15 both for overpressure and impulse. 
One of the parameters determining the damping properties of gas-contained medium, used 
for the SW attenuation, is the mass concentration of the condensed phase o. Changes of the 


wave duration were not revealed, as well as in [6.17]. 


At the fixed volume of protective envelope the efficiency of decreasing the SW 
parameters by multifold foams with o = 1,5... 2 kg/m? has appeared small in comparison 


with of the snow or the foams with o= 10... 20 kg/m’. 


For the SW parameter predicting it is necessary to know the relationship between 
the factor of attenuation and the properties of the medium used for localization of 
explosion, and geometrical sizes of protection screen. B.I. Palamarchuk analyzed the data 
[6.3,6.6] and established the criteria of SW attenuation at explosion in gas-contained 


environment using dimensionless parameter Z 


Z=(oVs/G)°*? 
(6.17) 


Where ois the mass concentration of the condensed phase in protective medium, 


kg/m’; 
Vs; is the volume of a protective medium, m*: 
G is the TNT equivalent of the explosion, kg TNT. 


The dependence of value AL, describing the decrease of the overpressure at the 
front of SW, on dimensionless value Z is shown in the fig.6.13 for the case of protection 
based on gas-contained media. The shaded area corresponds to the discrepancy of 
experimental data for such media as foams and snow. The same area contains the data of 
experiments on the attenuation of SW, transmitted through protective envelope from the 
fiber glass, the metal cotton wool, the foams obtained in [6.6]. The dependence presented 


in fig. 6.13 can be described by the formula: 


AL =3Z-4= 201g (AP, /AP 2), 
(6.18) 


Where AL is the difference between the level of sound pressure, dB, at the fixed 


distance from HE detonating charge without (AP) and with (AP?) protection. 
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Fig. 6.13. Dependence of SW damping efficiency on the properties of protective envelope 
at O, kg/m’: x - water foam (o = 6... 10); o- snow (o = 250... 500); — - sand (o = 1600) 


To determine the factor of explosion energy decrease we rearrange the equations 
(6.16... 6.18) in the form: 


7 = G/G* = 10° AL/m, m ~ 1. 
(6.19) 


The presented formulas (6.18 and 6.19) enable to estimate the efficiency of 
explosion localization by gas-contained media and to optimize the parameters of protection 
of technological metal processing by the explosion. For the practical applications of 
interest is to describe the experimental setup in [6.6]. The fig. 6.14 represents the 
container, in which the detonations of TNT charges were made (the masses of charges 


were 0,11 kg, 0,57 kg and 2,27 kg). 


Fig. 6.14. The scheme of experiment for determination of blast waves attenuation 
efficiency by foam: | - wooden frame; 2 - polyethylene film; 3 - HE; axbac~ (0,9... 1,2) m 
[6.6] 


Fig. 6.15. The scheme of microphones positions in experiments [6.6]: 1 - container; M), 
M>, M3, Mg, - microphones 

The HE charges were placed in the water foam or were buried by porous materials. 
The multifold of the foam was 30... 250, that corresponds to the concentration of two- 
phase medium, and equals accordingly o = 10... 4 kg/m’. The characteristics of other 
materials forming the two-phase medium are given in tab. 6.1. Two types of steel fabric, 
two types of gas balloons, two types of glass-wool, and straw were used. The experiments 
with the filling of the container by water-mechanical foam were done at the same setup. 
The measurements were performed by the scheme presented in the fig. 6.15, where the HE 
charge was in the polyethylene container, and M), M>, M3, Mz; are the pressure transducers 


(microphones). 


The results of measurements are shown in the fig. 6.16. Here 
1=0,5(abc)?=0,5??, 


a, b, and c are the linear sizes of the container, V is the volume of the container. 


Besides, p, = o is the volumetric density of two-phase filler, so ,V is the mass of the 


condensed absorber in the container, G is the mass of HE charge. 
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Fig. 6.16. Reduction of the sound pressure level at HE detonation in porous medium, O - 
steel fabric (thin), W - steel fabric (rough), 0 - gas balloons (small), @ - gas balloons 
(large), A - straw, © - glass-wool (thin), <I> glass-wool (thick) 


Table 6./ Parameters of materials for explosion damping 


Material Diameter, | Density, | Thermal Specific heat 
um kg/m? conductivity, capacity, J/kg K 


Steel fabric, fine 150 7750 
Steel fabric, rough 240 7750 
Air balloons, small 14500 Co 1430 


Air balloons, large 29000 } 9300: F038 1430 
50 0,039 1340 


Glass-wool, rough 10 2500 0,036 


The comparison of the shock waves attenuation data at the charges placement in 
explosion-absorption gas - containing envelopes shows, that the envelopes made from the 
water-mechanical foam, the snow and other elementary materials have the largest 
efficiency. The application of the envelopes made from the materials with large structure 
density does not give rise to appreciable advantages in efficiency, but introduces the 
serious technical difficulties at realization. Basic difficulty is connected with the 
inconveniences of formation of a heavy explosion-protective envelope around the 


explosive object, (as a rule, it is the HE charge with detonator). 


6.5. Basic directions of explosion protection 


It is possible to distinguish two ways of the explosion protection methods. The first 
direction is caused by the obvious attempt of HE charge shielding with the help of a firm 
impenetrable envelope. It is enough to mention a trivial case of placing HE charge in the 
closed or half-closed container [6.28]. The practice showed the insufficient efficiency of 
this way of localization. For improvement of operational properties the internal cavity of 
the container was filled by compressible medium. In [6.58,6.59] the liquid foam was used 
as filler, in [6.60] — liquid with gas bubbles, in [6.6] — granular or fibrous material. Each 
of the listed explosion absorbers can be characterized by the effective density o, differing 
from the density of the frame material . Thus, for the water foam o= 10... 30 kg/m’ at p= 
1000 kg/m’, for the water with gas bubbles o = 900. . 990 kg/m’ and for the fiber and 
granular fillers o> 1000 kg/m’? at 9 > 2000 kg/m’. 


The second direction of the explosion protection is characterized by principal 
decision to refuse of application of an impenetrable firm environment around the charge. 
As limiting here it is possible to examine the variant of the blast wave suppression in the 
volume of the sprayed liquid (more often water) [6.62, 6.75]. However, the impossibility to 
obtain the spray with high concentration of a liquid o> 1... 2 kg/m’ makes this way of 


attenuation inapplicable. 


To keep the concentration of explosion absorbing material guaranteeing reliable 
attenuation of the blast waves, the gas-contained covers made from water foam [6.58], 
polyurethane foam [6.64] were offered. For the water foam o = 10... 20 kg/m’, for the 


polyurethane foam o= 50... 100 kg/m? at p= 900 kg/m’. 


The practice of the application of the gas-contained envelopes for localization of 
the demolition loadings for HE charges up to 3 kg showed, that the basic part of the 
explosion products energy corresponds to the kinematics acceleration of a mass of 


protective envelopes [6.66]. 


Also essential is the transformation of the energy flows at the boundaries: products 
of explosion — explosion absorption medium — air [6.65]. Quite ineffective were the 
attempts of explosive loadings regulation by the energy losses due to the processes of 


evaporation or dispersion of the protective envelope material [6.86... 6.88]. 


In view of a dominant role of inertial properties of the protective envelopes for the 
blast waves attenuation, an effective way of explosion suppression is the change of the 
kinematics characteristics of the medium, which damps the expansion of explosion 
products. The major characteristic is the compressibility of the medium, surrounding a 
charge. The reliable suppression of blast loadings is reached, when the protective 
environment is the liquid guaranteeing (at least at early stages of expansion) against the 


outflow of explosion products into the protecting atmosphere. 


The measure of the compressibility of a liquid is the speed of sound. The regulation 
of the sound speed in the liquid can be achieved by the introduction of the gas bubbles 
[6.56,6.63]. The introduction of compressed inclusions in a liquid volume essentially 
influences the propagation of the blast waves of moderate intensity with the overpressure 
P,/Po < 1000. For the waves with larger intensity the presence of compressible elements is 
not so essential. The immersing of HE charges into compressible two-phase medium 
repeatedly changes the characteristic temporary scales of all wave processes and expands 
the opportunities for an effective exchange of energy between explosion products and 


explosion protective envelope [6.48]. 


The listed circumstances have allowed offering the new scheme of explosion 
suppression. Fig. 6.17a presents such a scheme for HE charge with the help of a layer of 


liquid to be dispersed, placed in volume between elastic envelopes [6.61]. 


The layer of a liquid is separated from HE charge by gaseous gap. The arrangement 
of the basic elements is different from the traditional scheme [6.66], shown in the fig. 6.17 
b. The basic role of gaseous gap consists in attenuation of the blast wave amplitude before 
its incidence with a liquid layer due to distinction of volumes of a charge and internal 


cavity of an envelope. 


Fig. 6.17. The possible schemes of explosion suppression with the help of liquid elastic 
containers. 1- gas (gap), 2- envelope (plastic), 3- liquid to be dispersed, 4- HE. 


On the basis of the known [6.66] and new experimental data it is possible to 
compare the efficiency of the air blast waves attenuation using the schemes shown in the 
fig. 6.17, a, b. The information on the blast waves attenuation with the help of a water 
layer around the HE charge, 1.e. by the scheme in the fig. 6.17 5, was taken from [6.66]. 
The information on the blast waves attenuation by the scheme in the fig. 6.17 a was 
obtained in a series of experiments, similar to [6.57]. The mass ratio of the water M to the 


HE charge G was equal to 50, G = 40 g in [6.56] and G = 100 g in [6.57]. 
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Fig. 6.18. Dependence of the attenuation factor for the blast wave by amplitude on the 
scaled distance for the various schemes HE detonation 
Le, kPa O33 
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Fig. 6.19. Dependence of the pressure impulse in a phase of compression on the scaled 
distance for the various schemes of HE detonation 


The change of blast wave attenuation factor G by amplitude is shown in the fig. 
6.18 depending on the scaled distance R:=R/G°”*. Distance R was counted from the center 
of explosion up to the place of measurement. Factor of attenuation 6 =AP(R+)/AP)(R:), 
where AP(R:) is the amplitude of the air wave at the explosion in the liquid envelope and 
APo{R+) is the amplitude at the explosion in the gas-filled envelope of the same size by the 
scheme in the fig. 6.17 c. It was made because the value AP» (shaded strip 1 in the fig. 


6.18) is 10... 20 % less, than AP;(R«) at a completely open charge. 


The strip 2 in the fig. 6.18 is obtained for the explosion by the scheme in the fig. 
6.17 5, and strip 3 is obtained for the explosion by the scheme in the fig. 6.17 a. 


Fig. 6.19 demonstrates the change of the impulse /* in the compression phase with 
the distance R» at the detonation of HE charge in the gas-filled air envelope (line 1 by 
[6.66]). The strip 2 for /* is obtained for the explosion by the scheme in the fig. 6.17 4, and 


strip 3 is obtained for the explosion by the scheme in the fig. 6.17 a. 


The diagrams (fig. 6.18 and 6.19) confirm the highest efficiency of the scheme 


shown in the fig. 6.17 a for the reduction of demolition consequences of explosion. 


The carried out studies are based on the rational use of the specific features of 
pressure and rarefaction waves propagation at HE detonation in a medium with elevated 


compressibility in comparison both with the liquid and the gas. 


Essential difference of the suggested scheme of explosion suppression is the 
preliminary correction of initial conditions of expansion of explosion products with the 


help of a gas layer between the HE charge surface and the compressible medium. 


6.6. Comparative characteristics of two ways of the air shock waves 


attenuation 


Because of various types of blasting operations and the design of the explosion 
protection structures, there are many works published concerning the problem of 


decreasing the overpressure and impulse of the air shock waves. 


As an example it is possible to list the cited above works on the theoretical and 
experimental study of the shock waves propagation in channels through the perforated 
screens, works on the suppression of the air shock waves at the spherical explosion of HE 
charges when placing the charges in the volume of the water-mechanical foam or in front 
of the perforated screens. Unfortunately, there are no reliable calculation models, capable 
to predict in advance the efficiency of particular method of the shock waves attenuation. 
Therefore when selecting the methods it is necessary to use the experimental data. A 
choice of the way suppression of air shock waves is difficult, because the developed 


methods and the approaches are not universal. 


Moreover, the transferring of results from one area of parameters on another can 
result in serious mistakes and obtaining the negative result. Let's present some examples 
based on the measurement of the shock waves parameters at HE detonation in the water- 
mechanical foam. As it was already mentioned, the mass of HE charges (TNT) in 
experiments was varied from 0,01 up to 3 kg. The density of th water-mechanical foam 
was o< 10 kg/m’. The dependence of the factor of shock waves attenuation by pressure 7 
=AP/AP» on the distance R/Ro is shown in the fig. 6.20. Here AP» is the shock wave 
intensity in the air at the distance RX, and APo is the overpressure in the water-mechanical 
foam at the same distance, Ro is the radius of HE charge. It is seen, that at R/Ro > 10 
AP>APo, but at R/Ro < 10 AP>APo. Thus, the zone of real attenuation of the shock waves 
is at rather large distance from a charge. The dependence of the factor of attenuation on the 
overpressure at the SW front is shown in the fig. 6.21. As it is seen, only in the far zone of 
explosion the limiting attenuation of the shock waves is reached # = 7...10. The 


overpressures on the abscissa axis correspond to the detonation of HE charge in the air. 


Fig. 6.20. Dependence of the factor of attenuation (amplification) by pressure on the 
scaled distance to the center of explosion (1- zone of amplification; 2 — zone of 
attenuation) 
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Fig. 6.21. Dependence of the factor of attenuation (amplification) by pressure on the 
overpressure in the front of a wave (1- zone of amplification; 2 — zone of attenuation) 
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Fig. 6.22. Dependence of the factor of attenuation n (amplification Kamp) by impulse on the 
overpressure in the front of a wave (1- zone of amplification; 2 — zone of attenuation) 


Though the information on the impulses of a compression phase in the air wave in 
porous medium is very limited, nevertheless, on the basis of some measurements it is 
possible as a first approximation to construct the dependence of the factor of attenuation by 
impulse 7 on the intensity of a shock wave. Such dependence is presented in the fig. 6.22. 
As it is seen, only at AP < 1 MPa, i.e. at the distance from the charge R/Ry > 20, the 
reduction of the impulse becomes appreciable. The reduction of the impulse in a pressure 
wave when the HE charge is immersed in the foam does not exceed the values 77 = 1,3... 


1,5 at R/Ro < 20, but will be increased with the distance. 


According to [6.8,6.43... 6.47,6.52... 6.54] another effective way of the shock 
waves attenuation is the installation along the waves path of the perforated partitions of 
various types (which are described in [6.43... 6.47]). The generalization of the 
investigations on the shock waves propagation through the perforated partitions has 
allowed obtaining some empirical dependences. According to [6.8] the pressure in an air 


wave behind a partition is determined by the expression 


AP =TAIZ (Ry ae) 220%: (6.21) 


Here Z = R/G°**, R*® = R/L, AP=P|-Po in MPa, P is the pressure at the wave front, 
G is the mass of a charge in kg, R is the distance in m, L is the characteristic size of the 
partition. The formula (6.21) is valid at 1,16 < Z < 8,45 m/kg®**, 0,015 < AP < 1 MPa, 0,69 
< R*< 4,55 and permeability of the partition 0,01 < a@ = Fo/F < 0,13. In the expression for 
a, there are designations: Fy - the total area of apertures, / — the area of cross-section of a 
protective partition. In [6.53] the expression for calculation the intensity of a weak wave, 


© Tt is rather close to 


passed through the perforated screen is presented in a form AP*( a) 
the conclusions from [6.51]. The change of the compression phase impulse can be 


determined by the equation for the estimation the impulse behind a barrier 


1G? =791Z°"* a" £19 %. 
(6.22) 


Here / in Pas. According to [6.51] the formula for calculation the impulse is valid 


at 1,16 <Z<5,95 m/kg’**, 0,03 < AP < 1 MPa, and 0,008<a, < 0,13. 


The analysis of the presented relationship shows, that for the passing waves the 


°° < 0.4. At the reverse 


dependence predicts the attenuation by impulse only at Z°”*a, 
ratio the dependence specifies the amplification of the wave transmitted through the screen 
by impulse. One of the possible reasons of the phenomenon observed is the probability of 
registration of the impulse from the waves, which have rounded the perforated plate. In this 
connection the performance of the measurements of impulse in conditions excluding the 
occutrence of diffracted perturbations is necessary. The estimations of the impulse 
attenuation to be presented below are made in the area, where the decrease of the wave 
impulse is marked. In view of these remarks the given dependences allow to determine the 
factors of attenuation for the shock waves by pressure f and impulse 77 due to the partitions 
with different degree of permeability. The dependencies of @ on the wave intensity are 
shown in the fig. 6.23 at a = 0.01 and 0.1 (the databases in zones I and II accordingly). 
The shaded zone reflects the distinction of # for the barrier with k* = 4,5 (top border of a 


zone) and with k*= 1,16 (bottom border of a zone). 


The dependences of factors of attenuation by impulse on the wave intensity are 
specified in the fig. 6.24 at a = 0.01 (line I) and @ = 0.04 (line II). The lines II in the fig. 
6.23 and 6.24 correspond to the immersing of HE charges in the foam. As it is seen, at @& = 


0.01 the perforated partitions are more effective. 
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Fig. 6.23. Dependence of the factor of attenuation by pressure on the overpressure in the 
front of a wave 
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Fig. 6.24. Dependence of the factor of attenuation by the impulse of blast loading on the 
overpressure in the front of a wave 


The Fig. 6.25 allows determining the level of permeability, at which the attenuation 
of the air shock waves by the partitions is more preferable, than at the HE charges 
immersing in volume of foam. The dependences I and II show the level of attenuation by 
pressure at different permeabilities @ and the wave intensity AP. The zone III corresponds 
to the shock waves attenuation by the foam at the concentration o= 10... 15 kg/m’. At a < 
0.02...0.03 the partitions are more effective, and at a > 0.04 the foam screens are more 


effective. 


As it is seen, the effective attenuation by the perforated partitions is achieved at 
their small permeability. Taking into account the information presented before, the 
pressure in the reflected wave at such partition is close the pressure of reflection at 
continuous wall. It means, that the partition will be subjected to the rather appreciable 


shock loadings, which are necessary to be parried by the appropriate technical decisions. 


0,06 0,08 0,1 


Fig. 6.25. Dependence of the factor of attenuation by pressure on the permeability of a 
partition 


6.7. Efficiency of the combined protection 


The reasonable choice of an acceptable combination of the protection methods from 
fragments and blast loadings can be made on the basis of the experimental data [6.59]. In 
this work the efficiency of the combined protective surrounding (layers of foam + metal 
envelope) and the efficiency of the foam layer without a metal envelope were compared. 
The first configuration is represented in the fig. 6.26, and second one is represented in the 
fig. 6.14. The basic geometrical sizes of the envelope 3 are the same, as in the fig. 4.1 and 


4.2. 


Fig. 6.26. The combined protective configuration: layers of foam + metal envelope: | - 
foam; 2 — HE charge; 3 - metal envelope; 4 - ground 

The HE charges 2 (TNT/RDX 9/91 — TNT equivalent is equal to 1,36) were placed 
at the height of 4 cm above the ground. The envelope was putted 6 cm deeper in the 
ground for prevention of the shock waves propagation through the gap between the 
envelope and the ground. All measurements were performed with the foamy layer 1, 
having the multifold ratio 30, i.e. the density about 33 kg/m’. Metal envelope 3 is made 
from a part of the cylindrical pipe with the diameter d and the height 4. We shall remind, 
that the reduction of the sound pressure at the explosion in the empty container (without 


the foam) does not exceed AP = | dB. 
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Fig. 6.27 Comparison of efficiency of various protective envelopes by sound pressure 
Z=(oVs/G)°” 

The comparison of the protective envelopes efficiency by sound pressure is 
presented in the fig. 6.27, where the line I corresponds to the explosion by the scheme 
shown in the fig. 6.26, and line II corresponds to the explosion by the scheme shown in the 
fig. 6.14. Evidently, that the combined protection (fig. 6.26) provides an effective 


attenuation of the shock waves. 


The scheme is simple for the realization and can be recommended for the practical 
use at the detonation of the charges with mass up to 1... 2 kg TNT. The preferable 
geometrical sizes of the container: d = 0,92 m; / = 1,22 m and d= 0,86 m; # = 1,16 m, 1. 


the volume V = 1 m’*. 


6.8. Attenuation of the air blast waves from HE detonation of charges in a 


liquid limited elastic envelope 


The obtained specific features of the blast attenuation at HE detonations by the gas- 
contained envelopes give additional opportunities to develop new types of the accident 


protection devices. It is urgent, as the operative application of the foamy envelopes is 


complicated sedimentation of the foam, and the preparation of the protective envelopes 


with other fillers is connected with a number of technical difficulties. 


Let's consider the problem on the increasing the compressibility of a liquid (for 
example, water) by capsulation in an elastic compressible envelope. The solution of this 
task is known, and as an example we shall consider the medium water + air bubbles in an 
impenetrable envelope (for example, in a tube) with flexible walls [6.56,6.63]. The sound 


speed in such medium can be determined by the relationship: 


-l 

, la 1-a@ 20r 

om =| ap, +(1-a@) p, | zs Pe + | (6.23) 
jie a 


Here @ is the gas volume fraction in a liquid; p, is the density of a gas; ~ is the 
density of a liquid; P is the pressure; / is the modulus of elasticity of the envelope material 
(wall of a pipe), © is the factor describing the conditions of attaching the envelope (O = 
1... 0,5 for a tube fixed at one single end, and © = | - s/ for a pipe jammed at both ends, 
Lis the Poisson coefficient), 6 is the thickness of the wall of an envelope, r is the radius of 


an envelope. 


Parametrical dependence c=f(@) is shown in the fig. 6.28 for representing the 
practical interesting case of the medium: water + air bubbles at P = 10° Pa and the 
parameter 207/6 = 10. The line | is calculated for the sound speed in the incompressible 
liquid (cy 0), placed in an incompressible envelope (E +). The line 2 is obtained for 
the water (cy= 1450... 1500 m/s) at E +. The line 3 is obtained for the water in a steel 
pipe at / = 200 GPa. The line 4 is obtained for the water in the envelope made from an 
elastic material with E = 1 GPa. As it is seen, at @< 10” the sound speed c is adjusted by 
the elasticity of the walls of an envelope. At a volume fraction of the gas inclusions in a 
liquid a > 10° or, more precisely at @ > AP/pecf the compressibility of the structure: 
envelope — liquid — gas bubbles does not depend on the properties of the filler and the 
material of the walls. Here AP is the amplitude of the pressure perturbation. The diagram 
presented in the fig. 6.28 allows to understand the basic features of the suppression of the 


HE generated blast waves by the envelopes considered in [6.1, 6.3, 6.4, 6.6, 6.17, 6.23, 


6.26... 6.30, 6.40... 6.42, 6.48]: 1) at the gas volume fraction in a liquid a@ > 10° the effect 
of the properties of a liquid and the parameters of the envelope material, transmitting a 
pressure wave, in definition of medium compressibility, is insignificant; 2) high 
compressibility of transmitting medium can be achieved not only by entering the gas 
inclusions into water, but also at ~@> 107 by the water capsulation into an elastic envelope 


with E < 10 GPa. 
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Fig. 6.28. Dependence of the sound speed in the system liquid - air bubbles — envelope on 
the volume fraction of gas bubbles [6.56,6.57,6.63] 


6.9. Checking of the efficiency of the blast loading reduction 


Study of the efficiency of the blast effects decreasing by placing the HE charges in 
a volume of the liquid limited by an elastic envelope, was carried out in a series of polygon 
tests. The parameters of the air blast waves were measured at the ground explosion of 
spherical TNT charges with mass 0,1 and | kg. To exclude the suppression of the blast 
waves due to the energy losses by the interaction with ground (formation of a crater, 


transfer of a part of energy by the wave inside the ground) the HE charges were mounted 


on the armour-steel plate with the area 1,2 x 1,2 m’. The thickness of the plate in the 


experiments changed from 30 to 60 mm. All explosions were executed in the open district. 


The parameters of air blast waves were measured by six piezoelectric pressure 
gauges located at certain distances from the charge RX in two perpendicular directions. In 
the case of 0,1 kg charges the gauges were established at the height of 0,1 m above the 
ground. In the case of 1 kg charges the gauges were placed at the height of 1 m to exclude 
the influence of the reflected waves. The distance of the measurement of the blast waves 
parameters was in a range 0,8...4,0 m. In the scaled variables of Sadovsky-Hopkinson 
R*=R (2G)? for the charges of 0,1 kg at the ground explosion 1,36 m/kg’**<R*<6,8 
m/kg°** and for the charges of 1 kg 0,63 m/kg?**<R*<3,2 m/kg””? accordingly. 


The results of the measurements of the pressure waves amplitudes from the open 
charges have coincided (with accuracy not less than 10 %) with the empirical values 
calculated by the Sadovsky formulas [2.1], used for the further comparison to the 
parameters of the waves from the protected charges. The capsulation of charges was 
achieved by their placing into a double cylindrical elastic envelope filled by the working 
liquid. The size of internal cylindrical capsule with the volume about 2,6 1: a diameter 0,15 
m and height 0,15m. Internal capsule is surrounded by a ring layer with thickness 0,04 m 
and is limited by the cover with the same thickness from upper side (double envelope with 


a liquid). Thus, the total volume of liquid filler is about 5,2 1. 


The explosion suppressors (envelopes) for HE charges with the mass 0,1 kg TNT 
were made from sheet polythene with the thickness 6 = 2 mm and with limiting axial 
lengthening ~ 400%. The mass of an envelope was 0,6 kg. A number of experiments with 
1 kg TNT charges was performed in a rubber envelope (the thickness of walls is 1,2 mm) 
filled by 50 kg of water. Limiting axial lengthening of the rubber was about 400 %. The 
restriction of the elasticity of the rubber envelope was reached at its covering by hardening 
polyurethane foam layer outside. The information on the explosion suppressors used is 


listed in tab. 6.2. 


The designations in tab. 6.2 are: cy— sound speed in the filler; ~- absolute density 


of liquid filler; (cp) - acoustic resistance of the filler; 14 — mass of the filler; (M/G)°”? - 


parameter of suppressor; 


y- factor of blast waves attenuation by amplitude 


AP 
= 20lg —-; 
a SAP. 


2 


£- factor of blast waves attenuation calculated by (6.18) 


No special measures on the degasation were undertaken while filling the envelopes 
by the liquid. The following working liquids are taken: the water, the water solution of 
calcium chloride, the glycerin. The solution of calctum chloride was used in order to 
change the density of the liquid [6.85]. Glycerin was taken as a sample of high-viscous 


liquid [6.88]. A number of experiments were performed with the sand as a filler [6.70]. 


Table 6.2 Parameters of explosion suppressors 


No | Material | Filler G 2p (cp), | (W/G)” | y, dB | B B' 
of m/s | kg/m? | 10° - 
envelope kg/m” 
S 

Polyethyle | Solution 2:09: 0,31...0,48 
ne of 25 

calcium 

chloride 
Polyethyle | Glycerin 1,96.. 0,32...0,51 
ne 1 
Polyethyle | Water 1,81.. 0,35 .....0,55 
ne 87 
Polyethyle Ps 0,26...0,42 
ne 12 
Rubber 1,79.. 0,35...0,56 

83 

Rubber + a bs Se 0,35...0,56 


polyuretha 


ne foam 


6.10. Influence of an elastic envelope on the parameters of air shock wave 


The general representation about the change of blast waves parameters at HE 


charges immersing in gas-filled envelope or in a liquid with an elastic envelope is shown in 


the fig. 6.29. The line | and the experimental points are obtained for the open charges in 


calculations [2.1] and in experiments [6.80] accordingly. The line 2 was obtained in 


[6.3,6.48] for TNT detonation in the water-air foam at the concentration of water o = 10 


kg/m*. The area 3 corresponds to the described experiments with HE charges, shielded by a 


liquid layer in elastic envelopes [6.94], and data [6.57]. As it is seen, the efficiency of the 


shock waves attenuation by the foam or by the liquid in an elastic envelope appears to be 


comparable. 


The effect of various explosion suppressors on the dependence of blast waves 
amplitude on distance is shown in the fig. 6.30. All liquid suppressors have approximately 
identical (within the experimental error + 10 %) efficiency. In the experiments with the 
sand the least degree of blast waves attenuation was obtained. Smaller compressibility of 
glycerin in comparison with the water (the speed of a sound for glycerin is 1920 m/s and 
1450 m/s for the water) compensates the influence of larger density [6.88]. The diagram of 


dependence of the attenuation factor for blast waves by amplitude on the distance 


6 HAP AP HGS (R’) is presented in the fig. 6.31. At HE detonation in the volume of 


a liquid the value f increases with distance in a range up to R* < 2,5 m/kg®*’. The 
comparison between calculated values (by formula for yw for gas-filled envelopes) and 
experimental data of limiting level of attenuation by amplitude fis presented in tab.6.2. It 
is possible to note the agreement between calculated values of / (taking into account the 
experimental error) and the values f. obtained in experiments. The results of 
measurements of amplitude of waves from explosion of 1 kg HE charge shielded by a 
water layer in hard (line 3) and not- hard (line 4) rubber envelopes are presented in the fig. 
6.32 in comparison with the data for an open charge [6.12]. A conclusion about the 
influence of envelope elasticity on the air waves formation partially proves to be true. It is 
useful to add the fig. 6.32 by the calculation results [6.55] for a system: HE charge (mass 
0,45 kg) + water bath + thin envelope from an aluminum alloy (/ ~ 70 GPa) - line 2. No 
appreciable decrease of the blast waves amplitude was found in calculations [6.55] for the 
ratio of masses of water and HE charge from 2,1 up to 16 and at R/Ro > 30. It once again 
reflects the appreciable influence of envelope elasticity at a choice of effective dampers of 
the blast loadings. The formulas for calculations of the attenuation degree ( by amplitude 
for gas-filled envelopes are inapplicable to the case described in [6.55], when the HE 


charge was immersed in a liquid and placed in a metal envelope. 
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Fig. 6.29. Dependence of air shock wave amplitude on scaled distance. Line | - by [2.1]; 
line 2 - by [6.48]. Points — by the measurements for an open charge [6.57]. The shaded 
area - for charges in elastic liquid explosion suppressors 
APs, kPa 
120 


100 
Bt 


60 f- 


0,7 10 13 1,6 smn 
1,20 1,71 222 2,74 R* puke? 


Fig. 6.30. Dependence of air shock wave amplitude on distance for various variants of 
explosion of 100 g TNT charge. Numbers of curves correspond to numbers of explosion 
suppressors listed in tab. 6.2 
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Fig. 6.31. Dependence of air shock wave attenuation factor on distance at ground 
explosion of 100g TNT charge. Numbers of curves see fig. 6.30 


To obtain the complete pattern of the blast waves attenuation from HE detonation, 
shielded by the protective liquid-elastic container, it is necessary to perform the 
measurements of the impulse of the static pressure and the duration of the compression 
phase. The construction of pressure - impulse diagrams, as well as the description of the 
blast wave profile transformation in coordinates pressure — time at any distance from the 
charge in comparison with the standard profiles for an open charge, are necessary for the 


further development of the methods of protection [6.67... 6.70]. 
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Fig. 6.32. Dependence of air shock wave amplitude on distance at ground explosion of 1 
kg TNT charge: 1 - by [2.1]; 2- by [6.55]; 3 - experiments with explosion suppressor No 6 
(tab. 6.2); 4 - experiments with explosion suppressor No 5. 

The results of the comparison between the pressure at the front of an air SW at free 
and localized explosions of TNT charge with the mass 0,8 kg are shown in the fig. 6.33. 
The lines of the pressure levels of 50 and 100 kPa, corresponding to the thresholds of ear 
damage and lungs damage, show the reduction of radius of the dangerous zone of 
demolition action in 2...2,5 times as a result of explosion localization. The experimental 


values of the impulses for free and localized 1 kg TNT are presented in the fig. 6.34. 


The obtained dependences of pressure and impulse on the distance and the mass of 
the open and localized TNT charges (the mass < | kg) allow to estimate the blast effect of 
the explosion, for example, by the construction of the diagrams in AP-/ coordinates. The 


example of the diagram is presented in the fig. 6.35. 
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Fig. 6.33. Pressure at the front of air SW depending on the distance at explosion of open 
and localized 0,8 kg TNT charges: 1 — lethality threshold; 2 -barotraumas threshold; 3 - 
open charge; 4 - localized charge 
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Fig. 6.34. SW impulse depending on distance at kg TNT explosion: | - open charge, 2 and 
3 - charge, localized by the 50 kg samples with various elasticity 
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Fig. 6.35. P —I diagrams for the open and localized TNT charges with masses 0,1... 1,0 
kg. Experimental points for R = 0,8m (Gryr = 0,1 kg) and R = 1,5m (Grnr = 0,4; 0,8; 
1,0m).: ° - data for the localized charges; J - data for the open charges 


The presented dependences show, that the blast effect of explosion exceeding a 
threshold of lungs damage in case of the open charges is reduced down to the safe level 


with the use of localizators. 


6.11. Transformation of the blast waves by the layer of two-phase medium 


Experimentally is proved [6.70, 6.71], that at the explosion of the isolated HE 
charge the blast wave is attenuated at crossing the liquid screen due to the loss of energy 


on the destruction of localizator and the kinetic energy of the liquid to be dispersed. 


The blast wave means the pressure perturbation with a_ continuously 
decreasing profile. Due to the energy loss and the transformation of the blast wave in 
the layered screen the amplitude of an air blast wave nearby localizator is reduced in 
several times in comparison with the explosion of the non-isolated charge. It seems, 
that the influence of loalizator is reflected not only in the peak parameters of the blast 
wave, but also on the frequency - temporary characteristics. The question on the 
transformation of the blast wave by its peak - frequency parameters under the 


influence of the gas-liquid screen is not completely resolved till now. 


The studies in [6.57,6.70,6.79] have revealed the main role of the compressibility of 
the liquid working fluid in localizator in the blast wave attenuation and transformation. 
Efficiency of the blast waves suppression is increased considerably at the replacement of a 
continuous liquid by the gas-liquid system. Such replacement has ensured the increase of 
the compressibility of the medium, transmitting the blast energy. Partial, but the not 
limiting increase of the compressibility of the liquid screen was achieved in [6.57] by the 
elastic envelope holding a liquid. Additional interest is the specific features of the 
processes accompanying the blast waves transformation at passage through the layer of a 
two-phase protective material. The introduction gas bubbles in a liquid changes 
compressibility of a liquid and drastically changes the speed of sound and shock-wave 
disturbances [6.73... 6.76]. In two-phase medium a liquid + gas inclusions the speed of a 


sound depends on the pressure, and it provides the effective blast wave attenuation. 


The multifold reduction of the speed of linear and nonlinear pressure disturbances 
in two-phase medium influences a structure of the pressure in the blast loading due to the 
stretching of the perturbation in time. The previous studies are limited to the shock waves 
propagation in gas-liquid media. The known investigations concern either to the very 
strong waves with the amplitude of several GPa [6.73], or to the rather weak waves with 
the amplitude up to 1 MPa [6.74, 6.77,6.78]. The duration of the compression phase was 
those that the statement of a question about the frequency characteristics had not sense. 
The amplitude of the pressure in the blast waves arising in gas-liquid localozators 
preventing from a HE charge detonation, typical for the explosive devices, was 0,1... 1,0 


MPa. 


The propagation of the blast waves in two-phase gas-liquid medium with bubbly 
structure was analyzed in [6.73,6.79,6.90], where the decrease of the blast waves speed and 
their effective attenuation by the subsequent rarefaction waves is confirmed. The inclusion 
of the gas bubbles in liquid localizator results in the increase of the time of destruction and 
changes all time intervals of interaction. Additional (to a cycle of tests [6.57]) experiments 
were carried out with the purpose to define the consequences of replacement of a liquid in 


explosion localizator by the transformation of the blast wave. 


6.12. Experimental technique and verification 


The TNT charges of mass 75 g equipped by the electrical blasting caps ED-8 were 


used. 


Model samples simulating the dispersed liquid layer in an elastic envelope, were 
produced in the form of containers made from a goffered cardboard with the thickness of 3 
mm and the following sizes: the thickness (the distance between the gauges) d= 115 mm, 
the width / = 300 mm, the height H = 340 mm, the height of liquid filling level / = 300 
mm. The tensoresistor gauges were pasted inside the opposite walls in the middle part of 
the container. The hermetical package with liquid to be dispersed was placed into the 
container tightly to the walls. The two-phase gas-liquid medium of bubbly structure with 


gas volume fraction 50...60% was used as the working fluid. 


The choice of a working body for filling the explosion suppressor is based on a 
principal property of two-phase gas-liquid media. It is known [6.73], that the equilibrium 


sound speed in such medium is determined by: 


ce =P pa(i-a)| 


(6.24) 


The minimal sound speed is realized at the volume gas fraction in a liquid a= 0,5, 
when the product a(1-@) reaches the maximal value 0,25. In a water-air medium at the 
initial pressure Py = 0,1 MPa we have C;;,= 20 m/s. Thus, the alternation of sound speeds 
during the blast wave propagation from HE charge proceeds as follows: 340 m/s (air) > 20 


m/s (two-phase screen) > 340 m/s (air). 


The tests were carried out in the open district. Samples were placed on a steel plate. 
The charge was suspended on the rack so, that the center of a charge was on an axis of the 
gauges and perpendicular to the wall of the container. The distance from the center of a 


charge to the nearest gauge was 1 = 0,1 m, and to a surface of a plate - 0,17 m. 
The scheme of tests is presented in the fig. 6.36. 


To register the moment of the blast wave arrival the tensoresistor gauges were used. 
The gauges were made from electrically conducting paper, which resistance is reduced at 
compression. Registration and subsequent processing of the signals were made using 16- 


channel computer oscilloscope with the discrete frequency 3 MHz. 
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Fig. 6.36. The scheme of tests 


The records of the signals from the gauges allowed estimating the time intervals As 
between the impulses at the first and the second gauges (by the start and the amplitude of 
the pulse). The average by the container thickness speed of the blast wave was determined 


by the relationship: D=6/Ar. 


Experiments with the samples filled by water were carried out for the verification 
of the proposed measurement technique. The typical example of a record in coordinates 
output voltage of tensoresistor gauge U/ — time / is presented in the fig. 6.37. The signals of 
the channels 1,2 correspond to the windward (nearest to HE charge) and leeward (behind 


the screen) gauges and are given in a uniform time scale. 


The results of reference tests are presented in the first row of tab. 6.3 by the 


following parameters: 


tg) and tg. are the durations of the output signal rise at the near and far gauge along 


to the direction of blast wave propagation. 


The listed parameters are given with the indication of the maximal disorder of their 


values. As it is seen, the average speed of the blast wave is D,, = 1380 + 30 m/s. 


The maximal overpressure AP,,, in the blast wave at the forward to a charge border 
air - screen falls in a near zone of explosion [6.80]. At the chosen configuration charge - 
screen AP, ~ 10 MPa. The reflection of this wave results in transmitting wave into the 


screen with overpressure AP, ~ 10 MPa. 
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Fig. 6.37. The records of the output signals of tensoresistor gauges in the experiments with 
liquid to be dispersed without gas inclusions 


The wave speed calculated using Tait equation of state [6.73,6.76] in water is D, = 
1580 m/s. The difference between the calculated and the measured values D, and D,, is 
12,3 %, that it is possible to consider satisfactory for the accepted method of the 
measurement. There are also objective reasons of an expected divergence because of 
placing the water layer in the container with the flexible walls. Besides, in experiences 


with the water no special degasation of the layer was performed. 


Table 6.3 Results of measurements 


Type of working} ‘1, us ti2, LS D, m/s 
fluid 

Water 24413 15405 1380 + 30 
Gas - liquid 18409 56414 134 + 16 


6.13. Results of experiments with gas-liquid screens 


Consider the experiments on the blast wave transformation by a layer of two-phase 
gas-liquid working fluid. Typical records of the output signals U/ - 7 of tensor-resistance 


gauges are presented in the fig. 6.38. 
The basic parameters are listed in the second row of tab. 6.3. 


Comparison of the records of the gauges in a water layer (fig. 6.37) and in gas- 
liquid screen (fig. 6.38) shows the qualitative distinction of the blast wave transformation 
process. The duration of the pressure rise at the wave front in the water does not grow. An 
essential (in ten times) stretching of the blast wave is always marked in gas-liquid screen. 
The speed of the blast wave in two-phase medium is almost 10 times lower in comparison 


with the water. 
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Fig. 6.38. Experimental records for gas-liquid tests. 


The output signals from the nearest gauge to a charge (1) and the distant gauge 
behind the screen for the water (2) and for the gas-liquid medium (3) are combined in the 
fig. 6.39. In addition, it is possible to note, that the signal amplitude at the distant gauge in 
water is almost in 3 times higher, than in two-phase medium. It is easy to note quantitative 


and qualitative distinction of output signals using record, presented in the fig. 6.40. 


It means, that the average Mach number of a blast wave in a working body of 
explosion suppressor between windward and back walls is about => 6,8. The 
overpressure in this case is about P|/PoxM*>46. The earlier determined value of pressure in 
a wave penetrating into two-phase medium, was equal to P;*~ 800, that corresponds to the 
Mach number Mx28,2. Thus, at the distance of 0= 115 mm in two-phase medium with @ » 
0,5 the blast wave speed decrease almost in 4 times, and the front overpressure - not less 


than in 16 times. 


Fig. 6.39. Comparison of output signals U -t at the blast waves transformation in water 
and gas-liquid layer: At; is the time of blast wave crosses the water layer, Atz is the time of 
blast wave crosses the gas-liquid layer 
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Fig. 6.40. The combined signals from the second gauge for the samples with water (1) and 
gas-liquid medium (2): t* is the time after blast wave crosses protective layer 


As a result, the layer of two-phase medium located on the way of the blast wave, 
not only reduced the amplitude of the wave past through a layer, but also transformed the 
form of the blast loading. Instead of the blast wave with abrupt front ahead of the screen 


we obtained the low-frequency compression wave behind the screen, which is stretched in 


time. 


It is necessary to emphasize the universal character of the marked phenomenon 
taking into account the data of the experiments in [6.82]. In this work the research of the 
transformation of underwater blast loadings generated by the TNT charges immersed in the 
water was performed. The transformation of underwater blast loading was made by the 
installation of the porous polyurethane plate (the density of the material is 27... 64 kg/m’) 


of the thickness d= 1... 10 cm on the way of the blast wave propagation. 


The distance from the porous screen to the charge L = 30... 110 cm. In the 
experiments performed in [6.82] the sound speed on the way of the blast wave propagation 
alternated as follows: 1470 m/s (water) > 60... 100 m/s (polyurethane foam) + 1470 m/s 


(water). 


The results of the experiments on the transformation of consequences of underwater 
explosion [6.82] are qualitatively and quantitatively close to the described results of the 


experiments on the transformation of the air blast waves. 


The amplitude of the underwater blast wave in [6.82] at the passage of a layer from 
a material with the lowered speed of a sound is reduced in 30... 50 times. The time of the 
pressure rise at the front waves has increased almost in 40 times and the blast loading was 
transformed into the compression wave, stretched in time. The similarity of suppression 
results for the air and the underwater explosions allows to assume an opportunity of 
creation of the universal accident protection devices based on the principles of the air blast 


waves attenuation, formulated in [6.67... 6.71]. 
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